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: A Note on the Position of the Centre of Flexure 
of a Beam of Isotropic Material Having a Section Bounded 
by a Parabola and a Straight Line 


By Amitava CHAKRAvortI, Calcutta, India 1) 


Nomenclature 
ax, V, Z rectangular co-ordinates; 
_lm,n direction cosines of outward drawn normal; 


 9,,%,,0, normal components of stress parallel to x, y, z axes; 
Try: Tyz Tzx Shearing stress components in rectangular co-ordinates; 
u,v, w components of displacement; 

_ T(x, y) stress function in rectangular coordinates; 

6 Polsson’s ratio; 

77 shear modulus; 

4 Bias Yce Ve, Shearing strain components in rectangular coordinates. 


i _Methadvot Solution 


Let a beam of uniform cross section with a straight edge have one end z = 0 
fixed and the other end z= 1 loaded by some distribution of forces that is 
_ statically equivalent to a single force W acting parallel to the straight edge. 
- At the end z = 0, axis of x is taken parallel to the load, y axis is taken perpen- 
dicular to the straight edge and z axis is taken along central line of the beam. 
Here x and y axes are parallel to principal axes of inertia at the centroid. The 

origin is at the middle point of straight edge. aie 
_ Following the semi-inverse method of DE ST. VENANT, we assume 


Ww 
Ce das t= 0) “and 6,=—= U—2)%, (1.1) 
where J = ‘i i. x dx dy over entire cross section. 
From equations of equilibrium we get | 
OT, = Cte =H (1.2) | : : rh. 


(ay a Oz 


= 1) Department of Mathematics, Brahmananda Keshab Chandra College. 
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and 


xz Oty, Aa Sir 1.3 
ox | 0 a al Os ui 


Equations (1.2) show that t,, and t,, are functions of x and y only. Equa- 
tion (1.3) is satisfied if we assume 


oT Ww oT 
Taz = > + Sly) — ee x, Tey BER ines (1.4) 


where S(y) is a function of y to be so chosen as to yield a simple boundary 
condition. Beltrami-Michell compatibility equations reduce in this case to 


il 0*o, =0;. V2 t%,=0. 


2 p , 
Ve tee + Leo 8 600x 02 


Therefore from (1.4) 


Hence, we have 


V#TR 9) = ee (1.5) 


where yp is y co-ordinate of centre of inertia of cross section. Local twist at 
point (x, y) of cross section is 


Ow 1 0 [ov Ou 1 (i) 0 

Gee 1t24" 0s tae ea) tae i Yue ~ Fy Vex) 
=a, (92-42) --a[r76+3] by aan 
=~ $7 [5- Fo) -2nal. iby (1.5)! 


Therefore mean local twist equals «. 

We define centre of flexure as position of load point for which the mean 
local twist « is equal to zero. Then we can take « = 0 in (1.5). Now let (z, y, ) 
be co-ordinates of centre of flexure. In our case since the load W acts parallel 
to a principal axis (x axis), hence x co-ordinate of centre of flexure is zero. 
y = y co-ordinate of centre of flexure is given by 


—y We |] ty, — 9 ty) de dy. 
Rk 
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Let C denote boundary of cross section and R, the region within C. Then from 
(1.4) 


2 y pe owe oT 4 ‘ 
—yW =f (= ae ) —¥ 35 + S(y) — hal xl dx dy | 
=I [2 T(x, v) + rs {= x T(x, v) — xv S(y) +2 x y} 
4, Ly LK: y) tI dx dy (1.6) 


an . 


I 


Z I| T(x, v) dx dy + | y T(x, y) dx 


R é 


( = W 
+ | [- x T(x, y) —xv S(y) + eT xP y| dy 
z 
(by GREEN’s theorem). 


Again we have the boundary condition 


i, + 2 t,=0, 


The function S(y) may be prescribed arbitrarily. We choose S(y) = 42 W/2 I 
on C so that above condition becomes dT /ds = 0. We take T = 0 on C since the 
constant term does not affect the stresses. Then we see that T(x, y) satisfies the 
differential equation 


o W dS 
VY? T(x, y) = Eiage ie (P= Vg) ss dy. (T = 0 on C) 
and S(y) = «2 W/2 I on C except that S(y) may take any value where dy/ds =0. 
Then from above considerations y is simplified and we get from (1.6) 
= 2 oe 1 F 
Af oe || IL ta, yvyedx dy + ees xy ay . CLN7) 
4 c 


ra 


R 


2. Section Bounded by a Parabola and a Straight Edge 


Here we assume (x + 7 y) = —7 (&+77)?. 


Therefore 
ea2eng, y= 17 —e. (2.1) 
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Then 7 = + a gives the parabolic curved boundary 
a 
a x? = 4 a® (a? — y) . (2.2)9 
ry ee Sk = The straight boundary is y = O or & = 27. 
In this case 
i. 1 Ox\2 (dy . 
oy = (2) + (GE) 48 +0), 
et a Ga 3 
x 2a 32 
ores creas ae he é 


In this case on parabolic boundary C 


W W 
Sy) aT = age 4 a? (a2 — y). 


Therefore T(x, y) satisfies the differential equation 


: W 2 W-a2 
YT, y= 4 : =o Re a Nal sie a 
2.4 
ers 2Wat (54+ 40) 2) 
“To fo oe Soe 


T =0 ony =a and also ony = 0, i.e. & = + 7. (2.4) is satisfied if we assume 
a particular solution of T, 


T=Ay(Pt+4a@y+By?+Cy, (2.5) 
where 
oW W a? 
A Feat to sa oe ee ee 


I) = 0ony=0. But T= 0 ony = + a also. 
Therefore we assume 


oe On 
m=0 
eS (2m+1)x7 (2m+1)aé 2m+1 é 2 ut 
x [cosh ie, a Se ae wee ee toe mat) 


This assumption of T satisfies differential equation (2.4) and gives T = 0 
ony=0,ie.€=+47. 
Since T = 0 on 7 = a hence it is necessary that 


O=Bé&—&(4Aat+ 2Bat+C)+a(4atdAt+ Bar+C) 


=> Cy gan COSk ek Pe cos ee (—a<&<a). 


m=0 2a 


“sen 


Te eT ee ae? a aT eT eee eT eT Ve 
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Putting 
Wai (54+ 40) 
he ay GY he peeled Te Oa pat 
C 4a4A—Ba s7 ae ; (2.8) 
we see that the above condition reduces to 


ja) 


Be Bist — a 2) +) Cag , cosh- Og A Rae ER (=e F< a). 


2 2 
m=0 a 


This relation gives 


C 32 Bat 48 sin (2m + 1) 2 
PENG aa oe aft a )3 5 3 — Fe TSE — ~ es 
(2m-+1)§a (2m+ 1)? Sieh aed En ' (2.9) 


(ree Dn Le 2eiey ale 


where B is given by (2.6). Thus 7 is obtained. Now from (1.7) substituting the 
value of T we get coordinates of centre of flexure. 


Here 
— a? (25 + 124) 512a4* 
y 35 (1 +o) a3 W p> Comes 

a Sa a ot 

sin (2m + 1) —-cosh(2 m + 1) = 3 sin(2m + 1) > sinh (2 m+ 1) > 

* (2m + 1)3 a? se a (2m-+ 1)4 j 
Substituting values of Cz,,., from (2.9) we get 
Eo a* (25.4 12 0) glee Genta) .5" | Oe ee | 
y 35(1 +0) a (ito) “|(2m+i1je  n®(2m-+ 1/8 


ieebieaaee 48 tanh (2 m+ 02 | 
bs 


4608 a? (5 — 60 ee 
sy (2m-- 1)? mot (2 


x? (1 + 0) 


Next using the results [2] ?) 


a 1 a8 ee 1 tb esas 
2 (2m+1)® — 960" Pe (2m+1je ~~ 2-7! 
we get 
3 (126 — 3) a? 4608 a* (5 — 6 a) a? | 
2 ae Bei ce Sr ae 
(1 + 9) mi ( ) (2.10) 
48 
“ee h(2 1 
oe atom py) anhiam + 1) > | 


The infinite series occuring in ohh are rapidly convergent. 


2) Numbers in brackets refer to References, page 338. 
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Let us consider those two series. 


co 1 Tu n [e.e) os we a 
oe (am.+ 1) tanh (2m + 1) = tanh> +d em ry tanh (2m +1)>. 


= 1 
tanh> > 0-9 and AS mate tanh (2 m + 1) 
m=1 


E14 
2 


is smaller than 


m=1 


approximately. Hence for practical purposes we may take only first term of 
series 


S (@m pay tanh@m+ 1) y. 
Similarly from the series 
uy ( mE tanh (2 m + 1) > 
we take first term. Then 
or F eA ent 2: (5 “a tanh > (approximately) : 
re on ah ae 2 a - RSs (approximately) . 


If we take o = 0-3 then y = 0-536 a? (approximately). It may be observed that 
x co-ordinate of centre of flexure is zero. 


In conclusion the author expresses his deep gratitude to Prof. B. SEN, 
D. Sc., F. N. I. for his guidance and helpful suggestions. 
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Zusammenfassung 


Die vorliegende Notiz behandelt das Problem der Auffindung des Biegungs- 
zentrums eines aus isotropem Material bestehenden Balkens, der einen durch eine 
parabolische und eine geradlinige Kante begrenzten Querschnitt besitzt. Die Last 
ist parallel zur geradlinigen Kante angenommen. Der Verfasser beniitzt zur Lésung 
der Aufgabe dasselbe Verfahren, das SoKoLNnikorF [1] fiir die Auffindung des 
Biegungszentrums eines halbkreisférmigen Balkens verwendete. 


(Received: November 13, 1958.) 
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Der Analog-Interpolator 


Von THEODOR ERISMANN, ETH, Ziirich 


Problemstellung 


Das ungeheure Zahlenmaterial, das moderne programmgesteuerte Rechen- 
automaten zu liefern vermégen, macht die Auswertung der gerechneten Resul- 
tate gelegentlich zu einem schwierigen Problem. Dies gilt besonders dort, wo 
die tibliche Ausgabe in Form von Zahlentabellen (allenfalls auf Lochkarten oder 
ahnlichen Tragern gespeichert) ungeeignet ist, weil eine anschauliche Darstel- 
lung des Resultats oder dessen Weiterverarbeitung in Analogform (etwa zur 
Steuerung von Werkzeugmaschinen) gewiinscht wird. 

Die bekannten elektrischen Digital-Analog-Umformer sind zwar zur Ver- 
anschaulichung der Rechenresultate durchaus geeignet, doch scheiden sie iiber- 
all dort aus, wo hohe Genauigkeiten erforderlich sind (Kurvenscharen in Rechen- 
blattern, Frasen von Schaufelprofilen usw.). Zudem erfordert ihre diskrete 
Arbeitsweise die Berechnung einer sehr grossen Anzahl von Punkten, sofern 
man nicht nachtraglich noch zum Kurvenlineal oder zu ahnlichen behelfsmas- 
sigen Interpolationsmitteln greifen will. 

Ein Gerat, das die geschilderten Aufgaben zu lésen hat, sollte also die fol- 
genden Eigenschaften in sich vereinigen: 

1. Umformung der verarbeiteten Zahlen aus der Codesprache aes Rechen- 
automaten in eine stetig sich 4ndernde Analogiegrésse. 

2. Méglichst geringer Genauigkeitsverlust der Analogdarstellung gegeniiber 
der Digitaldarstellung. 

3. Méglichst geringer Zeitverlust zwischen der Eingabe in Digital- und der 
Ausgabe in Analogform. 

Ein solches Gerat wird in der Folge als Analog-Interpolator bezeichnet. 


Das Interpolationsverfahren 


Der Ubergang von der diskreten zur stetigen Darstellung erfordert ein 
Interpolationsverfahren, das nach dem oben Gesagten folgenden Bedingungen 


geniigen muss: 
1. Die interpolierte Funktion muss mindestens mit der durch ihre weitere 


Verwendung bedingten Genauigkeit durch eine beliebige Anzahl von Punkten 
gelegt werden k6nnen. 
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2. Im Hinblick auf die praktische Verwendung (Oberflachenglatte usw.) 
muss die interpolierte Funktion stetig sein. Von der zweiten Ableitung an sind 
Spriinge normalerweise zulassig. 

3. Die interpolierte Funktion muss bei laufender Eingabe der Koordinaten 
fiir die aufeinanderfolgenden Punkte laufend gebildet werden. Zur Bestimmung 
ihres Verlaufes zwischen den benachbarten Punkten P;_, und P, sollte nach 
Moglichkeit die Kenntnis der Koordinaten vom Anfangspunkt P, bis P, ge- 
nugen. 

Auf den ersten Blick scheint eine aus quadratischen Parabelbogen zusam- 
mengesetzte Funktion den erwahnten Bedingungen bei bescheidenem techni- 
schem Aufwand zu geniigen: Ist der Verlauf der Funktion bis zum Punkt P;_; 
gegeben, so gibt es eine und nur eine Parabel, die in P;_, beztiglich des Funk- 
tionswertes und der ersten Ableitung mit dem vorhergehenden Teilbogen tiber- 
einstimmt und durch den nachstfolgenden Punkt P; geht. Einer Fortsetzung 
des Verfahrens ad libitum steht also nichts im Wege. 

Leider zeigt sich aber bei naherer Betrachtung, dass die beschriebene 
Parabelbogen-Interpolation eine gewisse Verwandtschaft mit dem Verhalten 
eines ungedampft schwingungsfahigen Gebildes aufweist: Liegen beispielsweise 
in einem gewissen Bereich die Stiitzpunkte der Funktion aquidistant auf der 
Abszissenachse, so geht die interpolierte Funktion nicht etwa nach einigen” 
Teilbogen in diese Achse iiber, sondern pendelt — wie man sich leicht iberzeugen 
kann — in gleichmassig wiederholten Parabelbogen beliebig haufig weiter um 
dieselbe. Es ist klar, dass ein solches Verfahren, dessen periodische Fehler 
offenbar sowohl von den Anfangsbedingungen als auch vom Verlauf der inter- 
polierten Punktfolge abhangen, in der Praxis genauigkeitsmassig unkontrollier- 
bar (und unter Umstanden sogar instabil) ist und dass nur bei Einfiithrung 
einer wirksamen «Dampfung» Aussicht auf Erfolg besteht. 

Der Grundgedanke einer solchen Dampfung geht aus Figur 1 hervor, welche 
die Umgebung eines Teilbogens darstellt. Hier wird die Parabel nicht durch die 
Punkte P;_, und P,, sondern durch zwei um den kleinen konstanten Betrag e 
nach links verschobene Punkte Q;_, und Q; gelegt. Zwar wird von Q;_, aus der 
durch P; gehende Parabelbogen angesetzt, doch wird dieser nur bis Q; gefiihrt, 
von wo der nachste Teilbogen nach P,,, bzw. Q;,, begonnen wird. 

Unter Verwendung des aus Figur 1 ersichtlichen Koordinatensystems é, 
und der ebendort abzulesenden Bezeichnungen Ax und yu erhalt man fiir den 
Parabelbogen (gegeben durch den Funktionswert Yoi—11m Punkt Q,_,, die nach 


5 baw. x gebildete Ableitung yp; _, im gleichen Punkt und den Funktionswert 
yp; 1m Punkt P,) die Gleichung 


= ; ae. pe { Ups ued vOu=4 
y YOi=1 4 E Voit oi | a ae a tee (1) 


Damit ist der Verlauf dieses Teilbogens vollstandig umschrieben. Insbesondere 


SVX 
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ist die Bestimmung der massgebenden Werte im Punkt Q; moglich, die sich mit 


we —1 ’ ie apy 
.. == VY. P Al . 1 =! i } 
“Qi “Qi-1 be r YOi—1? I mn F be (2) 
und 
pha ie o uw 2 2 
/Qi /Qi-—1 p® Ax bi YOi-1 ub 7 YPi wu? Ax (3) 


ergeben und als Ausgangsbasis fiir die Ermittlung des nachsten Teilbogens 
dienen. 


Figur 1 


Prinzip der gedampften parabolischen Interpolation. Erlauterungen im Text. 


Trotz ihrer Einfachheit sind diese Formeln fiir eine umfassende theoretische 
Analyse des Verfahrens iiber mehrere Teilbogen dusserst umstandlich, weil sie 
beim Ubergang von einem Teilbogen zum nachsten zu immer laéngeren und 
undurchsichtigeren Ausdriicken fiithren. Es erweist sich daher als zweckmassi- 
ger, eine Anzahl typischer konkreter Falle durchzurechnen. Wird dabei syste- 
matisch vorgegangen, so kénnen die Ergebnisse hinsichtlich ihres praktischen 
Wertes denen einer rein theoretischen Analyse gleichgestellt werden. 

Zunachst wird der bereits erwahnte Fall einer aquidistanten Punktfolge auf 
der Abszissenachse betrachtet. Dabei geniigt fiir einen vollstandigen Uberblick 
die Variation der drei Parameter , yy und yo. Wie aus Figur 2 ersichtlich, 
spielt yw die Rolle eines Dampfungskoeffizienten, der selbst bei ungtinstigsten 
Anfangsbedingungen ein rasches «Einschwingen» auf den Sollwert gestattet 
und etwa von yw = 1,5 an einen praktisch aperiodischen Verlauf sicherstellt. 
Damit sind alle Voraussetzungen fiir eine vollstandige Beherrschung der Ver- 
haltnisse gegeben. 

Die Durchleuchtung des vorliegenden Falles ist trotz seiner Einfachheit sehr 
wertvoll, da der Giiltigkeitsbereich der gewonnenen Resultate durch die fol- 


ZAMP X/22 
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genden Betrachtungen ganz wesentlich erweitert werden kann: Liegt die Punkt- 
folge nicht auf der Abszissenachse, sondern auf einer beliebigen quadratischen 
Parabel, so lasst sich leicht zeigen, dass die Differenzen zwischen den Ordinaten 
der interpolierten Funktion und der durch die Punktfolge gehenden Parabel — 


Figur 2 


Einschwingen aus verschiedenen charakteristischen Anfangsbedingungen fiir verschiedene Werte 
des Dampfungskoeffizienten ju. 


abgesehen von einer sinngemassen Transformation der Anfangsbedingungen — 
identisch sind mit den Ordinaten der interpolierten Funktion im oben behan- 
delten Fall. Mit anderen Worten: Das rasche Einschwingen auf vernachlassig- 
bare Fehlerwerte und die wirksame Dampfung bleiben auch dann in vollem— 
Umfang erhalten, wenn die Punktfolge in einem itber einige Teilstiicke gehenden | 
Bereich auf einer quadratischen Parabel liegt. 

Damit ist das Verhalten des Verfahrens fiir alle diejenigen Falle abgeklart, 
in denen die Fehler nur durch die Anfangsbedingungen entstehen, wahrend 
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bei richtiger Wahl dieser letzteren die interpotierte Funktion vollstandig mit 
dem durch die Punktfolge gehenden héchstens quadratischen Polynom iiber- 
einstimmt. Dagegen bleibt die Frage nach der Giite des Verfahrens bei der 
Interpolation beliebiger anderer Funktionen noch offen. Fiir die praktische 
Beurteilung geniigt dabei die Untersuchung des Falles einer kubischen Parabel. 
In der Tat lasst sich jede beliebige Funktion in einem gewissen Bereich stets 
um eine Gréssenordnung besser durch eine kubische als durch eine quadra- 
tische Parabel approximieren. 

Die etwas umstandliche Berechnung der Fehler bei der Interpolation einer 
kubischen Parabel kann im Rahmen dieser Abhandlung nicht in extenso be- 
handelt werden. Dagegen seien hier die Ergebnisse einer solchen Untersuchung 
in knapper Form zusammengestellt : 

1. Folgt die zu interpolierende Funktion tiber mehrere Teilstiicke einer 
kubischen Parabel, so streben die Fehler in den Punkten Q; rasch gegen einen 
konstanten Wert Ay,,, der durch die Beziehung 


m A’; 1 
Ay. = ¥ =e Blu — 5) (4) 


gegeben ist. Diese Fehler haben somit systematischen Charakter, was die Be- 
zeichnung «systematische kubische Fehler» rechtfertigt. 

2. Innerhalb jedes Teilbogens iiberlagert sich dem systematischen kubischen 
Fehler ein weiterer Fehler dy, der ebenfalls rasch einem fiir jeden Teilbogen 
gleichférmigen Verlauf zustrebt und daher als «periodischer kubischer Fehler» 
bezeichnet werden darf. Dieser Fehler verschwindet fiir € = 0 und € = Ax und 
erreicht zwei Extremwerte etwa bei = 0,21 Ax und & = 0,79 Ax. Der Absolut- 
wert dieses Fehlers an den bezeichneten Stellen kann durch die Formel 


| 6.5 max| Gz 0,008 y” ASx (5) 
ausgedriickt werden. 
Es iiberrascht keineswegs, dass beide Fehler — der systematische wie der 
periodische — der dritten Ableitung der interpolierten Funktion proportional 
sind; ist doch eben das Vorhandensein dieser Ableitung das Unterscheidungs- 
kriterium zwischen der quadratischen und der kubischen Parabel. Dagegen ist 
es bemerkenswert, dass der periodische Fehler vom «Dampfungsfaktor» « un- 
abhangig ist. Im iibrigen ist die Amplitude des periodischen Fehlers fiir die in 
Frage kommenden Werte von mu in der Gréssenordnung von 1,5 stets betracht- 
lich kleiner als der systematische Fehler. Das Pendeln um diesen letzteren dart 
als durchaus harmlos bezeichnet werden (fiir w= 1,5 ist My,, = 15,6 bOyepinneall 
Sehr wesentlich ist schliesslich die Abhangigkeit beider Fehler von der dritten 
Potenz der Schrittlinge Ax, die — bei wenigstens approximativer Kenntnis der 
dritten Ableitung von y — ein einfaches Mittel zur Beherrschung der Fehler- 
grossen darstellt. 
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Im Rahmen der Untersuchung von Punktfolgen, die in x aquidistant sind, 
bleibt nur noch die Frage offen, wie sich die Interpolation von Funktionen 
gestaltet, deren dritte Ableitung nicht konstant ist, und zwar insbesondere dann, 
wenn fiir die Wahl der Anfangsbedingungen nur die Kenntnis des Funktions- 
wertes und der ersten Ableitung im Anfangspunkt vorausgesetzt wird. Aus 
einer grésseren Anzahl von Funktionen, die zur Abklarung dieser Frage unter- 
sucht wurden, sei als besonders ungiinstiger Fall die Hyperbel y = 1/% im 
Bereich 1 < x < 2, also unweit der Singularitat in x = 0, herausgegriffen. Wie — 


LL 


Figur 3 
Fehlerverlauf bei Interpolation einer auf der Hyperbel 1/x liegenden Punktfolge. Schrittgrdsse 
Ax = 0,1, Dampfungskoeffizient w = 1,5. Anfangsbedingungen durch Funktionswert und erste 
Ableitung im Punkt *« = 0,9 gegeben. 


Figur 3 zeigt, erreicht auch hier der Fehler sehr rasch Werte, die annahernd im 
Rahmen der theoretisch vorausbestimmten Grenzen Ay, + | bo max| (gestri- 
chelte Kurven in der Abbildung) liegen. Ubrigens wird man beachten, dass die 
Fehler trotz der verhaltnismassig grossen Schrittlange Ax = 0,1 bereits ausser- 
ordentlich klein sind. 

Damit darf festgestellt werden, dass das Verfahren — wenigstens fiir in « 
dquidistante Punktfolgen — den eingangs formulierten Bedingungen entspricht 
und insbesondere folgende Eigenschaften besitzt: 

1. Weitgehende Unempfindlichkeit gegeniiber ungiinstiger Wahl der An- 
fangsbedingungen. 

2. Verhaltnismassig leichte Bestimmbarkeit der zu erwartenden Fehler. 

3. Einfache Mittel zur Reduktion der Fehler auf beliebig kleine Werte 
(durch Reduktion der Schrittlange). 

4. Moglichkeit zur Glattung zufalliger Unregelmissigkeiten (Streuungen) 
der Punktfolge (durch Wahl des Dampfungsfaktors ju). 

Ist die Punktfolge in x nicht dquidistant, so sind keine wesentlichen Uber- 
raschungen zu erwarten, wenn die Schrittlinge sich nach einem stetigen Gesetz 


wat 
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andert. Immerhin ist bei der Fehlerabschatzung einige Vorsicht am Platze, da 
die oben dargelegten Beziechungen nur unter gewissen Umstinden in vollem 
Umfang anwendbar sind. In der Praxis wird man bei zweidimensionalen Pro- 
blemen allerdings sehr oft eine Parameterdarstellung x(i), y(i) wahlen, in der 
die Ordnungszahl 7 als gemeinsame unabhangige Veranderliche auftritt und 
die Schrittlange mit 47 = 1 konstant ist. Diese Darstellung gestattet zudem die 
Beherrschung gewisser Falle, die der expliziten Form (x) bei parabolischer 
Interpolation nicht zuganglich sind (Durchgang durch vertikale Tangenten, 
geschlossene Kurven usw.). 


Die technische Verwirklichung 


Das Geradt muss offenbar aus zwei Hauptteilen bestehen, namlich einem 
Digital-Analog-Umformer und dem eigentlichen Interpolator, an den dann die 
entsprechenden Registrier- oder Arbeitsgerate angeschlossen werden. 

Die technischen Mittel der Digital-Analog-Umformung sind bekannt und 
bediirfen keiner weiteren Erklarungen. Als digitalen Eingang wird man nor- 
malerweise nicht den Rechenautomaten selbst, sondern einen geeigneten Zwi- 
schentrager (Lochkarte, Magnetband usw.), bei einfacheren Geraten sogar nur 
eine Tastatur bentitzen. Da es sich meist um die Steuerung mechanischer Vor- 
richtungen (Registriergerate, Werkzeugmaschinen usw.) handeln und in vielen 
Fallen eine sehr hohe Genauigkeit verlangt werden diirfte, kommt als Analogie- 
grosse wohl nur der Umdrehungswinkel ciner Welle in Betracht, der als einzige 
leicht verarbeitbare Analogiegrésse durch Wahl eines geniigend grossen Uber- 
setzungsverhaltnisses eine beliebige Genauigkeit auf beschranktem Raum ge- 
stattet. 

Der Aufbau des eigentlichen Jnterpolators ergibt sich in einfachster Weise 
durch Anwendung der Schaltungstechnik fiir Integrieranlagen: Durch zwei- 
malige Differentiation und anschliessende zweimalige Integration der Bezie- 
hung (1) erhalt man 


ep veer vogel. pemeein) \ ton 
y=2i]( p? A®x TP Ae ahs 


ve 


woraus ersichtlich ist, dass zur Darstellung von y zwei Integratoren 3 und 4 in 
der Schaltung gemass Figur 4 erforderlich sind. Das Schema bedarf nur weniger 
Bemerkungen zu seiner Erlauterung: Wahrend der Berechnung eines Teil- 
bogens ist y” konstant, und die Integratoren rechnen daraus bei kontinuier- 
licher Eingabe von x laufend die interpolierten y-Werte. Die Verbindung 
zwischen dem Summentrieb 2 und dem Integrator 3 ist wahrend dieses Vor- 
ganges unterbrochen (gestrichelte Linie). In einem beliebigen Zeitpunkt wird 
der nachste Wert fiir yp am Summentrieb 7 eingegeben, so dass beim Erreichen 


des nachsten Punktes der Q-Reihe am Ausgang des Summentriebes 2 der neue 


Wert fiir y” bereitsteht. In diesem Augenblick wird der Antrieb nach x unter- 
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brochen und durch eine (gemass der strichpunktierten Linie gemeinsam mit 
x gesteuerte) geeignete Nachlaufsteuerung die Verbindung zum Integrator 3 
hergestellt. Damit ist das Gerat zur Berechnung des nachsten Teilbogens bereit. 
Die Getriebe 5 und 6 dienen zur Einfiihrung der Konstanten w und Ax. 


Figur 4 


Interpolationsschaltung mit stufenlosen Getrieben (Integratoren). 


Am Rande sei erwahnt, dass das System von Teilbogen zu Teilbogen selbst- 
korrigierend wirkt, weshalb die Rechenelemente nur bescheidenen Genauig- 
keitsforderungen zu gentigen brauchen. Ferner ist zu bemerken, dass bei genii- 
gend rascher Arbeit der Nachlaufsteuerung fiir y” gegebenenfalls mit konti- 
nuierlich fortlaufendem Antrieb nach x gearbeitet werden kann. 


Anwendungsm6$lichkeiten 


Der Analoginterpolator kann in erster Linie zwei Aufgaben erfiillen: Die 
Veranschaulichung der Resultate digitaler Berechnungen und die Automati- 
sierung gewisser Produktionsvorgange. 

Im ersten Fall wird man das Gerat direkt mit einem geeigneten Registrier- 
apparat koppeln, wobei gegeniiber anderen Verfahren vor allem die hohe 
Genauigkeit ausgeniitzt werden diirfte, die es gestattet, hochpriazise Kurven- 
blatter fiir nomographische und ahnliche Rechen- und Auswertungsarbeiten 
herzustellen. 

Noch niitzlicher diirfte der Einsatz des Analoginterpolators fiir die Auto- 
matisierung von Produktionsvorgangen sein. Insbesondere bei der Herstellung 
kleiner Stiickzahlen (Prototypen und Nullserien) von Leitkurven, Nocken, 
Schaufelprofilen und ahnlichen Kérpern, wo die Anfertigung von Modellen fiir 


ti. 2 


oT. eer. ea Ce 
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Kopierfrasmaschinen sich nicht lohnt, steht ein wesentliches Betatigungsfeld 
offen. Dabei kann der Interpolationsvorgang entweder einmal durchgefiihrt und 
das Resultat in kontinuierlicher Weise (etwa mittels hochpraziser polydromer 
Funktionsdarstellung) auf einem-analog arbeitenden Zwischentriiger registriert 
werden, worauf es bei der Herstellung jedes Stiickes durch eine mit der Werk- 
zeugmaschinensteuerung verbundene geeignete Abtastvorrichtung abgetastet 
wird. Oder aber der Analoginterpolator wird selber zu einem Teil der Werk- 
zeugmaschinensteuerung und fiihrt die Interpolation — gesteuert von einem 
digitalen Zwischentrager (etwa einem Stoss Lochkarten) — bei der Herstellung 
jedes Stiickes von neuem durch. Angesichts der zahlreichen Argumente, die 
sich zugunsten beider Lésungen anfiihren lassen, wird die Praxis zeigen 
miissen, welche die grésseren Aussichten auf einen dauernden Erfolg besitzt. 


Summary 


The paper deals with a new method for transforming digital data tables into 
continuous shaft rotations. The method is characterized by a special type of 
parabolic interpolation incorporated into the analog part of the machine and 
conceived in a way to allow for perfect stability and smoothness of output. 
Formulae are given for predetermining the errors of the system for interpolating 
any desired function. 


(Eingegangen: 8. Dezember 1958.) 


The Dynamics of Turbo-Flow 
By Fausto G. Gravatos, Philadelphia, USA?) 


Foreword 


This work is essentially a dynamic analysis of the steady, axially symmetric 
flow of an inviscid, incompressible fluid through a turbo-machine. 

From the point of view of dynamics, the condition of axial symmetry makes 
it necessary to assume a field of forces physically equivalent to the presence 
of the blade row. This concept of force field appears — through rather vaguely — 
first in LoRENz’ [1]?). Von MisEs [2] defined the axial symmetric turbo-flow 
model by means of purely kinematic concepts; his equation has been quoted 
often in recent literature, but only MEYER [3] and MarBLE [4] show a true 
understanding of its meaning. 


1) Consulting Engineer, MSVD General Electric Company, Philadelphia, Pa., USA. (This 
work, finished in August 1953, was done while the author was at the Aircraft Gas Turbine Division 
of the General Electric Company, Cincinnati, Ohio.) 

2) Numbers in brackets refer to References, page 380. 
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From a physical point of view, the turbo-flow problem — under the above 
hypothesis — is reduced to the proper definition of the force field, F, so as to 
represent a blade row. 

From an analytical standpoint, the problem is to set down the conditions of 
compatibility between the equation of continuity and the three Euler equations. 

In the last years many attempts have been made to solve this problem. 
The writer showed [5] that if the component of F along the normal to the 
streamlines is zero, the problem reduces to the solution of STOKES’ stream 
function and the integration of aset of ordinary differential equations. Naturally, 
this hypothesis cannot represent the general case. 

What has been missing is a clear-cut formulation of the problem. This is 
accomplished in part 4; to do so, it is necessary to analyze the geometric 
structure and the dynamic characteristics of potential force flows. This is done 
in part 3. From this study, there result two possible generalizations of potential 
force flows: 

(1) von Mises flows; 
(2) quasi-potential flows. 


This classification results from the manner in which a ‘blade’ is defined. 
When the flow is Lorenz —i. e. when the constant in EULER’s integral does not 
vary with wy — both types of flows coincide. Therefore, the solution of quasi- — 
potential force fields given here implies the solution of von MIsEs’ equation. 

The study of flows with the initial assumption that F be equal to 
a(r, z) [grad 6 — f(r, z)] appears in GOLDSTEIN [6], who deals with the direct — 
problem for radial blades and in MONROE’s [7] thesis. These authors are, how- 
ever, not concerned with the problem as stated here. 

Parts 1 and 2 are of an introductory nature. Part 2, however, is essential 
to the understanding of the vorticity field. 

In part 5, where an illustrative example is worked out, the important role 
that numerical analysis plays in modern technology is brought to the surface. 


Nomenclature 
Symbol Dimensions 
a a constant Tee 
B blade surface 0 
Cc velocity vector Die 
6 speed vl ae 
ie average axial velocity ya he 
d diameter EL 
F force per unit mass y Bt Ba 
i blade shape function 0 


pres wm. BPMN 
: 
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Symbol Dimensions 
g angular momentum per unit mass (g = C ery LATA 
h blade height rs 
Betas | 
2 @ 
BERNOULLI’s constant (stator) E2T= 
ues =5 04 [Pw 
¢ 
BERNOULLI’s constant (rotor) irs 
k a constant lel 
K a parameter Leer 
p pressure ME T=" 
Yr radius, distance from axis Ib, 
Ry outer radius E 
R; hub or inner radius EE 
s arc length 1s 
t time a 
u a unit vector depending on use 
Ww relative velocity ET 
a load function zp 
‘g family of surfaces 0 
6. a parameter 0 
a v-component of vorticity ee 
n 6-component of vorticity ND 
60 angle or direction 0 
x a parameter 0 
gE z-component of vorticity (he 
% the number 0 
0 density. Ny ae 
@ a function id Nae 
QP direction or angular coordinate in rotors 0 
y stream function Ea 
@ angular velocity po 
2 vorticity vector i 
Q vorticity cp 


The subscripts on C, W, and F indicate components of C, W, and F along 
the corresponding directions. Use has been made of C,= W,, C, = W,. Sub- 
scripts on w indicate the direction of the unit vector. Subscripts on any other 
letter indicate the derivative with respect to the subscript, with the exception 


350 Fausto G. GRAVALOS ZAMP 


of the subscripts 0,7, 1, and 2, the meaning of which is specified in the text. 
The total time derivative, d/dt is sometimes referred to as material derivative. 
This unnecessary terminology has been avoided. 


1. Geometric Preliminaries 
1.1 Cylindrical Coordinates 


A stationary — 7, 0, z - and a moving — 7, @, 2 — system of cylindrical coor- 
dinates will be used. The space orientation is defined by the cyclic order 
rv, 0,z—or r, y, —so that from 7 the positive sense of rotation is from 6 to z; 
from 0, the positive sense of rotation is from z to 7, etc. 

These two systems of coordinates are related by 


6=m+o@t, mw =angular velocity ss oh 


and the homologous components of any vector A in space are the same in both 
systems. 

With this convention, the components of 2 = curl C are given, in the 
axially symmetric case, by 


1 

E= es 82> (1.2) 

‘i. OC ne Oe 

ieee peer Fo (1.3) 
1 

a (1.4) 


Since the choice of sign in the definition of the stream function 


dl 

Ch => We» (1.5) 
1 

CL=—F Vr (1.6) 


is essentially arbitrary and it determines the sign in the Stokes’ equation 
1 
Pre — Y Yr Vez = VY (1.7) 


obtained by expressing 7 in terms of y, it will be shown that this arbitrariness 
is only apparent. 

; At the outset the sense of flow, a physical data, from left to right or from 
right to left, determines uniquely the sign of the stream function yp. If the 
system is clockwise, yp is positive (Figure la). For the other case y is negative 


—— 


he ee 


7 ew Nt 
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(Figure 10). In both cases, the sense of increasing y, given by w,,, is obtained by 
a negative 2/2 rotation of the vector Cy, i. e., 


C C 
Fel tens ) r 
: ( > «9, ek (1.8) 


Figure 1 


(a) Clockwise; (b) counterclockwise. 


This can be proven easily. If ds, indicates the lineal element along the 
u,, direction, the differential of y along uw, is given by 


d, y = [y, cos (t,,, r) + w, COS (u,,, z)| ds, 


and use of (1.5), (1.6) and (1.7) yields, 
grad, p=7Cy (1.9) 


an essentially positive quantity. 
The positive sense corresponding to 7 in both cases, (a) and (0), is shown 
in Figure 2, and it is seen that the choice of orientation in space and the selection 


Positive 1 Positive 7 


a b 
Figure 2 


of sign in (1.5) and (1.6) are irrelevant and, therefore, (1.7) relates y and 7 in 
a unique manner, where both » and 7 have definite physical meaning as to sign, 
direction of flow and sense of rotation. 
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1.2 Natural Coordinates 


In this system the coordinate lines in the meridional planes, 6 = const, are 
given by y(r, z) = const and the family of curves M(r, z) = const normal to y. 
The third direction 0 is the normal to @ = const; unit vector uy: (0, 1, 0). 

Taking for positive sense along the p-curves that of the flow, the correspond- 


ing unit vector is 
C 


m(ce, 0, ct (1.10) 


and for positive sense along the M-curves that of increasing wy, the unit vector 
is u, of (1.8) and the cyclic order Uy, Uy, U, gives the same orientation of the 
space as 7, 0, z in all cases. 
In this system 
Cage Cae) a (1.11) 


C 
W: (Cy, W,, 0) (1.12) 


and the components of 22: (E, n, f) in this system are easy to obtain from 
(1.2), (1.3), (1.4), (1.8) and (1.10) 


&§=2-uy = — grad, g = —Cy g,s (1.139 
n=, (1.14) 

5 il 
f= 2 use bo (1.15) 


As it is shown in Figure 1, the selection of sense in a natural system deter- 
mines the sense that must be given to the corresponding system of cylindrical 
coordinates. 


2. The Force Field 
2.1 General Properties 


If an axially symmetric flow is to provide a model of the flow through a 
turbo-machine, it is necessary to assume the existence of a field of forces F, 
acting upon the fluid and replacing a given blade row. The angular component 
J, of F must naturally be different from zero. 

From a physical point of view the problem is then to define adequately the 
force field F'so as to represent a blade row. 

Since the fluid is assumed inviscid, the force field must be non-dissipative, 


i, €., it must do no work in the absolute motion through a stator-field, or in the 
relative motion through a rotor-field. 
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Since 
Cx 2-Wx2=opradg, (2st) 


the vectorial equation of motion is, for both stators and rotors, 
grad H* =Wx Q2+4F. (2,2) 
However, for the motion through a stator field it is more convenient to write 
gradH = Cx QiF. (2.2a) 


By dot multiplication in (2.2) and (2.2a) by W and C respectively, there results 
W - grad H* = 0 


and 
C grad 7 = 0. 
Consequently 
H* = H*(y) (2.3) 
and 
H = Hy). (2.3a) 


The converse is also true, and it may be proven by retracing the above steps, so 
that the following statement holds: 

_In the motion through a rotor, H* = H*(w), and in the motion through a 
stator, H = H(w) represent the necessary and sufficient conditions for the force 
field to be non-dissipative. 


2.2 The Force Field and the Vorticity Field 


First the case of rotors will be considered and (2.2) will be written in scalar 
form. 


O=W, 6+ Fy, (2.4) 


ae ae gy aa (2.5) 
grad, H* = Cy 4 —W,é+F,. (2.6) 


Then, take the tangent plane to a streamsurface at a point O. Evidently, this 
plane is determined by the tangent to the meridional streamline through O 
and the whirl direction, as it is shown in Figure 3 below. 

The vorticity field and the force field may be decomposed into two parts, 
along u,, and on the tangent plane, 


Q=2,42,, F=F,+F,. 
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In natural coordinates, 


REE OG a Oia 


Fe a Fe 0 ae 
From (2.4) and (2.5), it follows immediately that ii, = 0, then € + 0 and 
consequently Fy + 0, so that both components of &, are different from zero 
if ¢ is different from zero. This relation will be termed ‘physical equivalence’ 
since the presence of one implies the other, and conversely. /¢ is clear that F 
and Q,, are physically equivalent. 2, will be termed the normal component of 
vorticity. The relation between the other two vectors, F,, and 9,;, is not so 


Figure 3 
5 


simple. In the tangent plane, 2, (Figure 3) may, in turn, be decomposed along 
and in the normal direction to the relative velocity, 


2,=2,+ 2, 
and since W x 922, = 0, equation (2.6) may be rewritten, 
grad, 7* == WOFer ie. (2.6’) 


where 2, has a sign, positive or negative, depending on the sense of the trihedron 
u,, Uy, U, with respect to Uy, Uy, U,. The component 2, of the vorticity along 
the relative velocity does not enter into the equations of motion and will be 
termed The Beltrami Component of 2. The component Q, will be called the 
binormal component of 2. These directions are not parallel to the normal and 
binormal of the streamlines, except in exceptional cases. 

In regard to 2, and F;, two cases are possible. 


(I) If H* = 0, Lorenz flows, then Q, and f;, are physically equivalent. 
(II) If Hf + 0, a first binormal component of Q, Q,, is defined such that, 


WQ, = H* 
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which is ‘free’ vorticity and requires the presence of no forces, and a second 
binormal component of 2, Q,., such that 


b» 


WO, = —F- with Q, + 2, = 
physically equivalent to F,. 


This analysis may now be summarized, and free and bound vorticity 
properly defined: 

If the flow is a Lorenz flow, H* = 0, the field of forces F is physically equi- 
valent to a system of bound vorticity along the normal and binormal directions. 
In this case, the free vorticity is the Beltrami component of Q. 

If the flow is not a Lorenz flow, H* + 0, the field of forces F is physically 
equivalent to a system of bound vorticity made up of the normal and the second 
binormal components of 2. The free vorticity is the sum of the Beltrami and the 
first binormal components of Q. 

In both cases the free vorticity represents the initial values for HELMHOLTZ’ 
equation r 
ag ee (QV)C 

ruling the variation with time of 2 outside the force field. 

In a stator field the same analysis applies if the Beltrami component is 
defined as the component of 2 along C, H is substituted for H*, and the so 
called binormal components are along the normal to C, instead of W, on the 
tangent plane. 


3. Dynamic Classification of Turbo-Flows 
3.1 Potential Force Fields 


As in classical mechanics, turbo-flow fields may be classified as potentzal 
and non-potential. Taking the curl in both sides of equation (2.2) or (2.2a), there 
results , 

curl (C x 2) + curlF = 0 (3.1) 


in both cases. If F is a potential field, 
curl(C x 2) =0. (3.2) 


That (3.2) is the characteristic vectorial equation of potential force flows 
may be easily proven. The direct theorem, ‘if a fluid mass moves steadily under 
the influence of a potential force field, then (3.2) holds’, has been proven above. 

The converse, ‘if (3.2) holds and 0C/0t = 0, the fluid motion defined by the 
velocity C takes place under the influence of a potential force field’, will be 
proven immediately. 
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Taking the curl on both sides of the identicy 


Borns 2 26 + grad (5 Cc) Cx 
it follows that 
curl (+ C) =0, 


, the curl of the acceleration is zero. Hence, the resultant of the forces 
sete on the fluid is a potential force field. Since the pressure forces per unit 
mass are given by — grad ib adp/o, the remaining forces, i. e., the force field F, 
must also be potential. 

Consequently, from the three'scalar equations corresponding to (3.1), which 
in cylindrical coordinates read 


1 Oe foe 
=>, (ape) =. ae 
d 2 1 
aa (2) SEE eS curl, F , (3.4) 
1 0 d 
ee (=; s) = CULE, (3.5) 


d 

Pir g 0! 5 (3.6) 
d (y 2 

ae (F) = 78 BB Ca 


where @ is a spatial constant. Thus, (3.6) and (3.7) are the characteristic equations 
of potential force flows in kinematic terms. 

For the axially symmetric case, from the condition curl F = 0, the charac- 
teristic conditions in scalar terms will now be deduced. 


From curl,F = 0 it follows that there must exist a function, f(r, 2), such 
that, 


F=f, and = fy 
From curl, F = 0 and curl, F = 0 it follows that 


0 0 
ae (rf) =0 and ae (7 F,).= 0 


so that 
rF,=k (k= const). 
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Since, in cylindrical coordinates, grady J*= Iy/r, there exists a family of 
surfaces of the form, 


=k [6 — flr, 2)] (3.8) 
such that 
Fis pradl'. (3.9) 


Inasmuch as 6 = + wt, for rotors J’ will be taken in the form 
r= k[p— fr, 2)] (3.10) 


and (3.9) will be reserved for stators. The converse is also true, if the force 
field is given by (3.9) with J” as in (3.10), or as in (3.8), the force field is potential 
and axially symmetric, (3.10) and (3.9), or (3.9) and (3.8), are the characteristic 
conditions, in dynamic terms, for the flow to be a potential force flow. The k in 
(3.10) ts the same constant as the a in (3.6). 


3.2 Properties of Potential Force Flows 


To investigate the different properties of potential force flows, the scalar 
equations corresponding to (2.2) will be written in cylindrical coordinates. 


W, ; 
H* = —* g,—C,—aj,, See 
i d 
O=-—.—- g+<, (3.12) 
H* = uo) g+7C,—af,. (3.13) 
(3.12) simply says that 
wee const . 


Since for axially symmetric flows the value of 2 g at any point P of the field 
is the circulation about the circle through P normal to and with center at the 
axis, the difference in the time rate of circulation about any two such circles ts zero. 
These circles are circuits of the second kind, i.e., not reducible to a point 
within the field. For a circuit of the first kind, reducible to a point within the 
field, KELvrin’s theorem holds in its standard form, i. e., the time rate of circu- 
_ tation about a circuit of the first kind ts zero. 

By use of these results, it is proven in standard text books that, af the force 
field is potential, the vortex lines move with the fluid. 

It will be shown that to have the a of equation (3.12) equals a constant is 


a characteristic property of axially symmetric potential force flows. Other 


- 


properties of this type of flow will now be investigated. 


_ ZAMP X/23 
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(A) The I-Family of Curves: The equation of motion (2.2), with F from 
(3.9) and (3.10) reads 


erad H* = Wx 2 + agrad [p — f(r, 2)] . (2.2)9 
W grad H* = H* = 0 (2.3) 


therefore, dot multiplication by W in (2.2’) yields 
a 
“dt. [on z)] =0, 


i. e., the surfaces y — f(r, z) = const contain the relative streamlines. The whirl 
components of the relative and the absolute velocity are, therefore, given by 
the following expressions, 


We (Ca de Calas (3.14) 
Gre Ctet, Glance (3.15) 


Since, in potential force flows, the vortex lines move with the fluid, the 
I-surfaces will be made up of vortex lines if they contain one. From (2.2) it is 
clear that the surfaces 


B= H* — [’= const (3.16) 


contain both the vortex and the streamlines. However, only when H¥ = 0 the 
surfaces J’ will also contain the vortex lines. Evidently, the condition 


2 grad J a0) 
implies, from equation (2.2), that 
2 - grad H* = 2 - a grad [pm — f(r, 2)] « 


Since a4 = grad y-@ and grad H* = H*-grad y, it is clear that Q-a grad I" 
can only be zero when H* = 0. 
The structure of a potential force flow field may be described as follows: 
(i). the flow 1s Lorenz, there exists a family of surfaces I’ made up of vortex 
and streamlines. Furthermore, these surfaces are the equipotential surfaces of the 
force field F. 
(1i) I f the flow is not Lorenz, the family of equipotential surfaces of the force 
field F 1s made up of streamlines. However, the family of surfaces made up of both 
the vortex and the streamlines is given by the family B of (3.16) 


’ 
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For a stator field the same conclusions hold if in the above statements the 
absolute velocity is written instead of the relative, and the I-family is given by 
(3.8) instead of (3.10). 

(B) Analytical Properties: If in (3.7) the differentiation as to ¢ is carried out, 
and use is made of (1.5) and (1.6), there results 


1 1 Z 
y2 (1), Ys "s Yr) eS Pz = yt § &:- (3.17) 
Changing independent variables from 1, z, to 7, y, 
yp=(7,2), rar (3.18) 


a transformation’) with discriminant 


and, therefore, valid for all non-cylindrical flows, C, = y,/r = 0. Use of (3.16) on 
any function, say 7, gives 


Ne = Yr 7 Ny Van la Ny Ve 


so that, in the new variables z and y, (3.17) becomes 


the general solution of which is given by 


‘oe 2s Foe Th (0) SR eS 
Sia ae Oe ae (3.19) 


n= 7 |Kly) 


If to (3.13) is applied the transformation (3.16) — i.e., using (1.5) and 
dividing by y, — there results 


H* = aa g By wet ee * a AS hy : (3.20) 


VY 


It is easy to show that (3.19) and (3.20) are identical if K(y) = Hy. To 
prove this, it is enough to show that 


Hi pd ee Ne Re at ie (3.21) 


v2 Or 


3) This is a variation of a transformation first used by von Mises in viscous flow problems. See 
Bemerkungen zur Hydrodynamik, Z. angew. Math. Mech. 7 ( 1927). This type of transformation will 
be used later on and plays an essential part in the reduction of the equation for the blade. 
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Taking derivatives with respect to 7 in both sides 
1 ie x —_ = — 
maT (2 Spr Bae pl Oey a lage (3.22) 


Writing g in terms of 3 


g=f,C,r7+or 


and taking derivatives with respect to p 
ac, 


ae (3.23) 


By = Try C, Ae Le Phe 


Be may be eliminated between (3.22) and (3.23). There results, taking into 
account that a=, C,., 


a & oC, pao em Se Go 
& C, lz, * le t ae = C, lg Sry 1 Be Be 7) [g = Le Cy al C, Sry ° 


This equation may be reduced to 


which is an identity. Therefore, in an axially symmetric flow field with H* — or 
H — a given function of w, the condition that the time rate of the angular momentum 
be everywhere constant is sufficient for the force field to be potential. 

This characteristic property holds without any change for potential flows 
through a stator. 


3.3 Non-Potential Force Flows 


It was seen that if the flow is a potential force flow, there exist two one- 
parameter families of surfaces dynamically significant. The [-surfaces, 


I= a[g — f(r, z)] = const 
made up of streamlines, and to which the force field is normal; and the B-sur- 


faces, 
B = (H* — I’) = const 


made up of vortex and streamlines. If HY is zero, i.e., if H* is everywhere 
constant, then both families of surfaces coincide. 

Consequently, in the general case, H » + 0, for a potential force flow, one is 
at a loss as to what family of surfaces should be taken as the set of ‘blades’ in 


am 
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the axially symmetric flow pattern. Since these two properties are not indivi- 
dually characteristic, potential force flows may be generalized in two directions. 
(A) von Mises Flows: In 1909 — see reference [2] — von Misks defined the 
‘blade’ in the axially symmetric flow pattern by the condition that it should 
be made up of vortex and streamlines. 
The conditions for a family of surfaces 


B= 9 — f(r, 2) = const 


to exist and verify von Mises’ definition of ‘blades’ may be obtained — see 
reference [3] — from equation (2.2) as follows: 
Dot multiplication by W equals zero gives 


Ne ong eel Oe oe a 


iv 
and dot multiplication by 2 equals zero gives 
part ue = 
2 grad B = ¥ 8 fy + = == £4 = 0, 
two equations which may be solved for f, and /,. The condition 


i) 0 
Or in Oz f,-=0, 


from which is derived the name of ‘integrability condition’, using (1.5) and 
(1.6), yields von MIsEs’ equation for rotors 


0 (g = o#*) oF ¥* 7 0 (= @ 4") g. 4-77 9 Wy 
; 3.24) 
0z | Y (2, Yz — & cs) or | ii (85 Yz — 82 Yr) | ( ; 


von MisEs proved that the converse is true, if the angular momentum verifies 
(3.24) there exists a family of surfaces, y — f(r, 2) = const, containing both the 
vortex and the streamlines. Since (3.24) is an equation for the angular momen- 
tum g in terms of wy and ¢ there results an automatic ordering of the boundary 
conditions. These boundary conditions are well ordered in the sense that the 
meridional flow pattern, defined by equation (1.7), can be determined first. 
The boundary conditions are: 

(i) Casing shape. 

(ii) Ring vorticity. 
(iii) Distribution of axial velocity at the entrance and the outlet to the rotor. 
(iv) The function H*(y) — or H(y) — for stators. 

(v) 


An adequate distribution of g and its derivatives over one stream surface 4). 


4) Inasmuch as equation (3.24) is non-linear this problem is rather delicate and will not be 
discussed here. See von Miszs’ original paper [2]. 
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von Mises’ equation for the flow through stators can be obtained in the same 
manner, multiplying (2.2a) by C and 2. The final equation is 


ete [ B br Yr a | bet 7 Us (3.24a) 
Oz | (8, Vz — 82 Yr) or y. (Bn Yz — &z Yr) ‘ 


which may be obtained directly from (3.24) by setting @ = 0. 
It is simple to verify that if 


1 
y- (2, Wz &z Yr) 


d — — 
dies a 


equals a constant, both (3.24) and (3.24a) become equation (3.7), i. e., potential 
force flows are a subclass of von Mises flows. 

Using exactly the same reasoning as for potential force flows, it will be 
shown that the force field — which in von Mises’ theory is determined as a 
consequence of the knowledge of the flow pattern — can only be normal to the 
surface containing the vortex and the streamlines if H* — or H,, — is zero. 

(B) Quasi- Potential Force Flows: These flows will be defined by the condition 
that F be of the following form, 


F = « grad [p — f(r, 2)] (3.25) 


and represent, therefore, a generalization of potential force flows, whenever « 
is a function of y and z with «, and «, not simultaneously equal zero. The 
surfaces m — f(r, z) = const will represent the ‘blades’ for quasi-potential 
force flows. 

Equation (2.2) with F as in (3.25) shows that 


W - grad [p — f(7, z)] = 0 


and, therefore, the surfaces of the family » — f(r, z) = const contain the stream- 
lines. It follows then that W, and Cy, are given by (3.14) and (3.15). 

These ‘blades’ will contain the vortex lines only if H , 80! 

Furthermore, since 


grad {a [p — f(r, z)]} + « grad [p — f(r, z)] (3.25a) 
it is no longer possible to assume that there exists a family of surfaces con- 


taining both the vortex and the streamlines. 


The study of these quasi-potential force flows, the main goal of this paper, 
will be made in the next part. 


It should be noted that if H* = 0, the condition that F be as in (3.25) is 
equivalent to VON MIsks’ equation. Evidently, for F, and F, to be of the form 


Ea fi teeta 


Pa 
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Ww 


it isnecessary — and sufficient — that 


#2 (4), 2 


If in (3.26) «, &, F are substituted for their expressions in EULER’s equations, 
(3.26) becomes von MIsEs’ equation. 

Consequently, and this has some historical importance, the solution of the 
problem of quasi-potential force flows implies the solution of a very important 
class of von Mises flows. 

It will be seen that the equations for the flow through stators are a special 
form — w= 0 — of those for rotors; the flow field has the same geometric 
structure as in the rotor field, depending now on the form of H,, instead of H*. 


4. Quasi-Potential Force Flows 


It has been shown that the determination of an axially symmetric turbo- 
flow depends entirely, for non-potential force flows, on the definition of ‘blade’. 
The adequacy of this model to represent actual turbo-flows is therefore depen- 
dent on that definition. 

The ‘blades’ of a quasi-potential force flow, for rotor fields, are the surfaces 
of the family 

gy — f(r, 2) = const (4.1) 


containing the streamlines and such that the force field F is given by 
F= alr, 2) grad le — f(r, 2)). (4.2) 


EULER’s equations are (3.11), (3.12), and (3.13) with « instead of a. 

The hydrodynamic problem may now be stated: 

‘What are the set of equations and the corresponding boundary conditions which 
give a unique solution satisfying Euler’s equations, and the equation of continuity, 


divC =0.’ (4.9) 


Furthermore, ‘it will be required that the meridional flow pattern can be determined 
first’. 

With these requirements, (1.7) is the first equation to be solved and the 
corresponding, complete boundary conditions are: 


(i), (ii) and (iii) 


~ of (3.3). The components of the velocity are given by (1.5) and (1.6), and the 


equation of continuity is automatically satisfied. 
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It will now be proven that ‘if H* = H*(w) and two Euler equations are 
satisfied, the third equation holds also. Furthermore, if im these two equations, say 
the second and third, F.= —af, and F, =alr, where both « and f are functions 
of r and z (f, and f,, and a, and «, not simultaneously zero), then there exists a 
family of curves such as (4.1) which contains the streamlines and F is given 
by (4.2).’ 

If H* = H*(y), then d/dt H* = 0,1.e., 
dCo 
at 


dC 
ae fo Mp 


(Gate ihe Cou 
ia eG F 


3 == Oe 
dt 2 at oleae 


OC, Caiowr 


Multiplying (3.12) and (3.13) by W, and c, respectively, adding and sub- 
tracting the sum from the above equation, there results 


(On ; 
Be: 5 a pS San ey te, (FW, — af. CG.) 


dt 


which is the first Euler equation with 


= Wo as ip C, 
Eo = =o C (4.4) 
Since H* = H*(w) implies that 
FC, FW, -peb, Cg 0 (4.5) 


and W, = r dq/dt, using (4.4), (4.5) becomes 


a (Sf — Gr) =0 (4,6, +4.6,- 42) =0 


Or 


20 (3? — 1) =0 


and, therefore, the family of curves 
@ — { = const 


contains the streamlines and F, is given by —af,, i.e., F is given by (4.2). 
The statement is proven. This proof is independent of which two equations are 
selected. | 

Consequently, once the solution of (1.7) is known, a set of two Euler equa- 
tions should be chosen and the boundary conditions for the set be defined so. 


that the solution is unique and compatible with the meridional flow pattern 
already determined. 


_ 


“— 


A"). 
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In what follows, the variables z and wp will be used, i. e., the transformation 
y=ypr,z), 2=2 (4.6) 


will be performed on the chosen equations. This transformation is similar to 
(3.18) and valid wherever y, + 0, i. e., for all flows which are not purely radial. 
In what follows, however, the bars will be omitted. Naturally, it must be under- 
stood that the same symbol, say C,, does not indicate the same function in + 
and z as in z and y; of course, both expressions for C, become identical by use 
of (4.6). 7 will indicate the function of z and y, obtained by solving for 7 in the 
first equation of (4.6). 

It will be shown in part 5 that, paradoxically as it may seem, z and y are the 
most adequate variables to express analytically the solution of (1.7) from 
numerical data. 

Equations (3.11) and (3.12), in the variables z and yp, read 


W, y 
Hate Mg) p testa f 5 (4.7) 
a=C,*.(Pf,C, + Po). (4.8) 


This is a system of two equations for the unknown functions « and /. Elimination 
of « between (4.7) and (4.8) gives the following equation for /: 


He 7 =C. (te op — fy Ge) Cet Pe | 
PCr Tye te led th Cale Cay ty OCs) | (4.9) 
1. C, alae sao ty(?”) 1 C 


The above equation for / is a non-linear partial differential equation of second 
order, the coefficients of which are all known, either from the original data, as 
H* and y, or from the solution of equation (1.7). The determination of «, to be 
made in equation (4.8), is done directly in terms of f and known quantities. 

In physical terms, this may be interpreted as follows: 

For a predetermined meridional flow pattern, the differential equation for the 
function f of (4.1), giving the blade shape, is obtained by elimination of the function x 
— load function — between the equations of motion along the w®) and 6-directions. 
The working load®), given by the function a, corresponding to that blade 1s then 
given directly in the equation of motion along the 0-direction. 


5) Equation (4.7) is the equation of motion along the w-direction, changed to the derivatives 
with respect to y-from those with respect to Sy — by use of (1.9). 

6) The term ‘working load’ applied to & is justified inasmuch as the change in angular momen- 
tum per unit mass is given in terms of « only, i. e., does not depend on F’,, and F,. See reference [4], 


Chapter I. 
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Consequently, the following statement is proven: In the axially symmetric 
model of turbo-flows, with quasi-potential force frelds, the meridional flow pattern 
is independent of the motion along the whirl direction. And, evidently, this 1s no 
longer true for potential force flows; for, the « im (4.8) will be a constant only 


‘accidentally’. Ste 
The manner in which (4.7) and (4.8) are obtained is simple: 


(4.7) is obtained from (4.11) by division by ¥y,. 
From (4.1) it follows that 


i d 
gah, 8=PhG+Po, Ge=aG, 


and substitution in (4.12) gives (4.8). 

The remaining question is, then, what are the boundary conditions giving 
a unique solution for (4.9). 

The differential equation of the characteristics of (4.9) is 


f, dz dp + f, dy? =0 (4.10) 
or 
dy (f, dz + f, dy) = 0 


and, therefore, the equation is of hyperbolic type and the characteristics are 
y=const, f=const. — (4.11) 


Consequently, the reduction of (4.9) to normal form is accomplished by — 
interchanging /, z, y, so that z be the dependent variable, i. e., (4.9) is, up to a 
factor, the differential equation in normal form defining the function 


z= 2(f,y). (4.12) 


Taking (4.12) as a transformation, the following formulae hold: 


ak 28 ht Bn Urges 
ers (ip) COnST) A Bois (f = const) . (4.13) 
From either one of (4.13) 
faa f —_ fzyplz 
254 = ts B fav te (4.14) 


Using (4.13) and (4.14), (4.9) becomes 
ty = (Lb — Ht)? C2 8 4 TC), & GAC rip 
fy Y y a / Zp (7 ey z,| Cc; KE SF 
(4.15) 
+ [(7*), + (7), Zp Cem eae 
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an equation which will be written symbolically, 


Shy aa D(z, Y, ey 24) ps (4:15") 


In this form — see reference [7] — a set of boundary conditions, complete from 
a mathematical standpoint and admissible from the point of view of appli- 
_ cations, is the following: 


(I) f=, (yp) at- z=0, 
(4.16) 
(II) f=%,(z) at p=0. ) 
As shown in Figure 4, a cartesian representation of z= 2(f, p), the first 


condition above is the equation for the leading edge (assuming the entrance to the 
rotor field at z = 0). The second condition gives the blade form at one streamline. 


y 

~ y rh, (y)= 0 

e 
leading edge \. 
as Second system of 
\ characteristics 
ce 
aN 
Lx 

fh ret 

ak _*d First system of 

Be eG characteristics 

er ai ok 
Ze ~~ a 
Y;3 ral <= 2 a a 4 
of ae Fh, (Z)=0 
Wy - Blade sage at y-0 ~~ 2% inv (f)] 
Y; s > SK 
N 
~ 
N 
6-f i 


Figure 4 
Cartesian representation of the characteristics and the boundary conditions. 


_ In mathematical terms, these conditions give the value of z over one characteristic, 
Bs — hin’(f) at y=0, and the value of z, z=0, over a non-characteristic line, 

f = hy(p). Picarp proved that for these conditions the solution is unique and 
_ may be obtained by repeated iterations, if D(z, yp, z;, 2,) is adequately bounded. 
_ See again reference [7] for a precise definition of these terms. 
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From (4.15), taking into account (4.16), condition II, integration with 
respect to y gives’) 


Z= hie’ (f) +[@ [2(f, Po), 2r(f, Po), Zy(f, Wo)» Yo! Wo - 


Integrating ie with respect to /, and taking into account condition I, 


” il hae: in + fort [2(f, Wo); 2,(f, Yo)» z,(f, Yo)> Pol apo af, - (4.17) 
hy(v) 
Since 
f : 5 
Hit” (fg) dy — H(A) — Ha yy) 
Iy(y) 


Calling hi"(f) = A(f) and Ay” [h,(w)] = Bly), (4.17) may be rewritten 


f y 
z= Alf) — Bly) 4 /| [© le 20% Vol M0 | Aho (4.18) 
0 


hy () 


where the first approximation, 2,, 1s given by 


z= A(f) — Bly), 


i. e., by the boundary conditions. Because this process is convergent, the non- 
linearity of the equation is no longer relevant. (If the iteration is carried out 
locally, the non-linearity of the equation restricts the validity of the process. 
This will be seen in part 5.) 

For stators, the stream function is the same as for rotors. Equation (4.1) 
reads 


6 — f(r, 2) = const . 


(4.2) is the same — with « instead of a, of course. (3.11) reads 


Co 
= spores fied re 1 (3.11a) 


(3.12) is also the same. Applying the transformation (4.6) to (3.11a) and (3.12) 
there results the system 


H, = mt n 0p, 
ner mis rams Dia R= Ort eG | 


7) hinv j is the inverse function of fy, that giving z as a function of if him’ is simply the deri- 
vative oh ie with respect to f. 


The subscript 6 indicating the variable of integration. 
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and the result of eliminating « is 


Wore 0 ae) ' 
Bote elie mip) orca, 


i. e., the equation for the stator blade is obtained from (4.9) by setting w = 0, 
and changing H* to H. 

In the form of an answer to the formulation of the problem given at the 
beginning of this part, a summary of these results is given below. 

“The hydrodynamic problem for quasi-potential force flows is reduced to the 
solution of STOKES’ equation for the stream function y, and the hyperbolic equation 
for the function f giving the blade shape that results from elimination of the load- 
function « between the equations of motion along the wy and 6-directions.’ The 
boundary conditions defining uniquely the flow are: 

(I) The distribution of energy, H*(y) — or H(y). 
(II) », the vorticity along the tangential direction. 
(III) The casing shape. 
_ (IV) The distribution of axial velocity at the entrance and outlet of the blade row. 
(V) The leading edge of the blade. 
(VI) The intersection of the blade and one stream surface. 


5. An Example 
5.1 The Meridional Flow Pattern 


Boundary conditions II, III, and IV determine uniquely the meridional 
flow [9]. 
In regard to III, the casing shape was selected with Ry = const and the inner 
radius from the general expression 
Y x—1 


= Gas eRe Bee ean 


cn 


36 —2) 22+ (412.6%) 8+ (65-3) 24] (5.1) 


giving a family of curves with zero tangent at R;, and R;,; z varying from zero 
to one; and such that z=0, r= R,, and for z=1, r= RR; *#= Raliee: 
_ éis the z-coordinate of the inflection point. 
Figure 6 is a casing for the following specifications: 


Pei Sa 05h, 22) Rees 

2 The geometric configuration is determined by the ratio of any one radius to 
the axial length. This ratio was chosen, 

3 gs 3 Bape 18 

4 ae, Oe ete ee 


» a beef aed bo She 
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That the concrete specification of this ratio is necessary is seen by inspection 
of equation (5.5). 

n was selected from experimental work in the following manner. From 
HowELL’s measurements [8] of the axial velocity distribution through an axial 
compressor, taking differentials in 


(ca) (cs) 
pam is transformed into ee 
(7) 


By graphical differentiations, and assuming that HowELt’s five curves of 
Figure 5 correspond to z = 0, 0-25, 0-5, 0-75, 1, there results a distribution of 


1-0 


—_— 


8 
6 


] 

0 

0-4 [ 

0-2 A 

005 ZZ 


00 02 04 + 08 88. GIO. 42 
Axial velocity/mean axial velocity ——— 


Blade height —— 


Figure 5 


Distribution of axial velocity through a four-stage compressor (from Howe tt). 
Curve O after entry guide; J after first stage; 2 after second stage; 3 after third stage; 4 after fourth 
stage. 


Figure 6 


Case shape and net code arrangement. 


a. i ee ee 
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meridional vorticity, 7, throughout the casing as shown in Figure 9. This vor- 
ticity distribution was obtained under the assumption 


Physically, 7/r has more meaning than , for it represents the strength of 
the vortex ring [11]. Consequently, three cases are possible: 
Case I: H=0. 
Case IT: y/r = function of 7 and z. 
Case III: »/r = function of y. 


To obtain 7/7 as a function of y, curve 2 of Figure 5 was used under the 
assumption (5.2). This function is shown in Figure 7. 


0-2 


Figure 7 
Vortex ring intensity corresponding to curve 2 of Figure 5. 


Figure 8 
Entrance and exit-velocity distribution corresponding to the vorticity above. 
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Although (5.1) is given in dimensionless z, it is necessary to determine the 
ratio of one radius to the axial length, L. This chosen ratio, R;,/L, was set equal 
to 4/3, giving L = 16/3. 

For 7 = 0, the axial velocity was assumed constant with 7, at the entrance 


and the outlet. 
For cases II and III, the axial velocity distributions at the entrance and 


the outlet are given in Figures 10 and 8, respectively. 


Figure 9 
Vorticity distribution as derived from Figure 5. 


0 7 2 3 


os aa 
ifs 


Figure 10 
Entrance and exit-velocity distribution corresponding to the vorticity above. 


For the net code arrangement shown in Figure 6, the numerical values of p 


are given in Table 1. These values are used : : 
. as a means to obta 
expressions for the solutions, in analytical 
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Table 1 
Case I Case II Case III 
Pt. No. p-values Pt. No. p-values Pt. No. p-values 
1 0-0356 1 0-0318 1 0-0270 
2 0-0381 2 0-0326 2 0-0292 
3 0-0418 3 0-0346 3 0-0329 
4 0-1023 4 0-0965 + 0-0892 
ia) 0-1094 5 0-1017 5) 0-0966 
6 0-1199 6 0-1116 6 0-1081 
vi 0-1629 ‘A 0-1576 7 0-1510 
8 0-1736 8 0-1686 8 0-1624 
9 0-1895 9 0-1856 9 0-1801 
10 0-2173 10 0-2137 10 0-2098 
11 0-2303 11 0-2282 11 0:2240 
12 0-2499 12 0-2504 iy 0-2455 
13 0-2655 13 0:2637 JT) 0-2625 
14 0-2793 14 0-:2796 14 0-2780 
ED 0-2999 15 0-3037 15 0-3004 
16 0-3073 16 0-3073 16 0-3087 
1y; 0-3209 a 0:3226 Li 0-3236 
18 0-3477 18 0:3512 18 0-3505 
19 0-3427 19 0:3442 19 0-3460 
20 0-3570 20 0-3589 20 0-3595 
21 0-0245 21 0-0190 21 0-0199 
22 0-0284 faye 0-0215 22 0-0237 
23 0-0337 23 0-0248 23 0-0289 
24 0-0388 24 0-0270 24 0-0336 
79) 0-0710 25 0-0600 25 0-0617 
26 0-0833 26 0-0700 26 0:0747 
27 0-0990 Bip 0-0817 27 0-0912 
28 0-1142 28 0-0879 28 0-1067 
29 0-1154 29 0-1048 29 0-1051 
30 0-1358 30 0-1235 30 0-1266 
Sil 0-1617 31 0°1473 31 0-1541 
32 0-1868 32 0-1663 32 0-1808 
33 0-1575 33 0-1498 33 0-1478 
34 0-1856 34 0-1786 34 0-1780 
35 O-2217 35 0-2176 35 0-2168 
36 0-2572 36 0-2543 36 0:2546 
29) 0-1972 37 0-1933 37 0-1895 
38 0-2326 38 0-2325 38 0-2278 
39 0-2790 39 0-2858 39 0-2773 
40 0-3257 40 0-3308 40 0-3268 
41 0-2343 41 0°2342 “isl 0:2293 
42 0.2766 42 0:2832 42 0:2748 
43 0-3336 43 0-3409 43 0:3348 
44 0-2687 44 0-2710 44 0-2665 
45 0-3174 45 0-3252 45 0-3186 
46 0-3003 46 030177, 46 0-3009 
47 0-3548 47 0-3585 47 0-3565 
48 0.3290 48 0-3316 48 0-3315 
49 03549) 49 0:3571 49 0-3570 
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By cross-plotting the values of p(7, 2), a set of streamlines — or stream- 
surfaces — are obtained for all three cases, Figures 11, 12, and 14. 


Figure 11 Figure 12 
Irrotational meridional flow pattern (case I). Rotational meridional flow pattern (case IT). 


Case I! Case II 
y=0 y=0- 

0.0960 
0:1956 
0:2988 
0:3691 

a 

Figure 13 Figure 14 
Comparison of meridional flow patterns Comparison of meridional flow patterns 
(cases I and II). -—~—case I; — case II. (cases II and III). — case Il; ——-— case III. 


For the first two cases, along these streamlines 7 as a function of 2 is 
expressed in Fourier series. Then, the coefficients are functions of wy, and are 
expressed as polynomials. These analytical expressions are given below. : 


Case I 
r = (5-002 — 4-203 p — 0-0125 p® — 8-857 pi) : 
0-004 — 1-879 p + 2-121 wy — 11-571 side! | 
tak 1-879 p + 2-121 py? — 11-571 w’) cos > — (5.3) 


Se 


+ (— 0-136 p + 0-494 p? — 0-553 wp) cos, 
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Case II 
ry = (5-0 — 5-148 p + 4-011 py? — 12:84 p) 


+ (— 1:365 y — 0-5067 p® — 8-484 3) cos = 


B w . 2 f =a a S08 
+ (0-1740 p — 3-414 wy? + 7-959 p>) sin * 5.4) 
+ (0-3709 p — 3-996 p? + 8-096 wp) pos 


+ (— 0:1363 p + 3-170 y — 7-613 #) ne 
pes L=5 + unit 
2 = 4 units , 


ry, z real coordinates, C, average axial velocity at the entrance. 
In reality (5.3) and (5.4) are the solutions of the equation into which 


1 
Pre — Y Pr + Yaz = 7H (5.5) 


transform when z and w are taken as independent variables, i. e., (5.3) and (5.4) 
are the solutions of 


2 hee 2 Ay 2 73 == 
ak a a te ee a OP Pe ee ht hae A/a (5.6) 


for the specific boundary conditions. 
The advantages of such conversion of the numerical data to analytical form 
_ are manifold: 
It avoids graphical differentiations which, tf of second order, can be so inac- 
curate as to be almost meaningless. 
It provides the solution in the needed form for the blade equation. 
Finally, (5.6) may be interpreted — after division by r?r}—as an error 
_ function. This error function has dimensions of cycles per second per unit of 
average axial entrance velocity. It should be noted that, having only the 
numerical values, it would be impossible to estimate the error. A map of this 
error function is shown in Figure 19. 
In Figure 13, the flow patterns for cases I and II, as obtained from (5.3) and 
_ (5.4), are compared. All the velocity distributions, Figures 15, 16, 17 and.18; 
_ are derived from (5.3) and (5.4). 
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Figure 17 Figure 18 
Axial velocity corresponding to flow pattern Radial velocity corresponding to flow 
of Figure 12. pattern of Figure 12. 
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Figure 19 
Map of error function with the same dimensions as IC ay L?. 
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5.2 The Blade Shape 


To illustrate the method and to simplify computations, it is only necessary 


to consider a stationary ‘blade’. The boundary conditions V and VI of part 4 
are specified as follows: 


Over the hub, i. e., in the (f, z)-plane of Figure 4, the distribution of turning 


is given by 


»" 


Xx.AT 


» scl tne sade Tht LoS the yes tS 


f = 0-1713 2? — 0-009532 28. (5.7) 


The leading edge is taken along the y-axis at the entrance, which for case I 
coincides with the radial direction. 

Table 2 gives the values of z over the (f, y)-plane as obtained directly, by 
iterative use of equation (4.18). This map shows why the blade, as given in the 
final drawing, Figure 23, does not go over the whole axial length of the hub. 


Table 2 
Values of z on the (wp, f)-Plane as Computed Directly 
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From the numerical values obtained, the curves of { versus z for different 
y-values are given in Figure 20. To illustrate the results, a ‘blade’, assembled 
from these data, is shown in Figure 23. 
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Figure 20 Le 
Angle variation through the blade. Figure 22 


Angular momentum distribution. 


Once the ‘blade’ is known, the distribution of Cy and of the angular momen- 
tum, C,7, can be calculated. This distribution is given in Figures 21 and 22, 
having been calculated as explained in part 4. 


The flow pattern then is completely determined. 
The pressure distribution depends, naturally, on the physical conditions at 
the entrance as described by H(y), which in the example under consideration 


is assumed constant, i. e., H,,(p) = 0. The actual values are to be obtained by 
substitution of C,, Cg, C, in 


1 
2 Cee rig const . 
@ 
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In Figure 23 the ‘blade’ is shown as assembled from all the computed data. 
The medium heavy lines are the lines of constant z and constant wy making up 
the blade. The heaviest line is the vortex line corresponding to 


drawn as if the computational errors, evident in Figure 20, were absent. The 
vorticity along this line is made up of free and bound vorticity. If the ‘blade’ 


Figure 23 


Trimetric projection of blade. 


were to be cut along this line, the ‘free vorticity’ shed would be the Beltrami 


component, i. e., the projection of 2 on the streamlines. 


Finally, these computational results show the inadequacy of choosing a 
radial entrance edge for the distribution (5.7). That is, the interrelationship 


_of the different boundary conditions defining the flow through a turbo-machine 


s 
y 
mn 
ee 
J 
4 
a 
~ 
2 


4 
~ 


is made evident. 
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Zusammenfassung 


Ausgehend von der Tatsache, dass eine Schaufelreihe 4quivalent ist einem 
Kraft- und damit auch einem Wirbelfeld, wird der Begriff der Potential-Kraft- 
Stré6mung erweitert. Diese schliesst den «von Mises-Fall» ein, wenn das Euler- 
Integral im ganzen Feld konstant ist. 

Die Struktur des Wirbelfeldes, sowohl des freien als auch des gebundenen, 
wird analysiert. Das Problem der Turbinenstrémung reduziert sich auf die Lé- 
sung von zwei partiellen Differentialgleichungen, fiir die ein vollstandiger Satz 
von Randbedingungen gegeben ist; die Rechnung wird auf bekannte numerische 
Methoden zuriickgefiihrt. 

Wahrend die iiblichen Rechnungsmethoden nur Eintritts- und Austritts- 
diagramme liefern, erhalten wir hier die vollstandige Geschwindigkeitsverteilung 
im ganzen Kraftfeld. Die Schaufelform wird angegeben, und es ist moglich, das 
Zusammenwirken verschiedener Entwurfsparameter zu untersuchen. 


(Received: September 20, 1958.) 
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Der Hagelversuchskanal 


Von Roranp List, Weissfluhjoch-Davos!) 


1. Einleitung 


Das Problem jeder Hageltheorie besteht darin, samtliche bekannten Phano- 
mene der Hagelbildung in ein mehr oder weniger geschlossenes System zusam- 
menzufassen und zu erklaren. Als Musterbeispiel sei hiefiir die noch heute 
giltige Theorie erwahnt, die TRABERT 1899 aufstellte [1]2). Die darin gemachten 
prinzipiellen Angaben geniigten den Bediirfnissen bald nicht mehr, wollte man 
doch verschiedene Einzelprobleme mathematisch darstellen. Mit der Entwick- 
lung der Aerodynamik und der Thermodynamik wurden daher grissere Studien 
uber Warmeiibergange und Strémungsverhiltnisse [2,3] angepackt. FosTER 
und Bates [4], Mason [5] und Luptam [6] — um nur einige zu nennen — 
ubernahmen die zum Teil weitgehend idealisierten Grundlagen und bauten die 
Theorien weiter aus. — Im Schosse der Eidgenéssischen Kommission zum Stu- 
dium der Hagelbildung und der Hagelabwehr reifte indessen der Entschluss, 
das physikalische Grundproblem unabhangig von den bestehenden, mehr theo- 
retischen Arbeiten von der experimentellen Seite her anzupacken. Man hoffte, 
durch geeignete Messungen bessern Einblick in die Vorgange des Wachstums 
der Eisbildungskerne zu Hagelkérmern zu erhalten. Als Hauptzweck ist dabei 


die Abklarung der Bedingungen in einer Gewitterwolke zu nennen, die zu 


schwerem Hagel fiithren. Diese Problemstellung ist, abgesehen vom reinen For- 
schungsinteresse, speziell fiir die Hagelabwehr von Bedeutung, soll doch damit 
abgeklart werden, wo am ehesten Ansatzpunkte fiir eine Hagelbekampfung zu 
finden sind. 

Als Untersuchungsmethode fallt die direkte aviatische Messung an Ort und 
Stelle der Hagelbildung infolge der ausserordentlichen Gefahrlichkeit solcher 
Fliige sowie der ungeniigend entwickelten Messmethodik dahin. Es soll daher 
versucht werden, auf Grund von Untersuchungen an natiirlichen Hagelkérnern 
Riickschliisse auf die urspriingliche Gewitterwolke zu erhalten. Studien an 
atmospharischen Eispartikeln [7-9] zeigen die ausserordentliche Vielfalt der 
Hagelkornstrukturen. Diese stellen prinzipiell das Produkt ihrer Entstehungs- 
geschichte dar. Eine Interpretation derselben ist jedoch héchstens in qualita- 
tiver Richtung méglich; sichere und genaue Angaben sind nur iiber verglei- 


1) Forschungsstelle der Eidgendssischen Kommission zum Studium der Hagelbildung und der 
Hagelabwehr am Eidgendéssischen Institut fiir Schnee- und Lawinenforschung. 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 413. 
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chende Experimente zu erwarten. Zur Erklarung des Aufbaues natiirlicher 
Hagelkérner sollen daher bei messbaren, den atmosphiarischen Verhaltnissen 
analogen Bedingungen kiinstliche Eiskérner erzeugt werden. Die Apparatur, die 
diese Nachbildung gestattet, ist der Hagelversuchskanal, ein im Messteil verti- 
kaler, kiinstlich klimatisierbarer Windkanal mit geschlossenem Kreislauf. 


2. Bedingungen zur Schaffung einer «atmospharischen» Umégebung fur 
ein Test-Eiskorn 


2.1 Zur Berechtigung der Betrachtung eines einzelnen Evskornes 


Die Méglichkeiten zur Nachbildung der Atmosphare im Laboratorium sind 
in erster Linie raummiassig beschrankt; sobald aber nur ein Hagelkorn als fiir’ 
sich isoliert betrachtet werden kann, eréffnen sich Wege zur Verwirklichung 
der Nachbildung im Experiment. Diese Annahme scheint gentigend gerecht- 
fertigt zu sein, wenn man die nach schweren Hagelfallen meist geringe Bedek- 
kung des Bodens mit Eisschlossen in Betracht zieht. In die gleiche Richtung 
weisen auch die vielfach zu beobachtenden, erstaunlichen Symmetrien der 
Hagelk6rner [9]. 

2.2 Wahl der Experimentparameter 


Das Wachstum eines Hagelkornes wird bestimmt durch seine Umgebung. 
Diese lasst sich charakterisieren durch Lufttemperatur, Luftfeuchte, Luftdruck, 
den freien Wassergehalt in Form von Tropfen oder Eisteilchen, deren Gréssen- 
verteilung und Temperatur sowie das Eisbildungs- und Kondensations-Kern- 
spektrum. Unter Vernachlassigung allfalliger Beschleunigungen kann weiter 
eine Relativ-Geschwindigkeit zwischen Korn und Umgebung von der Grésse 
der freien Fallgeschwindigkeit der Eispartikel angenommen werden. Die Um- 
gebungsverhiltnisse kénnen sich dadurch zeitlich rasch andern; mit Gleichge- 
wichten von Temperatur und Wasserdampfdruck zwischen Korn und Umge- 
bung ist nur in wenigen Fallen zu rechnen. 

Von verschiedenen Autoren wurde vermutet, dass auch elektrische Effekte 
das Wachstum von Hagelkérnern wesentlich beeinflussen. Da aber die Hagel- 
bildung nicht unbedingt an Gewitter mit Blitz und Donner gebunden ist (siehe 
Hagelschlag vom 29. Juli 1958 in Luzern mit Hagelk6rnern von tiber 10 cm 


Durchmesser), soll diese Méglichkeit als sekundar aufgefasst und vorlaufig nicht. 
als Parameter vorgesehen werden. 


2.3 Dimensionierung der Experimentparameter 


Die Dimensionierung der Experimentparameter soll so geschehen, dass die 
Versuche im Laboratorium auch bei noch extremeren Bedingungen durchge- 
fiihrt werden kénnen, als sie in der Natur zu erwarten sind, Da die Wolken- 
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gipfelhdhe selten itber 10000 m iiber Meer liegt, ergeben sich folgende Para- 
meterbereiche : 

Der Druck variiere zwischen 760-und 200 mm Hg. 

Die Temperatur liege zwischen + 25° und — 40°C, 

Der maximale freie Wassergehalt betrage 10 g/m, 

Die Tropfengrisse bewege sich von 1 p bis 200 pu. 

Unter der Beschrankung, im Experiment an einem zur Beobachtung giin- 
stigen Ort Eisteilchen von maximal 2,5 cm Durchmesser zu untersuchen, 
diirfte eine Luftgeschwindigkeit von maximal 25 m/s angemessen sein. 

Ein Ausschweben der Testobjekte wire erwiinscht, wird aber nicht zur Be- 
dingung gemacht, zeigten doch die Versuche am Aerodynamischen Institut der 
ETH die ganze Schwere dieses Unterfangens [10]. BLANCHARD, der schon frither 
Schwebeversuche mit Wassertropfen im offenen Windkanal machte [11], ge- 

- langte infolge der Selbststabilisierung durch Anpassung der Tropfenform und 
allfallige innere Ausgleichsstr6mungen des Wassers leichter zum Erfolg, als es 
fiir feste K6rper zu erwarten ist. 

Auf die Festlegung der restlichen Gréssen sei zum vornherein verzichtet, da 
nur das Experiment die im Laboratorium mit verniinftigen Mitteln erreich- 
baren Eisbildungskern-Charakteristiken und Unterkiithlungs- sowie Ubersiatti- 
gungsgrade aufzeigen kann. Gewisse Vorversuche zeigten aber bereits, dass 
Wassertropfen in Form von Nebel vorlaufig bis —17°C unterkiihlt werden 
konnten [12], und dass eine Variation der Gefrierkernkonzentrationen innerhalb 
weitester Grenzen méglich ist [13]. Unter diesen Voraussetzungen ist zu hoffen, 
dass auch hier atmospharische Bedingungen angendhert werden kénnen. Die 
Gewissheit hierzu wird jedoch erst durch einen Vergleich zwischen kiinstlich 
erzeugten und natiirlichen Eispartikeln erbracht werden kénnen. 


2.4 Der Verzicht auf die Druckvariation im Experiment 


An und fiir sich ware es vorteilhafter, die Druckvariation im Experiment zu 
verwirklichen. Infolge des damit verbundenen, unverhdaltnismassig grossen 
finanziellen Aufwandes wird aber darauf verzichtet und versucht, die Verhalt- 
nisse fiir alle méglichen Drucke mit Hilfe von Ahnlichkeitsgesetzen auf Grund 
des Laborluftdruckes (entsprechend 2665 m tiber Meer) zu berechnen. 

Wie aus der Strémungslehre bekannt ist, besteht die Méglichkeit, aus Mes- 
sungen an einem Versuchsobjekt auf die Verhaltnisse an einem andern geo- 
metrisch ahnlichen Korper zu schliessen — die Bedingung ist nur, dass die Rey- 
nolds-Zahl Re in beiden Fallen dieselbe ist. 

Solche Ahnlichkeitsgesetze bestehen ebenfalls fiir Warme- und Stoffiiber- 
gange — Vorginge, die das Wachstum von Hagelkornern wesentlich mitbeein- 

- flussen — und werden durch die Nusselt-Zahl Nu, beziehungsweise durch die 
Analog-Nusselt-Zahl Nw’ charakterisiert. 
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Die Druckabhangigkeit dieser Kennzahlen kann nach den zusammenfassen- 
den Arbeiten von Ranz und MARSHALL [14, 15] fir glatte Kugeln wie folgt 
implizit dargestellt werden: 


Nu = 2,0 + 0,60 Pril8 Ret?, (1) 
Nu’ = 2,0 + 0,60 Scll8 Ret?. (2) 


Dabei sind: 


Nu =ad/d, Nusselt-Zahl; 
Pr =»/a, Prandtl-Zahl; 
Re =v d/y, Reynolds-Zahl; 


Nw’ = pL ee , Nusselt-Zahl zweiter Art oder Analog-Nusseltzahl; oder 
CL 


Nu’ = Balk; 

Se =v/k, Schmidt-Zahl; 

«a Warmeiibergangszahl [kcal/m? h °]; 

d  Léange [m]; 

2  Warmeleitzahl von Luft [kcal/m h °]; 

y  kinematische Zahigkeit [m?/h]; 

a  Temperaturleitzahl [m2/h]; 

v __ Luftgeschwindigkeit [m/h]; 

k,, Masseniibergangszahl [kg/m? h Ap M]; 

Ap Partialdruckdifferenz des Wasserdampfes langs dem Transferweg [mm Hg}; 

M_ Molekulargewicht des titbergehenden Stoffes; 

M mittleres Molekulargewicht des Gasgemisches im Transferweg; 

p  Mittelwert aus der Differenz Totaldruck minus Partialdruck des diffundie- 
renden Dampfes im Transferweg [mm Hg]; 

k  Diffusionszahl von Wasserdampf in Luft [m?2/h]; 

0, Dichte des Luft-Gas-Gemisches im Transferweg [kg/m*]; 

B  Stoffiibergangszahl [m/h oder kg/m? h AC]; 

AC Konzentrationsanderung [kg/m3}. 


Wenn wir aerodynamische Versuche machen mit einem Testobjekt, das 
einen Durchmesser dy aufweise, so kénnen wir Schliisse ziehen auf ein geo- 
metrisch ahnliches bei anderem Druck. Soll die Reynolds-Zahl gleich bleiben, 
so wird dadurch der Durchmesser d, des letzteren bestimmt. Da fiir die ge- 
wiinschte physikalische Ahnlichkeit der Warme- und Stoffiibergange die Kenn- 
zahlen Nw und Nu’ konstant bleiben miissen, hat dies zur Folge, dass im einen 
Fall die Prandtl-Zahl Pr, im andern die Schmidt-Zahl Sc durch eine Druckvaria- 
tion gemdss Abschnitt 2.3 nicht beeinflusst werden diirfen. — Die Berechnung 
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dieses Einflusses auf die Kennzahlen Py und_Se ergibt Schwankungen, die 
kleiner als 1°, sind. Das heisst, dass unsere Forderung praktisch erfiillt ist — die 
Berechtigung fiir Extrapolationen von Experimenten aus dem Labor in die 
nattirliche Atmosphire ist damit fiir glatte Kugeln gegeben. Es sei hier noch 
speziell vermerkt, dass es sich hierbei nur um Luft-Luft-Transformationen 
handelt! Aus Form-Betrachtungen (siehe zum Beispiel [9, 16]) wissen wir aber, 
dass Hagelkérner héchst selten durch glatte Kugeln zu charakterisieren sind; 
wir miissen also auf die allgemeinen Gleichungen fiir beliebige Kérper zuriick- 
greifen. Diese lauten wie folgt: 


Nu =}(Pr, Re) , (3) 
Nu’ = f'(Sc, Re) . (4) 


Oben erwahnte Uberlegungen sind damit ebenso giiltig fiir den allgemeinen 
Fall; die Nusselt-Zahl Nw und die Analog-Nusselt-Zahl Nu’ bleiben gleich bei 
gleichem Re und Druckvariationen in dem von uns geforderten Bereich. 

Die physikalische Ahnlichkeit gilt damit auch fiir den allgemeinen Korper. 

Neben den Problemen der Str6mung, des Warme- und des Stoffaustausches, 
die durch die Kennzahlen Re, Nu und Nw’ charakterisiert werden, spielt auch 
das Einfangen freien Wassers durch das Hagel- oder Testkorn eine Rolle bei 
der Betrachtung des Wachstums einer Eispartikel. Wie weit diese Einfang- 
wahrscheinlichkeit druckabhangig ist, kann nicht entschieden werden, es be- 
steht jedoch die Hoffnung, dass sie bei gleicher Reynolds-Zahl ebenfalls gleich 
bleibt. 


3. Standort des Hagelversuchskanals 


Im Programm der Eidgenéssischen Kommission zum Studium der Hagel- 
bildung und der Hagelabwehr war von allem Anfang an die Forschungsstelle 
am Eidgendéssischen Institut fiir Schnee- und Lawinenforschung zur Bearbei- 
tung der Probleme iiber die kiinstliche Erzeugung hagelbildender Verhaltnisse 


- vorgesehen. Eine langere Tradition dieses Institutes in den Untersuchungen an 


festen Niederschlagspartikeln sowie gréssere praktische Erfahrungen tber 


_ atmosphirische Vereisungen in der Natur und im Windkanal waren die Griinde 
hierzu. 


Durch eine wesentliche Vergrésserung des Institutes auf Weissfluhjoch 


~ wurde der Laboratoriumsraum zur Aufnahme der geplanten Anlage geschaffen ; 
_ Kalterdiume zur Eisuntersuchung standen bereits mehrere zur Verfiigung. 


- 
‘a 


Die Wahl des Standortes ist auch giinstig, da zufolge der Hohe dieser For- 


_ schungsstation (2665 m tiber Meer) mindestens zehnmal pro Jahr Graupelfalle 
zu beobachten sind. Die Untersuchung der Reif- und Frostgraupeln — Ur- 


~ sprungspartikeln der meisten Hagelkérner — brachte im Laufe der vergangenen 


_ Jahre eine wesentliche Klarung dieser speziellen Entwicklungsstufen des hau- 


a 
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figsten Hagelkornwachstums. Die Ansetzung des Versuchsprogrammes mit dem 
Hagelkanal wurde damit wesentlich erleichtert. 


4. Aufbau und Wirkungsweise des Hagelversuchskanals 


4.1 Ubersichtsschema 


Die Hauptcharakteristik des Hagelversuchskanals ist der Darstellung in 
Figur 1 zu entnehmen: Ein Ventilator 7 treibt die Luft durch den Temperie- 
rungs- und Filterteil 2 in ein Winkelstiick, dem die 2m lange Befeuchtungs- 
und Ausgleichsstrecke 3 folgt. Dieser angebaut ist die eigentliche MeBstelle 4, 
in der die einzelnen Experimente stattfinden. Der Luftstrom gelangt hernach 
durch einen Diffusor in die eigentliche Riickfiihrung 5 und zuriick zum Venti- 


J 
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Figur 1 
Der Hagelversuchskanal (1:90). 7 Ventilator: 2 T, 
und Ausgleichsstrecke; 4 MeBstelle; 5 Riickfiihru 
keitsabscheider; 9 Ammoniakkondensor; 10 Ko 


emperierungs- und Filterteil; 3 Befeuchtungs- 
13 Kompensationsschreiber; 14 Te . 


ng; 6 Steuerpult ; 7 Kaltekompressor; 8 Fliissig- 
mpressorkiithler; 71 Olabscheider; 12 Motorventill 
mperaturregelanlage; 15 Rostbéden. | 
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lator, womit der Kreislauf geschlossen ist. Die gewiinschten Feuchtebedingun- 
gen werden durch Einspritzen von Dampf oder Wasser im Anfang der Aus- 
gleichsstrecke 3 erzeugt; am gleichen Ort geschieht auch die Zugabe der Eis- 
bildungskerne, die die Vereisungs-Charakteristik der an die Versuchsobjekte 
angestromten Feuchte bestimmen. Die Steuerung der Anlage erfolgt voneinem 
Pult 6 aus, wo auch samtliche wichtigen Gréssen angezeigt und registriert 
werden. 7 stellt den Kaltekompressor dar, 8 bis 72 dessen Zubehérteile. 


Figur 2 


Der Temperierungs- und Filterteil einschliesslich Befeuchtungsstrecke (1:50). 1 Geblase; 

2 Ammoniakverdampfer; 3 Heizelement; 4 elektrostatisches Luftfilter; 5 Eisbildungskerneinftih- 

rung; 6 Wassereinftihrung; 7 Wassereinfthrung; 8’ Dampfeinfiihrung; 9 Temperaturkontrollmess- 
stellen; 10 Temperaturfiihler fiir Kompressorsteuerung. 


Figur 2 zeigt den Temperierungs- und Filterteil einschliesslich der Befeuch- 
tungsstrecke etwas ausfithrlicher. 7 bedeutet hier das Geblase, 2 den Ammo- 
niak-Verdampfer oder Luftkihler, 3 ein Heizelement und 4 das elektrostatische 
Luftfilter. Durch die Leitung 5 werden die Eisbildungskerne eingespritzt, wah- 
renddem 6 und 7 fiir Wasser-, 8 fiir Dampfeinfiihrung bestimmt sind. Die mit 
9 bezeichneten Punkte bedeuten TemperaturkontrollmeBstellen, wahrenddem 
am Orte 70 ein Temperaturfiihler fiir ein Regelgerat liegt. Es sei noch vermerkt, 
daB Anlage und Labor nach dem Baukastenprinzip errichtet sind und leicht 
demontiert oder umgebaut werden kénnen. 


4.2 Zur Konstruktion des Windkanals 


Die Bedingungen an den Kanalquerschnitt wurden gegeben durch die For- 
derung eines MeBstellenquerschnittes von 25 x25 cm; dies fiithrte zu einer lich- 
ten Weite von 50 x50 cm, die nur im Temperierungsteil verdoppelt wurde, um 
durch Vergrésserung der Luftkithler-Ubergangsflachen sowie der Verweilzeiten 
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bessere Warme- bzw. Kalteaustauschwerte zu erreichen. Die Hohe des Kanal- 
systems wurde bedingt durch die notwendige Lange der Ausgleichsstrecke, soll 
doch das darin eingespritzte Wasser bis zur Mefstelle die Temperatur der um- 
gebenden Luft praktisch angenommen haben. 

Von der Kanalwandung wurde gefordert: gute Isolation, geringe Warme- 
kapazitat sowie durchgehend praktisch gleicher Ausdehnungskoeffizient. Die 
erste Bedingung resultierte aus dem Wunsch, auch bei tiefsten Betriebstempe- 
raturen von — 60°C nur geringe Warmeverluste decken zu miissen, die zweite 
aus dem Bediirfnis, rasche Temperaturanderungen vornehmen zu k6énnen, 
wahrenddem die dritte allzu grosse Spannung innerhalb der Kanalwand aus- 
schliessen soll. Nach verschiedenen Material- und Bauversuchen wurde auf 
10 cm dicke Zublanit-Kunststoffplatten gegriffen. Sie bewahrten sich in der 
Folge ausgezeichnet. Die mechanischen Eigenschaften liessen es zudem zu, die 
Platten zu einer selbsttragenden Konstruktion zu verleimen. Zum Schutz gegen — 
Beschadigungen sowie zur Schaffung einwandfreier Verbindungsstellen wurden 
allerdings die Kanten mit Holzleisten geschiitzt und aus Holz bestehende 
Flansche aufgeklebt. 

Zur Glattung der Oberflache sowie zu ihrem weiteren Schutze wurde der 
Kanal in- und auswendig mit Araldit-Spachtelharz 880 AB bearbeitet. Dieser 
Uberzug, der speziell fiir diesen Zweck von der CIBA entwickelt wurde, zeigte — 
auch bei den tiefsten Temperaturen noch eine gewisse Elastizitat, die der sonst 
allgemein zu beobachtenden Rissbildung entgegenwirkt. Allerdings scheint eine 
extrem starke Strapazierung mit standigen Temperaturwechseln und Ver- 
eisungen eine Durchweichung des Anstriches und eine Verminderung seiner 
Haftfahigkeit an den Platten zur Folge zu haben. Die Unterteilung des Kanals 
erfolgte meist in 1m lange Stiicke, die mit aufklappbaren Schrauben leicht 
getrennt und wieder zusammengesetzt werden kénnen. Zur Einfiihrung von — 
MeSsonden und Einspritzdiisen oder auch zu Beobachtungszwecken sind ver- 
schiedene Deckel und Klappen angebracht. 


4.3 Das Geblise 


Die Erzeugung der verlangten Luftgeschwindigkeit geschieht mittels eines 
Sulzer-Radialgeblases, das durch einen Gleichstrom-Nebenschlussmotor an-_ 
getrieben wird. Die variable Feldspannung wird tiber eine Tyratrongruppe pro- 
duziert. - Zur Vermeidung von Druckstéssen wurde eine Drehzahlprazision 
von besser als 1% iiber den ganzen Drehzahlbereich von 100 bis 3000 U/min 
gefordert, unabhangig von der Belastung des Motors. Diese Bedingung konnte 
mit einem elektronischen Regelverstarker eingehalten werden, der den Befehl 
des Prazisionspotentiometers (zur Sollwert-Drehzahleinstellung) vergleicht mit 
dem geblasegekuppelten Tachodynamo. Eine allfallige Abweichung zwischen 
Soll- und Istwert wird sofort auskorrigiert (Lieferfirma Sprecher & Schuh AG) 
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4.4 Die Kiailteanlage 


Der prinzipielle Aufbau des Hagelversuchskanals hat zur Folge, dass die 
geforderte Feuchte der Kanalluft stetig beigegeben werden muss. Im Extrem- 
fall ist daher anzunehmen, dass simtliches eingespritzte Wasser in fester oder 
fliissiger Phase wahrend eines Kreislaufes im Kanal oder am Luftkiihler nieder- 
geschlagen wird. Der Wasserdampf selber diirfte bis auf den Dampfdruck iiber 
dem kaltesten Teil kondensieren oder sublimieren. — Die aufzuwendende Kalte- 
energie ist aus diesem Grunde verhiltnismassig gross und steigt proportional 
mit dem Wassergehalt pro Kubikmeter Luft und der Luftgeschwindigkeit. — 
Der Vorteil der Methode liegt im einfachen Aufbau der Versuchsanlage sowie 
der grossen Variabilitat der Versuchsatmosphare. 

Nachstehend in Abschnitt 4.7 besprochene Eigenheiten der Tropfenerzeugung 
verlangen eine Heissdampfeinspritzung in den Kanal, die allein gewisse kleinste 
Tropfengréssen bei noch einigermassen verniinftigen Apparatebedingungen 
garantiert. Der Kaltebedarf wird aber damit gegeniiber einer direkten Wasser- 
einspritzung um einen Faktor von rund 7 vermehrt! Daher riibrt das Verlangen, 
mit der Kalteanlage Temperaturen von — 60°C zu erreichen, damit das resul- 
tierende Wasser-Luft-Gemisch eine tiefste Temperatur von — 40°C erreiche 
(gemass Bedingungen in Abschnitt 2.3). 

Auf Grund des dadurch vorgeschriebenen Kaltebedarfes — er wird spater 
explizit berechnet — plante die Firma Sulzer die Kalteanlage. Als Maschinentyp 
konnte sie einen 2stufigen Serientyp CV 180 (Figur 3) verwenden, der durch 
Aufbau einer speziellen dritten Stufe ein besonders auch bei tiefsten Tempera- 

-turen vorteilhaftes Leistungsdiagramm aufweist. Als Arbeitsmedium wurde 
Ammoniak gewahlt, das direkt im Luftkithler verdampft. Der Kaltekompressor 
wird durch einen 25-PS-Motor wahlweise auf 360 oder 720 U/min betrieben. 

Eine gewisse Komplizierung der Kalteanlage ergab sich durch den Standort 
und die damit verbundene Unméglichkeit einer Wasserkiihlung des Ammoniak- 
kondensors sowie des Kompressors selber. So wurde auf eine Aussenluftkiih- 

‘lung gegriffen, die sich auch bei leicht stiirmischem Wetter bewahrte (Figur 1, 
Positionen 9 und 70). 

Um eine allzu grosse Uberhitzung der Ammoniakdampfe in den einzelnen 
‘Verdichtungsstufen zu verhindern, wurden zusitzlich zwei Zwischenkthler 
-eingebaut, die nach Bedarf automatisch zum Einsatz kommen. Der Luftkiihler 
‘oder Verdampfer (Figur 2, Position 2) im eigentlichen Windkanal ist als 
_Rippenrohrsystem ausgebildet, durch das die Luft mit einer rund 8mal kleine- 

‘ren Geschwindigkeit als durch die MeBstelle streicht. 

Die Steuerung der Kalteleistung kann auf zwei Arten geschehen: die erste 
Art ist von Bedeutung, wenn tiefere Temperaturen auf schnellstem Wege er- 
‘reicht werden sollen. Sie besteht in einer Leistungssteuerung, die eine Einstel- 

dung eines Motorventils (Figur 1, Position 72) in der Ammoniaksaugleitung so 
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Figur 3 
Der Kaltekompressor. 


wahlt, dass die Kompressormotorleistung mit dem eingestellten Sollwert iiber- 
einstimmt. Insbesondere kann damit bei Maximalleistung durch das Tempera- 
turgebiet gefahren werden, bei dem der Kompressor mehr Energie aufnehmen 
kann, als der Motor abgibt. Dies kann der Fall sein bis zu Verdampfertempera- 
turen von rund —5° C bei Tourenzahl I und —30° C bei Tourenzahl II; diese 
Werte werden aber auch durch die Kondensationstemperatur des Ammoniaks 
beeinflusst, das heisst durch die Aussenlufttemperatur. 

Die zweite Steuerart erlaubt ein Erreichen und die Einhaltung einer ge- 
wahlten Verdampfertemperatur. Der Regelfiihler befindet sich dabei aller- 
dings nicht direkt am Verdampferelement, sondern im Luftstrom gemass Posi: 
tion 70 von Figur 2. Die Regelung geschieht wiederum tiber das gleiche Motor: 
ventil in der Ammoniaksaugleitung, das bei Erreichen der Luftkiihler-Sollwert 
temperatur die Kalteleistung des Kompressors der Kalteabgabe am Ver 
dampfer anpasst. Diese Steuermethode ist speziell von Vorteil bei der Ein 


haltung gleicher Ausgangsbedingungen einzelner Versuchsreihen. Als Regel 
system beliebte eine Rheoflex-Anlage der Firma Sauter AG. 


| 
: 
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4.5 Steuerung und Regelung der Mefstellentemperatur 


In der Natur muss damit gerechnet werden, dass Eispartikeln in einer Ge- 
witterwolke mit mehr oder weniger Aufwind laufend in kiltere oder wirmere 
Zonen gelangen. Die Versuchsbedingungen beziiglich der Temperatur miissen 
daher auch im Hagelversuchskanal relativ rasch variiert werden kénnen. Dieser 

_Bedingung wirkt aber die Warmekapazitat der rund 1,2 t schweren Ver- 
dampferelemente entgegen, die den Hagelversuchskanal in dieser Hinsicht sehr 
trage gestalten. - So wurde denn das energetisch etwas unsinnige Prinzip an- 
gewendet, am Verdampfer mehr Kilte bei tieferen Temperaturen zu erzeugen, 
als im jeweiligen Versuch gebraucht wird. Der Kalteiiberschuss wird mittels 
einer im Kanal eingebauten 19-kW-Heizung vernichtet (Figur 2, Position 3). 

Da die durch den Luftstrom des Kanals gefiihrten Heizdrahte praktisch 
tragheitslos eingesetzt werden kénnen, ist damit zu rechnen, dass die Tempe- 
ratur gentigend rasch variiert werden kann, wenn auch der Verdampfer seine 
Temperatur praktisch beibehalt. Da jedoch das Ausmass der Temperaturande- 
Tung wesentlich durch die Tragheit des Kanalstiickes zwischen Heizung und 
MeBstelle beeinflusst wird, kénnen nahere Daten nur durch entsprechende 
Messungen ermittelt werden (siehe Abschnitt 6.3). 

Die Heizleistung kann auf jeden Zwischenwert eingestellt werden, nachdem 
die Lieferfirma (Sprecher & Schuh), ahnlich wie bei der Ventilatorsteuerung, 
auf eine Gleichstromspeisung gegriffen hat. Je nach den Impulsen eines elektro- 
nischen Steuerdispositivs werden Schalttyratrons eingesetzt, die in der Folge 
drei Ignitronréhren ziinden. Die damit produzierten Gleichstréme werden 
tiberlagert und der Kanalheizung zugefiihrt. 

Der Einsatz der Heizung kann auf folgende Art geschehen: 

a) Handsteuerung. An einem Potentiometer kann die gewtinschte Heiz- 
leistung direkt in Kilowatt eingestellt werden. 

b) Temperaturregelung I. Wird eine bestimmte MeBstellentemperatur ge- 
wiinscht, so kann diese an einem Ferngeber eingestellt werden. Uber einen 
Temperaturfiihler sorgt ein elektronisches PID-Regelsystem (Fabrikat Hart- 
mann & Braun, Frankfurt) fiir raschestes Erreichen des Sollwertes durch ent- 
sprechende Befehlsabgabe an das Steuerdispositiv der Kanalheizung. Infolge 

der vorteilhaften Konstruktion des Reglers «Regelux 3» kann durch eine do- 

-sierte, separate Einstellung der Proportionalregelung, des Vorhaltes und der 

_Nachhaltezeit die optimalste Regelung entsprechend den Anlage- und Betriebs- 

_bedingungen erreicht werden. 

c) Temperaturregelung II. Anstatt auf die MeBstellentemperatur kann aut 
die Verdampfertemperatur geregelt werden, was speziell beim Abtauen des 
-Verdampfers, wo die Kalteanlage stillgelegt wird, von Vorteil ist. 

d) Temperatur-Programmsteuerung. Ein Programmgeber ersetzt den Fern- 

geber der Temperaturregelung I und gestattet damit ein zum voraus gewahltes 


392 ROLAND LIsT ZAMP 


Temperaturprogramm ablaufen zu lassen. (Die Planung des gesamten Tempe- 
raturregelteiles wurde in Zusammenarbeit mit der Firma Camille Bauer AG, 


Basel, ausgefiihrt.) 
4.6 Die Filterstrecke 


Die Reinigung der Luft, insbesondere von den eingespritzten Eisbildungs- 
kern-Substanzen, obliegt dem elektrostatischen Luftfilter, Fabrikat Trion AG, 
Ziirich. Dessen vorteilhafte Eigenschaften wurden bereits frither beschrieben 
[13], so dass auf eine Wiederholung verzichtet werden kann. 


4.7 Die Befeuchtungs- und Ausgleichsstrecke 


Die Befeuchtungs- und Ausgleichsstrecke iibernimmt trockene und gerei- 
nigte Luft einer bestimmten Temperatur aus dem Temperierungs- und Filter- 
teil und schafft durch Zugabe von Feuchtigkeit fliissigen oder dampfférmigen 
Charakters innerhalb der Ausgleichsstrecke das gewiinschte Klima. Die Lange 
der Ausgleichsstrecke wird, wie schon erwahnt wurde, bestimmt durch die Zeit, 
in der eine Anpassung der Temperaturen der erzeugten Wassertropfen an die 
der Umgebungsluft erfolgt. Eine theoretische Berechnung dieser Verhaltnisse 
zeigte in der Folge, dass eine Ausgleichsstrecke von 2 m Lange geniigt, sofern 
nur Tropfen eingespritzt werden. Praktische Versuche zeigten sogar auch 
stationare Mischtemperaturen bei Dampfeinspritzung und nachfolgender Kon- 
densation und Abkihlung. 

Die Produktion von Tropfengesamtheiten verschiedener Tropfendurch- 
messer und Konzentrationen stellt eines der schwierigsten Probleme der Erzeu- 
gung atmospharischer Verhaltnisse dar. Im Gegensatz zu den grossraumigen 
Gewitterwolken komplizieren sich die Verhaltnisse in einem Laboratoriums- 
Experimentierraum durch das Vorhandensein von Wanden. Meist ist es da 
nicht méglich, Wasser in fliissiger Phase durch Expansion nahezu gesattigter 
Luft zu erhalten. Die entsprechenden apparativen Bedingungen sind zu kompli- 
ziert, auch spielt der Zeitfaktor eine ungiinstige Rolle, miissten wir doch die 
Kondensation auf kurzen Wegstiicken produzieren. Daher ist es nicht verwun- 
derlich, dass auf verschiedene andere denn beim Wettergeschehen vorkommende 
Tropfenproduktionsvarianten gegriffen werden musste. 

Als Hauptbedingungen an die Erzeugungsmethodik sind zu erwahnen: 

1, Keine merkliche Beimischung kanalfremder Luft. 

2. Die eingespritzten Tropfen diirfen keine allzu hohe Geschwindigkeit 
haben; es soll der Impuls der eingespritzten Feuchte nur ein Bruchteil des Im- 
pulses der dazugehérigen Luftmenge sein. (Damit sollten die Stérungen, die 
sich auf die Strémung iibertragen, auf ein ertragliches Mass reduziert werden.) 

Diese beiden Einschrénkungen liessen verschiedene sonst gebrauchliche 
Methoden als unbrauchbar erscheinen, so zum Beispiel jegliche Tropfenerzeu- 
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gung mit Beimischung von Pressluft an der Diise»Die zur Produktion kleinster 
Wasserpartikeln ohne Pressluft notwendigen hohen Drucke (einige 100 kg/cm?) 
mussten ebenfalls fallengelassen werden. So ergaben sich schliesslich noch fol- 
gende brauchbaren Erzeugungsarten: 

a) Fiir Tropfen mit Durchmessern der Gréssenordnung 1 uw: Einfiihrung von 
kondensierendem Heissdampf in die Kaltluft des Kanals. 

b) Fur Tropfen von 5 bis 20 w: Einfiihrung eines in einem mechanischen Zer- 
stauber erzeugten Tropfen-Luft-Gemisches (die beigemischte Luftmenge ist hier 
relativ klein). 

c) Fir Tropfen mit Durchmessern von 20 "und mehr: Einspritzung von 
Wasser durch konische Schleuderdiisen, 

Die entsprechenden Apparate sind: 

ad a) Dampfkessel: Der Heissdampf wird in einem Elektrodendampfkessel 
mit einer Stundenleistung von 20 bis 40 kg bei einstellbarem, konstantem 
Druck erzeugt. Die Beimischung zur Kanalluft erfolgt iiber ein fernsteuerbares 
Motorventil, das geeichte Dampfmengen ausstrémen lasst. 

ad b) Zerstduber: Eine energetisch interessante Produktion kleiner Tropfen 
kann iiber eine rotierende Schleuderscheibe geschehen. Durch stetige Wasser- 
zufiihrung wird auf dieser Scheibe ein abgleitender Wasserfilm erzeugt, der 
laufend auf ein umrahmendes Kammsystem aufprallt. Die Wasserhaut zer- 
platzt darauf zu nicht benetzenden Trépfchen, die in hohen Konzentrationen 
von einer geringen, automatisch mit angesogenen Luftmenge mitgenommen 
und der zu befeuchtenden Kanalluft beigemischt werden. Der Wassernachschub 
geschieht vollautomatisch iiber ein schwimmergesteuertes Magnetventil. 

adc) Schleuderdiisen: Durch einen Luftkompressor, dessen Vorratsgefass 
zum Teil mit Wasser gefiillt ist, kann letzteres tiber ein Magnetventil in eine 
Schleuderdiise gefiihrt werden. Das Wasser wird in dieser tangential eingeleitet 
und zerspritzt direkt in die Kanalluft. Durch Variation des aufgewendeten 


_ Druckes kann die Wassermenge dosiert werden. Fiir speziell konstante Bedin- 


gungen wird anstelle des Kompressors eine Pressluftflasche tiber ein Reduzier- 
ventil auf das Luftkissen im Wasservorratsgefass geschaltet. — Die Einspritz- 
diise ist zudem heizbar und kann ebenfalls fernbedient werden. 

Der Aufbau dieser Anlagen geht aus Figur 4 hervor, die als Grundriss die 


_ Anordnung von Dampfkessel 7, Motorventil 2, Zerstauber 3, Luftkompressor 4 
und Befeuchtungsstrecke 5 darstellt. 


Der Charakter der kiinstlich erzeugten Atmosphare wird aber ausser durch 


t Temperatur und Feuchte (und dem umgebenden, aber nicht variablen Luft- 
_ druck) weiter bestimmt durch die Art und Konzentration der Luftverunreini- 
- gungen, die speziell fiir die Einsatztemperatur des Gefrierprozesses verantwort- 
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lich sind, das heisst durch die Eisbildungskerne. Die Erzeugung derselben erfolgt 


nach der bereits friiher beschriebenen Methode [13] durch Verdampfen spe- 


zieller Substanzen in einem Hochtemperatur-Rohrofen (Figur 4, Position 6). 
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Figur 4 
Die Aufstellung der Pevenero ee (Ansicht des untersten Stockwerkes von oben) 
(1:90). 1 Dampfkessel; 2 Dampfmotorventil; 3 Zerstauber; 4 Luftkompressor; 5 Befeuchtungs- 


strecke; 6 Hocktesieratinoten: 7 Ventilatormotor; 8 Temperierungs- und Filterteil; 9 Flussig- 
keitsabscheider; 10 Kaltekompressor; 11 Kompressormotor; 12 elektrische Hauptverteilung; 
13 Wasserzuleitungen. 


Eine Sublimation dieser Dampfe in durchstrémender Luft ergibt Teilchen mit 
Durchmessern von 0,02 bis 0,05 p, die der Kanalluft beigegeben werden 
konnen. 

Uber die Zusammenhange zwischen Substanz und Gefriereinsatztemperatur 
sei zum Beispiel auf die Veréffentlichungen von PRUPPACHER und SANGER 
{[17, 18] verwiesen. 

Aus ali diesen Ausfiithrungen geht hervor, dass die MeBstellentemperatut 
nichts anderes darstellt als die Mischtemperatur, die sich ergibt, wenn der 
trockenen, kalten Luft des Temperierungsteiles Feuchte irgendwelcher Form 
beigegeben wird. Aber auch die Eisbildungskerne kénnen diese resultierende 
Temperatur beeinflussen, sofern sie einen Phaseniibergang fliissig—fest ganz 
oder teilweise vor der MeBstelle zur Folge haben. Die Untersuchung dieses 
letzten Falles ist besonders interessant fiir einen Vergleich des Wachstums eines 


Objektes in einem Aufwind mit angestrémten kleinern Partikeln fliissiger 
fester oder gemischter Phase. 


4.8 Die Mefstelle 


Der Raum, in dem die eigentlichen Experimente stattfinden, wird als MeB 
stelle (Figur 5) bezeichnet. In ihr wird die Vereisung von aufgehangten oder fre 
im Luftstrom ausgeschwebten Objekten studiert. Im Gegensatz zu den andert 
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Kanalstiicken, wo gleiche Elemente Luftstrombegrenzung und Isolation iiber- 
nehmen, sind diese Funktionen bei der MeBstelle getrennt. Die direkte Begren- 
zung des Luftstromes ist, um diese Vorgéinge sichtbar zu machen, aus Plexiglas 
7 gehalten. Ein allfalliges Vereisen desselben kann durch Einschalten und 


50 cm 


Figur 5 


Schnitt durch die MeB stelle (1:15). 1 Plexiglasrahmen; 2 Strahlungsheizung; 3 Lautsprecher;, 
4 Mikrophon; 5 Feuchtesonde; 6 Temperaturfiihler; 7 MeBstellenrahmen; 8 auswechselbare Abdeck- 
platten; 9 Fenster; 10 Durehfithrungen. 


Leistungsregulierung der angebrachten zwei Strahlungsheizungsrohren 2 ver- 
_hindert werden. Prinzipiell dienen zwei Gegenseiten zur Beobachtung, wahrend- 
dem die gekreuzten Gegenseiten verschiedene MeBsonden aufzunehmen haben. 
Ein dusserer Rahmen, unterteilt in einzelne, vertauschbare und teilweise 
“mit Glasfenstern versehene Platten reduziert die Warmeverluste. Dies hat den 
Vorteil, dass die MeBsonden nur durch das diinne Plexiglas eingefiihrt werden 
“miissen, wobei zwischen diesem und der Isolation noch geniigend Raum zum 
_Auf- und Ausbau derselben besteht. 
_ Zur Erleichterung einer Ausmessung des Luftstromes innerhalb des geraden 
'Plexiglaskanals sind auf zwei einander gegentiberliegenden Platten zwei Raster- 
“netze eingeritzt, die eine rasche Einstellung der Mefisonden auf die bezeichneten 


, 
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Positionen erméglichen. Zur Einfiihrung der Sonden sind auf den gekreuzten 
Gegenseiten Schlitze im Abstande von 3,5 cm angebracht, so dass im ganzen 
sieben Ebenen parallel zu den Rasterplatten erfasst werden konnen. Ein derart 
abgetasteter Luftstrom diirfte innerhalb der MeBstelle als gentigend bekannt 
gelten. Das gerade Plexiglassttick der Luftstromfiihrung — das Kernstiick der 
MeBstelle — ist auf Schienen senkrecht zum Luftstrom beweglich. Dies ist ein 
grosser Vorteil, wenn darin aufgehangte, vereiste Objekte aus dem Kanalwind 
entfernt werden mussen. 


4.9 Die Strimungsverhalimsse 


Bei der Konstruktion des Windkanals wurde grosser Wert auf méglichst 
einwandfreie aerodynamische Gestaltung gelegt. Dabei ist speziell zu bemerken, 
dass eine Vernichtung grésserer Turbulenzen, herriihrend vom Radialgeblase, 
bereits im folgenden Kanalstiick mit dem Luftkiihler erfolgt. Auch das aus einer 
Reihe paralleler Platten bestehende elektronische Luftfilter wirkt sich giinstig 
auf die Gleichrichtung des Luftstromes aus. — Die Umlenkung in die Aus- 
gleichsstrecke erfolgt unter Kontraktion durch vier Umlenkschaufeln, die eine 
gleichmassige Geschwindigkeitsverteilung im Vertikalteil gewahrleisten. Vor 
der MeBstelle wird der Luftstrom nochmals in einem Effusor beschleunigt. 

Eine Ausmessung der Geschwindigkeitsverteilung innerhalb der MeBstrecke 
ergab die Tatsache, dass, abgesehen von einer unwichtigen Randzone von 
2,5 cm Machtigkeit, die Abweichung von der mittleren Geschwindigkeit im 
Bereich von 1 bis 10 m/s 1% nicht wberstieg; im Bereich von 10 bis 28 m/s 
betrug die maximale Schwankung 1,5%%. 

Der auf die MeBstelle folgende Diffusor wurde mit einem totalen Offnungs- 
winkel von 6° konstruiert und biirgt damit fiir eine méglichst vollstandige 
Riickgewinnung der Druckenergie. 

Was die grésste erreichbare Luftgeschwindigkeit anbelangt, so kann er- 
wahnt werden, dass diese ohne weiteres durch Verengung des MeBstellenquer- 
schnittes erhoht werden kénnte. Zur Gleichmassigkeit der Luftgeschwindigkeit 
ist zu vermerken, dass diese nur gewiinscht wird bei Experimenten mit fixierter 
Objekten. Frei ausgeschwebte Partikeln verlangen andere Verhiltnisse mit 
ganz bestimmten Geschwindigkeits- und Druckverteilungen tiber den Quer: 
schnitt und die Lange der MeBstrecke [11], wie sie bis jetzt nur mit Einbauter 
von Gittern und Stauhauben erreicht werden konnten. Die einzige Ausnahme 
bilden die Versuche, wie sie am Aerodynamischen Institut der ETH mit Kurz 


diffusoren unternommen worden sind [10]. Doch soll an dieser Stelle nicht weitel 
auf diese Punkte eingegangen werden. 


4.10 Das Steuerzentrum 


Das Nervenzentrum des Hagelversuchskanals ist das Steuerpult (Figur 6) 
In ihm laufen alle Steuer- und Messleitungen zusammen. Der Aufbau geschal 
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Figur 6 
Das Steuerzentrum. 


auch nach dem Baukastenprinzip, so dass Umbauten oder Erganzungen der 
Verdrahtungen mit Leichtigkeit geschehen kénnen. Die Messungen der verschie- 
denen Gréssen erfolgt mit Ausnahme der Strommessungen iiber Réhrenvolt- 
meter-Schaltungen. Zehn solche Einheiten mit teilweise separater Basis und 
teilweise getrennter Speisung sind in einem Kasten untergebracht, wo die zu 
messenden Gleichspannungen angeschlossen werden kénnen. Wechselspannun- 
gen werden zuerst iiber Steckeinheiten gleichgerichtet. Diese Messgerate kénnen 
an alle Spannungen zwischen 0,5 und 500 V so angepasst werden, dass jeder 
Zwischenwert auf den Anzeigeinstrumenten des Pultes als Maximalausschlag 
von 100% registriert werden kann. 

Das Steuerpult ist fahrbar und wird iiber 16 Zehnfachkabel mit den tibrigen 
Apparaten verbunden. 

Die Anordnung des Messbodens ist aus Figur 7 ersichtlich, die die Lage des 
Steuerpultes 7 direkt vor der MeBstelle darstellt. Uber dem Steuerpult in 
Augenhohe liegen auch die drei Kompensationsschreiber 2, die die Protokoll- 
fiihrung tibernehmen. Der eine dient zur Registrierung der MeBstellentempera- 


398 RovanpD List ZAMP | 


6 


RAQAAVAAAY 
oa 


| 
| 
| 
| 
| 
7 
| 
I 
| 
| 
| 
| 
| 
| 
| 


Figur 7 
Apparateanordnung auf dem Messboden (1:50). 1 Steuerpult; 2 Kompensationsschreiber; 3 Ver- 
teilerkasten mit Temperaturregelung; 4 elektronische Umformergruppe zu Ventilator; 5 elektro- 
nische Umformergruppe zu Heizung; 6 Hangebahn; 7 Temperatur-Schallmessanlage; 8 Feuchte- 
und Eisbildungskernmessung; 9 elektrische Druckwandler; 10 Rost; 11 MeBstelle. 


tur; er ist ausgefiihrt als Linienschreiber. Die beiden andern Registriergerate 
sind 12-Punkt-Schreiber (Fabrikat Hartmann & Braun, Frankfurt), wobei der 
eine Temperaturen druckt, der andere in Gleichspannungen von 0 bis 10 mV 
umgewandelte elektrische Gréssen. Die Umwandlung dieser Gréssen und ihre 
Reduktion erfolgt im Verteilerkasten rechts, wo auch die Gerate zur Tempera- 
turregelung und Programmgebung eingebaut sind. Die beiden Geratekasten 4 
und 5 beherbergen die Umformergruppen fiir Ventilator- und Heizungssteue- 
rung, wahrenddem die gestrichelten Linien eine Hangebahn 6 andeuten, an det 
Gestelle mit der Temperaturmessung iiber die Schallgeschwindigkeit 7 (siehe 
Abschnitt 7.2), der Feuchte- und Eisbildungskernmessung 8 hangen. 9 weist 
auf die elektrischen Druckwandler zur Luftgeschwindigkeitsmessung hin. 
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Als weitere Erleichterung ist im Messzentrums-Boden ein Rost 70 ein- 
gelassen, der die Uberwachung des Kompressors gestattet; tiber ein Spiegel- 
system lassen sich dessen Arbeitspunkte anhand von Manometern ablesen. 

Diese Anordnung und dieser Aufwand hat zur Folge, dass einfachere Mes- 
sungen gut von einem einzigen Mann ausgefiihrt werden kénnen. Fiir kompli- 
ziertere Experimente ist eine zweite Kraft einzusetzen. 


5. Die Leistung des Hagelversuchskanals bei stationairem Betrieb 


Die Bedingungen, die an den Hagelversuchskanal gestellt wurden (siehe 
Abschnitt 2.3), betreffen eigentlich nur einzelne Arbeitspunkte. Es diirfte daher 
interessant sein zu verfolgen, erstens wie sie eingehalten werden kénnen und 
zweitens, wie die Leistungsfahigkeit in den Zwischengebieten verlauft. Die 
Kenntnisse der Leistungsfahigkeit sind zudem von ausserordentlicher Bedeu- 
tung fiir die Gestaltung des Versuchsprogrammes itiber die Wachstumsbedin- 
gungen von Eispartikeln. Es wird auch auf Grund dieser Messungen méglich 
sein, abzuschatzen, ob und wie weit extremere Verhaltnisse erzeugt werden 
k6énnen, als sie bei der Bildung von schwerem Hagel wahrscheinlich sind. 

Die Leistungsmessungen zerfallen prinzipiell in zwei Gruppen: 

1. Das Verhalten des Hagelversuchskanals bei stationaéren Betriebsbedin- 
gungen. 

2. Das Verhalten des Hagelversuchskanals bei instationaren Betriebsbe- 
dingungen. 

In diesem Abschnitt soll nur von der ersten Gruppe gesprochen werden, 
wahrenddem die zweite dem folgenden Abschnitt Inhalt gibt. 

Wenn wir eine allgemeinere Energiegleichung fiir den stationdren Fall auf- 
schreiben, so lautet sie prinzipiell wie folgt: 


= Or + Qn + OQ - (5) 

Dabei sind: 

_K KaAltemenge pro Stunde, zugefiihrt durch den Kompressor; 

QO, Leistung, zugefiihrt durch die eingespritzte Feuchte ; 

Qy stiindliche Warmemenge, zugefiihrt durch die Heizung; 

Q, Leistungsverluste durch Warmezuleitung durch die Kanalwand zusatzlich 
der vom Geblise und den Wirbelverlusten erzeugten Warme. 

- Die Betrachtung der einzelnen Gréssen zeigt: 

a) Die Leistung des Kompressors K ist nur theoretisch bekannt und kann 

leider nicht einfach gemessen werden. Sie ist zudem abhangig von der Konden- 

sationstemperatur des Ammoniaks, die ihrerseits infolge der Luftktihlung von 

- der Aussenlufttemperatur abhangig ist. Eine Rolle spielt ferner der Betriebs- 

_ gustand, das heisst angesammelte Olmengen in den Zwischenkihlern der Kom- 

_ pressorstufen usw. 
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b) Die durch Einspritzen von Feuchtigkeit zugefiihrte Leistung Q; kann 


genau errechnet werden. 
Die entsprechende Formel fiir Wassereinspritzung lautet: 


On = ws [by — ty) Cw + S] Um Fy - (6) 
Die Dampfeinspritzung ergibt sich zu: 
Or = Wp [tp — te) Ow + EL + S] ng Py. iG 


Dabei sind: 
w, Menge des eingespritzten fllissigen Wassers {kg/m3]; 
Wp Menge des eingespritzten Dampfes [kg/m*] ; 
t,; Temperatur des eingespritzten Wassers [°C]; 
ty Temperatur des eingefiihrten Dampfes [°C]; 
ty Temperatur an der MeBstelle [°C]; 
C» spezifische Warme von Wasser [kcal/kg °]; 
S  Gefrierwarme des Wassers [kcal/kg]; 
L Verfliissigungswarme des Wasserdampfes [kcal/kg]; 
vy Luftgeschwindigkeit an der MeBstelle [m/h] ; 
Fy, Querschnittsflache der MeBstrecke {m7}. 


Die mit Formel (6) ermittelten Daten sind fiir Luftgeschwindigkeiten von 
10 und 25 m/s in Figur 8 eingetragen. Als Temperatur des eingespritzten 
Wassers wurde 10°C eingesetzt, die MeBstellentemperatur wurde mit 0°, — 20° 
und -- 40°C gewahlt (Jy, = 0,0676 m?). Analoge Werte fiir Dampf zeigt gemdss 
Formel (7) Figur 9, wobei als Dampftemperatur 100°C gilt. 

Ein Vergleich dieser beiden Darstellungen weist darauf hin, dass bei Was- 
sereinspritzung die MeBstellentemperatur einen relativ grossen Einfluss hat, 
wahrenddem dieser bei Dampfeinspritzung nahezu verschwindet. 

c) Die durch Heizung zugefiihrte Warmemenge Q,, kann durch Messung 
des Stromes und der Spannung ermittelt werden. 

d) Die Verluste Qy zerfallen in zwei Teile. Der erste umfasst die Warmeent- 
wicklung durch den Ventilator — sie sind nicht allzu gross, da sich der Antriebs- 
motor ausserhalb des Windkanals befindet — sowie die Warmeentwicklung der 
Luftstr6mung an und fiir sich. - Dem Kanal wird auch Warme von aussen zu- 
gefithrt durch Warmeleitung durch die Kanalwande hindurch. Dieser Effekt 
ist abhangig von der Isolationsfahigkeit der Kanalwandung sowie auch von der 
Luftgeschwindigkeit, die die Warmetibergangszahl Kanalwand—Kanalluft be- 
stimmt. Eine Hauptrolle spielt hier natiirlich die antreibende Temperaturdiffe- 
renz zwischen Labor- und Kanalluft. 

Aus diesen Betrachtungen ist zu ersehen, dass eine allgemeine, hinreichend 
genaue Erfassung der Energiebilanz des Kanals viel zu kompliziert ware. Wir 
haben jedoch die Méglichkeit, durch giinstig ausgelegte Messreihen auf einfache 
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Figur 8 


Dem Kanal zugefiihrte Leistung Qp als Funktion der Wassereinspritzung w,;, mit den Parametern 
Luftgeschwindigkeit und Me8stellentemperatur. 


0 1 2 8 4 g/m 
Figur 9 
Dem Kanal zugefiihrte Leistung Qf als Funktion der Dampfeinspritzung wp, mit den Parametern 
Luftgeschwindigkeit und MeBstellentemperatur. 


Weise ein zuverlassiges Bild der Leistungsfahigkeit des Hagelversuchskanals zu 
_bekommen. Da es fiir die Kalteanlage gleich ist, ob Warme iiber Feuchtigkeit 
_ oder Heizung — nur diese sind auch im allgemeinen Experiment von Interesse — 

zugefiihrt wird, kénnen diese Gréssen ohne weiteres vertauscht werden. Es ist 

hierbei in erster Linie an die Ersetzung der Feuchteeinspritzung durch die ein- 
_ facher zu handhabende Heizung gedacht. Die Kanalleistungsfahigkeit kann nun 
_ direkt gemessen werden, wenn bei einer gegebenen Luftgeschwindigkeit an der 
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MeBstelle die Heizung einer bekannten Leistung eingeschaltet wird und gleich- 
zeitig der Kaltekompressor auf Volleistung lauft. Jetzt braucht man nur ein 
Gleichgewicht mit konstanten Temperaturen abzuwarten, wo sich die Kalte- 
leistung mit der Heizleistung und den Verlusten die Waage halt. Auf diese 
Weise kann unter Umgehung einer Messung und Abschatzung der Kanalver- 
luste direkt die dem allgemeinen Experiment bei einer bestimmten Temperatur 
und einer gegebenen Luftgeschwindigkeit zur Verfiigung stehende maximale 
grésste Leistung gemessen werden. Wird nun die Heizleistung verandert, so 
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Figur 10 


Grésste Kompressorleistung Lx bei verschiedenen Luftgeschwindigkeiten vag und den beiden Touren- 
zahlen KTZ I und KTZ II als Funktion der Verdampfertemperatur ty. 


ergeben sich andere Leistungswerte bei andern resultierenden Temperaturen. 
Figur 10 zeigt so ermittelte Daten bei der hohen und der niedrigen Kompressor- 
tourenzahl fiir Luftgeschwindigkeiten von 10, 15, 20 und 25 m/s. Die Verschie- 
bung der einzelnen Leistungskurven resultiert aus den unterschiedlichen Ver- 
lusten. 

Es muss hier noch vermerkt werden, dass die fiir die Kalteleistung mass- 
gebende Temperatur immer die Ammoniak-Verdampfertemperatur und nicht 
diejenige der MeBstelle ist ! Wir méchten jedoch ein Bild, das zeigt, wie es sich 
mit der letzteren verhalt. 

Die Beziehungen zwischen der Temperatur an der MeBstelle und am Ver- 
dampfer kénnen am einfachsten tiber die zwischen dem Verdampfer und der 
MeBstelle zugefiihrte Energie erfasst werden, wenn wir diese aufteilen pro durch- 
gestromten Kubikmeter Luft und mit der zugehérigen spezifischen Warme und 


einer mittleren Luftdichte dividieren. So resultiert die Temperaturerhohung 
auf Basis der Verdampfertemperatur. 


Die entsprechende Forme! lautet: 


YOKE er Fyumc, * (8) 
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Die neu eingefiihrten Symbole bedeuten: 
ty Yemperatur des Verdampfers [°C}; 

Q; bis zur MeBstelle durch Feuchte zugefiihrte Leistung [kcal/h]; 
Qu Leistung, zugefithrt durch Heizung [kcal/h]; 

Qy Verluste zwischen Verdampfer und MeBstelle [kcal/h); 

¢, spezifische Warme der Luft [kcal/kg°]; 

o, mittlere Luftdichte {kg/m}. 

Die bis zur Mefstelle durch Feuchte zugefithrte Energie Qj}. ist nicht mit der 
totalen Energie Q;, gemiss Gleichung (6) oder (7), gleichzusetzen. Normaler- 
weise sind ja die Wassertropfen bis zum Testobjekt nicht vereist. — Da jedoch 
ein Teil des eingefiihrten Wassers bereits in der Ausgleichsstrecke an den 
Wanden abgelagert und ausgefroren wird, kann trotzdem naherungsweise Q}, 
durch Q; (Figur 8) ersetzt werden. In Wirklichkeit sind aber meist giinstigere 
Werte zu erwarten. 

Werden der aktuelle MeBstellenquerschnitt von 0,0676 m2 und eine mittlere 
Luftdichte von 1 kg/m (¢ = — 20°C, £,,,. = 550 mm Hg) in Rechnung gesetzt, 
so ergeben sich Darstellungen der Funktion ¢,, gemass den Figuren 11 und 12. 
Die Figur 11 zugrunde gelegte Luftgeschwindigkeit betragt 10 m/s, wahrend- 
dem Figur 12 fiir 25 m/s gerechnet ist. Die Verluste Qj, sind direkt aus der 
Messung ermittelt worden bei 0; = Qy = 0, ebenso die Grenzleistung des Kom- 
pressors bei Tourenzahl II, die aus Figur 10 tibernommen worden ist. Das mit 
der Anlage durchfahrbare Gebiet befindet sich links der Grenzleistungskurve. 

‘Der Vergleich von Figur 11 mit Figur 12 zeigt, dass die Differenz von Ver- 
dampfer- und Me8stellentemperatur bei gleicher zugefiihrter Energie mit klei- 
nerer Luftgeschwindigkeit grésser wird. Dies ist verstandlich, wird doch eine 
kleinere Luftmenge durch eine gleichbleibende Warmezufuhr starker erwarmt 
als ein grésseres Luftvolumen. 

Es muss hier noch auf zwei Punkte hingewiesen werden, die auch fiir das 
_ Folgende Giiltigkeit haben: 

a) Die Grenzleistungskurve kann sich je nach Witterungsverhaltnissen ver- 
schieben. Da aber bei einer mittleren Aussenlufttemperatur von + 2°C gemes- 
sen wurde, diirfte es sich um eine iiber das Jahr hindurch mittlere Lage handeln. 
Im Winter sind durchschnittlich gréssere Leistungen zu erwarten als im Sommer. 

b) Der Fithler zur Ermittlung der Verdampfertemperatur ist am Rippen- 
rohrsystem ausserhalb der Kanalstrémung befestigt. Da die Temperaturver- 
teilung am Luftkithler inhomogen ist und sich je nach Betriebsbedingungen 
etwas andern kann, braucht der MeSwert nicht unbedingt mit dem effektiven 
- Wert iibereinzustimmen. Bei allen folgenden Darstellungen wird aber der ge- 
messene Wert als Grundlage genommen. 

Eine der wichtigsten Fragen, die bei stationdrem Betrieb an den Hagelver- 
_ suchskanal gestellt werden, geht nach dem Ausmass der grésstméglichen Was- 
sereinspritzung. Als Bezugstemperatur ist dabei die MeBstellentemperatur er- 
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Figur 11 


Abhangigkeit der MeBstellentemperatur ty von der Verdampfertemperatur ty als Funktion der — 
zugefiihrten Leistung Q77+ Op, zugrunde gelegte Luftgeschwindigkeit vjg = 10 m/s. Die Kurve Lxg 
stellt die Kompressorgrenzleistung bei Tourenzahl II dar. 


0 -10 -20 -30 iii) o 


Figur 12 


pati neces der MeBstellentemperatur ty von der Verdampfertemperatur ty als Funktion der 
zugefihrten Leistung Qy + Op, zugrunde gelegte Luftgeschwindigkeit vj = 25 m/s. Die Kurve Lx 
stellt die Kompressorgrenzleistung bei Tourenzahl II dar. ‘ 
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wiinscht. Diese Beziehungen kénnen anhand der Grenzleistungskurven. des 
Kompressors gemass Figur 10 sowie der Formeln (6) und (8) fiir Wasserein- 
spritzung errechnet werden. Es resultiert damit bei der Bedingung Or = Of 
eine Darstellung gemiass Figur 13. 

Umrechnungen auf Dampfeinspritzung kénnen iiber Figuren 8 und 9 ge- 
schehen oder zur Ermittlung von Richtwerten durch Division der Wassermen- 
gen gemass Figur 13 durch den Faktor sieben. 


We 
gim? 


20F 


0 -10 -20 -30 -40 -50 °C 
Figur 13 
Groésstmogliche Wassereinspritzung w,; bei verschiedenen Luftgeschwindigkeiten vg als Funktion 
der MeBstellentemperatur ¢y7 (Kaltekompressortourenzahl IT). 


Aus all diesen Daten ist ersichtlich, dass den Anforderungen an den Hagel- 
versuchskanal bei stationarem Betrieb vollauf Geniige geleistet wird. Im 
Gegenteil kénnte sogar die Frage aufgeworfen werden, warum die Kalteanlage 
speziell um 0°C eine so grosse, nicht geforderte Leistungsfahigkeit aufweist. 
Der Hauptgrund zu dieser Tatsache liegt darin, dass der Kompressor ausgelegt 
wurde, um auch Temperaturen von — 50°C und tiefer zu erreichen. Da seine 
Leistungskurve jedoch gegen tiefere Temperaturen stark abfallt, resultierte 
schliesslich eine entsprechende Mehrleistung bei relativ hohen Temperaturen. 
Dies ist jedoch sehr niitzlich, wenn rasch zu tieferen Experimentausgangstem- 
peraturen gefahren werden soll, auch erméglicht es, wie der folgende Abschnitt 
zeigen wird, gréssere Méglichkeiten beim dynamischen Ausfahren des Hagel- 
versuchskanals. 


6. Das dynamische Verhalten des Hagelversuchskanals 
6.1 Allgemeines 


| Im Unterschied zum stationaren Verhalten des Hagelversuchskanals spielt 
_ beim Betrieb mit wechselnden Verhiltnissen die Warmekapazitat der verschie- 
denen Apparateteile eine entscheidende Rolle. Die Leistungsbilanz gemass 
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Gleichung (5) ware damit insofern abzuandern, als die Verluste durch Warme- 
zuleitung Q, ihrem Sinne nach erweitert werden mussen und auch die stiind- 
liche Energie, die gebraucht wird, die verschiedenen Kanalteile in ihrer Tem- 
peratur an die Lufttemperatur anzugleichen, zu umfassen hat. Das Vorzeichen 
dieses zusatzlichen Betrages hangt jeweils von den Betriebsbedingungen ab. 
Die Hauptursache fiir die Warmetragheit des Hagelversuchskanals bildet der 
Verdampfer ~ er besteht aus rund 1,2t Eisen — wahrenddem Heizung und 
Filter inklusive Ubergangsstiicken aus etwa 50 kg Aluminium bestehen. Die 
eigentliche Kanalisolation aus Zublanitplatten dirfte infolge der geringen 
Dichte und kleinen Warmekapazitét keinen bedeutenden Einfluss auf die 
Tragheit der Temperaturanpassung haben. Von grésserer Bedeutung ist die 
direkte Kanalbegrenzung gegen innen, der Aralditharzanstrich. 


6.2 Der Einfluss der Kanaltraghett bei stetiger, gleichsinniger Temperatur- 
dnderung 


Die Warmetriagheit des Windkanals wird zuerst sichtbar, wenn der Kalte- 
kompressor eingeschaltet wird, um zu tieferen Kanaltemperaturen zu fahren. 
Die dabei zu beobachtende Abhangigkeit der Verdampfertemperatur von der 
Abkiihlungszeit ist fiir Luftruhe v,,=0 sowie fiir die Geschwindigkeiten v,, = 10 
bzw. 25 m/s in Figur 14 aufgetragen (durchgezogene Linien). Durch das Grésser- 


0 25 50 75 700 min 
Figur 14 
Abkuthlungs- und Aufheizkurven, gefahren bei héchster Leistung und verschiedenen Luftgeschwin- 
digkeiten, 


werden der Verluste bei ansteigender Luftgeschwindigkeit wachst auch die | 
Zeitspanne, in der eine bestimmte Verdampfertemperatur erreicht werden kann. 
(Diese Versuche wurden alle bei grésster Leistung auf Kompressortourenzahl II 
durchgefiihrt.) 


Nach Erreichen einer Verdampfertemperatur von — 50° bzw — 47°C wurde 
der Kompressor jeweils ausgeschaltet, um bei einem Heizeinsatz von 18 kW 
wieder gegen den Ausgangspunkt zuriickzufahren. Die Darstellung dieser Ver- | 
haltnisse erfolgt ebenfalls in Fi 14, j i i inl | 

n Figur 14, jedoch durch die gestrichelten Linien. — 
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Im Gegensatz zur Abkiihlung wird die Aufheizung durch eine grosse Geschwin- 
digkeit begiinstigt, da die mit ihr ansteigenden Kalteverluste im gleichen Sinne 
wie die Heizung wirken. 

Eine Differenzierung gewisser Kurven von Figur 14 gibt uns ebenfalls fiir 
zwei Luftgeschwindigkeiten den eigentlichen Abkithlungsgrad gemass Figur 15. 


Aty 
AT 
%min 
2 


Figur 15 


Abkuhlungsgrade in Funktion der Verdampfertemperatur ty: a bei grésster stetiger Temperatur- 

senkung fiir eine Luftgeschwindigkeit vyg = 10 m/s; b bei grésster stetiger Temperatursenkung fiir 

eine Luftgeschwindigkeit vjyg = 25 m/s; c maximale Abktthlung, ausgehend von stationaren Be- 
triebszustanden (praktisch geschwindigkeitsunabhangig). 


Das Hauptkennzeichen dieser Experimente liegt darin, dass hier die Verz6- 
gerung der Temperaturanpassung aller Kanalteile am extremsten ist! 

Arbeitet jedoch die Anlage seit langerer Zeit bei der gleichen Verdampfer- 
temperatur und konstanten Betriebsverhaltnissen, so ist anzunehmen, dass sich 
die Temperaturen der iibrigen Kanalteile mit der Zeit ebenfalls auf konstante 
Werte einstellen. Wollen wir nun von einem solchen stationaren Betriebspunkt 
mit angepassten Temperaturen mit Volleistung nach tieferen Temperaturen 
fahren, so waren eigentlich gréssere Abkiihlungsgrade zu erwarten als im weiter 
oben beschriebenen Falle, wo die Anpassung verzdgert stattfindet. Infolge der 
Anpassungszeit des Kompressors an neue Betriebsbedingungen ist dies nicht 
der Fall; es ergeben sich fiir diesen Fall praktisch geschwindigkeitsunabhangige 
Daten, die im Mittel durch die gestrichelte Kurve von Figur 15 dargestellt 
~ werden. 
Wird beim Aufheizen aus dem zeitlichen Verlauf der Verdampfertemperatur 
_gemass Figur 14 auf den Erwarmungsgrad umgerechnet, so ergeben sich die 
Kuryen a und } von Figur 16. Was dabei sofort in die Augen fallt, ist der Um- 
stand, dass der Erwirmungsgrad bei tiefsten Temperaturen 0 betragt. Schuld 
 daran diirfte die voriibergehende Weiterproduktion von Kalte durch Verdamp- 
fen des Ammoniaks sein, wenn der Kompressor abgeschaltet wird. Erst von 
einer Verdampfertemperatur, die héher ist als — 30°C, stimmt das Bild, dass 
_der Aufwarmungsgrad mit steigender Temperatur abnehmen soll, mit unserer 
Erwartung tiberein. 
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Vergleichsweise zeigt Figur 16 mit den Kurven ¢ und d auch den Erwar- 
mungsgrad ohne Heizeinsatz. Fiir diesen sind also nur die Kalteverluste sowie 
die durch den Ventilator erzeugte Warme verantwortlich. 

Hierbei muss noch erwahnt werden, dass die Differenz der Erwarmungs- 
grade mit Heizung nicht der Differenz ohne Heizung zu entsprechen braucht. 
Die Verzogerung des Kanals in der Temperaturanpassung ist ihr nicht ahnlich. 


Figur 16 


Erwarmungsgrad in Funktion der Verdampfertemperatur: a Heizeinsatz 18 kW, Luftgeschwin- 

digkeit vyy = 25 m/s, ausgehend von — 47°C; b Heizeinsatz 18 kW, Luftgeschwindigkeit vyg = 10 m/s, 

ausgehend von —50°C; c ohne Heizung, Luftgeschwindigkeit vjg = 25 m/s; d ohne Heizung, 
Luftgeschwindigkeit vzz = 10 m/s. 


6.3 Die Anderung der Mefstellentemperatur bei gleichbleibender 
Verdampfertemperatur 


Wiahrend der vorhergehende Abschnitt die Kanaltragheit und ihren Einfluss’ 
auf die experimentbestimmende Verdampfertemperatur beschrieb, sollen an 
dieser Stelle die Variationsméglichkeiten der MeBstellentemperatur behandelt 
werden. 7 

Wie wir aus den Beobachtungen an natiirlichen Hagelkérnern wissen [9], 
sind diese vielfach schalenartig aufgebaut mit einem Wechsel klarer und opaker 
Schichten. Sobald nun solche Typen im Hagelversuchskanal reproduziert wer- 
den sollen, sind periodisch andernde Wachstumsbedingungen vonnéten. Die 
Planung sah vor, dass die notwendigen Temperatur- oder auch Feuchteschwan- 
kungen bei gleichbleibender Verdampferblocktemperatur nur durch Variation 
des Heizeinsatzes erzeugt werden sollten. Zu diesem Zwecke kann der Kalte- 
kompressor mit Vorteil durch die in Abschnitt 4.4 beschriebene Steuerung be- 
trieben werden, die die Kalteleistung so reguliert, dass der Verdampfer stets die 
gewiinschte konstante Temperatur aufweist. 

Die durch die Heizung erreichbare grisste Temperaturerhéhung entspre- 
chend 18 kW Leistung ist in Abhangigkeit der Luftgeschwindigkeit als Kurve a 
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in Figur 17 eingetragen. Diese theoretisch mégliche Differenz stellt sich jedoch 
nicht augenblicklich ein. Die tragheitsbedingte Anpassung mit der Zeit ergibt 
sich, wie die durchgezogenen Kurven derselben Figur fiir Zeitintervalle von 
1, 2 und 10 min zeigen, (Die letzte Zeitangabe stellt die Zeitspanne dar, in der 
mit einer praktisch vollstandigen Temperaturanpassung gerechnet werden 
kann, sofern das experimentelle Vorgehen mit dem hier beschriebenen iiber- 


arn 
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30 t—. 


20 


nach 10min | 
~~ —nach 10min 
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0 | "1 
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Figur 17 
Zeitabnangige, durch Heizeinsatz von 18 kW erzeugte ErhGhung der MeBstellentemperatur ¢jy in 
Funktion der Luftgeschwindigkeit vjyg (bei jeweils konstant gehaltener Verdampfertemperatur). 
Kurve a, theoretisch berechnete Werte; Kurve b, Werte, die erreicht werden mit im Moment der 
Aufheizung abgeschaltetem Kaltekompressor. 


_einstimmt.) — Wenn der theoretisch mégliche Wert nach 10 min iiberschritten 


wird, so riibrt das von der Berechnung der umgewalzten Luftmenge auf Grund 
der zentralen MeSstellengeschwindigkeit, ohne Beriicksichtigung der am Rande 
auftretenden reduzierten Werte, her. Kurve } von Figur 17 zeigt die Tempera- 


_turen an, die sich bei jeweils abgeschaltetem Kompressor einstellen wiirden. 


Die Abhangigkeit all dieser Werte von der Verdampfertemperatur ist, wie 
verschiedene Messungen im Bereiche von — 5 bis — 40°C ergeben haben, nicht 


_ von allzu grosser Bedeutung, so dass in Figur 17 Mittelwerte eingesetzt wurden. 


Es besteht jedoch die erwartete Tendenz, dass die zeitliche Verzégerung bei 
tieferen Temperaturen zu kleineren Werten strebt. 

Da gemiss den vorgenommenen Messungen die Riickkithlung auf den Ur- 
sprungswert die gleiche Zeit braucht wie die Aufheizung, ist das Ausmass 


_ erzwungener periodischer Temperaturschwankungen beziiglich ihrer Amplitude 


im Zusammenhang mit der Periodendauer und der Luftgeschwindigkeit be- 


_kannt. Es muss jedoch die Einschrankung gemacht werden, dass die Kom- 


-pressorleistung den jeweiligen Verhaltnissen angepasst werden kann. Die Ver- 


haltnisse komplizieren sich weiter, wenn auch Feuchtigkeit in den Kanal ein- 
gespritzt wird. Diese beeinflusst jedoch, sofern sie konstant bleibt, die Tempe- 
raturvariationsmoglichkeiten gemass Figur 17 praktisch nicht, es sei denn, dass 
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bei zeitweisen Temperaturen tiber 0°C eine zusatzliche Tragheit in Rechnung 
gezogen werden muss, die ihre Begriindung in allfalligen Eisansatzen an der 
Kanalwandung hat. Da aber in Abhangigkeit des Feuchteaufwandes der Ar- 
beitspunkt der Kalteanlage angepasst werden muss, ergeben sich Verhiltnisse, 
die in Figur 18 abgeklart werden sollen. Sie zeigt die grésste mdgliche perio- 


+Aty 
4min 
°6/min 

20 


Y= 20m/S 


0 + ee Ge hee ee °C 
Figur 18 


Darstellung der gréssten, alle vier Minuten erzeugbaren periodischen Erhéhung der Mefstellen- 
temperatur in Funktion der MeBstellentemperatur, bei zwei Luftgeschwindigkeiten und verschie- 
denen durch Feuchte zugetiihrten Energiemengen Qf (Heizeinsatz maximal 18 kW). 


dische Temperaturerhohung bei einer 4 min dauernden Periode in Funktion 
der tiefsten MeBstellentemperatur. Als Parameter tritt die durch Feuchte ein- 
gespritzte spezifische Energie in Megakalorien auf, die MeBstellengeschwindig- 
keit wurde im einen Fall mit 10 m/s, im andern mit 25 m/s angenommen. 
Figur 18 demonstriert, dass bei gleicher Luftgeschwindigkeit oberhalb einer 
Mefstellentemperatur von — 20°C das Ausmass der periodischen Temperatur- 
erhéhung einzig durch den Einsatz der Heizung (jeweils 18 kW) bestimmt wird, 
wahrend unterhalb von — 30°C die Kompressorleistung und die Feuchte die 
Amplitude beschranken. 

Weiter kann festgestellt werden, dass sich bei gleicher Heizung die Ampli- 
tude nicht umgekehrt proportional zur Luftgeschwindigkeit verdindert. Die 
Tragheit des Kanals wirkt sich bei grésseren Temperaturveranderungen ver- 
haltnismAssig starker aus; die Zeit zur Anpassung an den theoretischen Grenz- 
wert steigt. 

Aus all diesen Daten sollte nun die Leistungsfahigkeit des Hagelversuchs- 
kanals hinreichend bekannt sein, um die Méglichkeiten, die dem Experiment 
zur Verfiigung stehen, erfassen zu kénnen. Ob allerdings die erzeugbaren perio- 
dischen Schwankungen der Nachbildung der Natur geniigen, lasst sich disku- 


: 
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tieren. Die Entscheidung auch dieser Frage wird aber erst ein Vergleich ent- 
sprechender natiirlicher mit kiinstlich erzeugten Hagelkérnern ergeben. 


7. Messmethoden 
7.1 Grundsatzliches 


Die Messungen der Experimentparameter sind mit Schwierigkeiten ver- 
bunden, sobald wir den Charakter einer Atmosphare mit unterkithlten Wasser- 
tropfen ermitteln miissen. Diese Tropfen gefrieren bei den vorkommenden Be- 
dingungen nicht ohne Anwesenheit von Eisbildungskernen. Das gilt jedoch nur 
fiir den Fall, wo sie ganz von Luft umgeben sind. Sobald sie auf Oberflachen 
auftreffen, setzt meist sofort der Gefriervorgang ein. — Das heisst, dass jegliche 
Sonden, die wir als Messfiihler in den Windkanal einfiihren, sich sogleich mit 
einem Eisansatz bedecken. Dieser Vorgang kénnte wohl bei Temperaturen, die 
nur wenig unter 0°C liegen, etwas eingedammt werden, bei extremeren Tempe- 
raturen ist jedoch kein Gegenmittel bekannt. Wie sich diese Tatsache auf die 
einzelnen Parameter auswirkt, soll in den folgenden Abschnitten dargestellt 
werden. 


7.2 Die Temperaturmessung 


Der schlimmste Einfluss des Eisansatzes auf Sonden wirkt sich auf die 
Temperaturmessung aus. Infolge Freiwerdens der Gefrierwarme wird der Fihler 
eine gegentiber der Umgebungstemperatur meist zu hohe Temperatur anzeigen. 
— Es muss hier noch vermerkt werden, dass bei diesen Bemerkungen still- 
schweigend angenommen wird, dass die Luft zumindest mit Wasserdampf ge- 
sdttigt, eventuell sogar iibersattigt ist, was auf Grund der Tropfenerzeugungs- 
methodik der Fall sein sollte. 

Als beste Abschirmmethode ist der Schutz eines feinen Thermoelementes 
durch den Windschatten eines diinnwandigen Plexiglasrohres zu erwahnen. 
Aus der Menge des angelagerten, gefrierenden Wassers sowie der Grenzschicht- 


-dicke und der Luftgeschwindigkeit kann die Erwarmung der das Thermoele- 


ment umstromenden Luft abgeschatzt und ein entsprechender Korrekturwert 
angebracht werden. 
Diese Methode erlaubt es aber prinzipiell nur, zu einzelnen Zeitpunkten 


_ Messungen vorzunehmen, da sonst die Verhaltnisse an der MeBstelle allzu stark 
- beeinflusst wiirden. Sie ist also speziell brauchbar bei Experimenten mit statio- 


, 


a ek 


+ 


naren Betriebsbedingungen. 

Als weiterer Weg zur Temperaturmessung wird eine tiber den Kanalquer- 
schnitt integrierende, sondenlose Methode verwendet, die auf der Basis der 
temperaturabhangigen Schallgeschwindigkeit arbeitet. Ein Ton wird durch 
einen Lautsprecher quer durch die MeBstelle geschickt und itber ein Mikrophon 
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empfangen (siehe Figur 5). Die Phasendifferenz der Schwingung zwischen 
Lautsprecher und Mikrophon ist direkt ein Mass fur die Laufzeit des Tones 
und kann in Temperatur umgerechnet werden. Da der Einsatz dieser in Erpro- 
bung stehenden Anlage noch nicht ganz befriedigt, sei nicht naher darauf ein- 
getreten. 

7.3 Die Feuchtemessung 


Als genaueste Methode zur Ermittlung des totalen Wassergehaltes wird ein 
Absaugen von Kanalluft unter praktisch augenblicklichem Verdampfen der 
Tropfen betrachtet. Die dazu verwendete spezielle Sonde (siehe Figur 5) ist 
heizbar ausgefiihrt und erlaubt es, durch Einbau von DruckmeBstellen an deren 
Innen- und Aussenwand eine zur Kanalgeschwindigkeit passende Absauge- 
geschwindigkeit zu wahlen. Abgesaugt wird mittels emer Vakuumpumpe. Vor 
derselben befindet sich ein Absorptionsgefaéss mit Phosphorpentoxyd. Dieses 
nimmt allen Wasserdampf auf, so dass sich bei Kenntnis der durchgesetzten 
Luftmenge der totale Wassergehalt pro Kubikmeter wagen und errechnen lasst. 

Bei Kenntnis der Lufttemperatur ist es weiter méglich, unter Annahme 
einer Wasserdampfsattigung der Kanalluft auf den freien Wassergehalt zu 
schliessen. 

7.4 Die Messung der Luftgeschwindighest 


Die Geschwindigkeitsmessung, insbesondere ihre Eichung, erfolgt normaler- 
weise mit Prandtl-Rohren. Diese sind aber bei vereisenden Bedingungen nur zu 
Punktmessungen zu verwenden, wahrenddem zur Ermittlung laufender Daten 
der Druckabfall an den Enden des vor der MeBstelle sich befindenden Effusors 
beigezogen wird. Die Vereisung dieser Stellen ist normalerweise stark reduziert. 
Es muss aber darauf verwiesen werden, dass infolge von Querschnittsverande- 
rungen durch Eisansatz hin und wieder Eichpunkte mit dem Prandtl-Rohr 
eingemessen werden miissen. Als MeBsysteme werden ein Prazisions-Wasser- 
mikromanometer sowie zwei ebenso genaue Druckwandler mit elektrischer An- 
zelge verwendet. 


7.5 Die Eisbildungskernzahlung 


Die Aufnahme der Eisbildungskerncharakteristik der Kanalluft wird eher 
nur summarisch vorgenommen. Es ist vorlaufig vorgesehen, dass als Fremd- 
kerne nur Silberjodidpartikeln verwendet werden sollen, deren Konzentration 
Null, wenig (100/1 Luft) oder viel (105/l) betragen soll. Die Messung geschieht 
nach der bereits frither beschriebenen Methode [13], wobei noch erwahnt werden 
konnte, dass die Proben der zu untersuchenden Kanalluft ebenfalls iiber die 
Feuchtesonde abgezapft werden. 

Wenn auch verschiedene, hier in Abschnitt 7 beschriebene Messmethodent 
noch nicht befriedigen, so ist doch zu erwarten, dass im Laufe der Experimente 
eine weitere Entwicklung zu leistungsfahigeren Apparaten stattfindet. 


; 
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Mit dem Hagelversuchskanal wurden bereits Vereisungsversuche vorgenom- 
men. Die ersten Experimente zeigten dabei, dass die Anlage den gestellten 
Anforderungen im grossen und ganzen-geniigt und dass durch deren Einsatz 
tatsachlich erwartet werden darf, die Kenntnis iiber die Hagelbildung merklich 
erweitern zu kénnen. 


8. Nachwort 


Der Hagelversuchskanal wurde im Rahmen des Versuchsprogrammes der 
Eidgendssischen Kommission zum Studium der Hagelbildung und der Hagel- 
abwehr erbaut. Ihrem Prasidenten, Herrn Professor Dr. R. SANGER, bin ich ver- 
ptlichtet fiir seinen Einsatz zur Verwirklichung des Projektes; Herrn Professor 
Dr. J. ACKERET méchte ich fiir seine stete Hilfsbereitschaft und die wertvollen 
Ratschlage herzlich danken. 

Erwahnt sei an dieser Stelle auch die Mithilfe von Herrn Dr. H. SPRENGER 
vom Aerodynamischen Institut der ETH sowie von Herrn Dr. F. FURRER 
vom Institut fiir Hochfrequenztechnik der ETH. 

Zum Gelingen des Projektes haben durch ihren grossen Einsatz auch die 
Firmen Gebriider Sulzer AG, Winterthur (ihr verdanken wir auch die Unter- 
lagen zu den Figuren 1, 2 und 4); Camille Bauer AG, Basel; und Sprecher & 
Schuh AG, Aarau, beigetragen. Ihnen wie auch den andern Firmen, die kleinere 
Auftrage ausfiihrten — es seien dabei erwahnt Trion AG, Ziirich; CIBA, Basel; 
Chemische Werke Neuallschwil; E. und H. Huber AG, Pfaffikon; J. Huber AG, 
Wiirenlos; Sauter AG, Basel; H. Siegrist & Co., Winterthur; Elcalor AG, 
Aarau —, gebiihrt fiir die tadellose Arbeit Anerkennung. 

Zu grésstem Dank bin ich jedoch meinem Chef, Herrn Dr. M. DE QUERVAIN, 
verpflichtet, der die beschriebenen Arbeiten leitete und in allen Phasen der 
Planung und des Aufbaues wesentliches zur Abklarung der einzelnen Probleme 
beigetragen hat. Gedenken méchte ich auch meiner beiden Mitarbeiter, Herrn 
F. Peter und Herrn P. RAss, die tatkraftig bei den experimentellen Vorarbei- ° 
ten und der Montage des Kanals mitwirkten. 

Der Bau des Hagelversuchskanals wurde erméglicht, indem der Schweize- 
rische Nationalfonds die bendtigten bedeutenden Kredite bereitstellte; das 
Laboratorium wurde durch die Abteilung fiir Landwirtschaft des Eidgendssi- 
schen Volkswirtschaftsdepartements finanziert. Der Eidgenéssischen Inspektion 
fiir Forstwesen, Jagd und Fischerei sei bei dieser Gelegenheit dafiir gedankt, 
‘dass sie die Aufnahme der Hagelforschung in das Arbeitsgebiet des Schnee- 
forschungsinstitutes erméglicht hat. 
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Summary 


In the past the problem of fundamental research into hail formations has been 
dealt with in a more or less theoretical fashion on the basis of meteorological 
information. The attempt will now be made to reach an interpretation of how 
hail forms by an experimental physical method. To this end a wind tunnel has 
been built in which experiments can be conducted for studying the growth o: 
suspended or floating particles of ice in terms of environmental conditions. Sub. 
sequent comparison with natural hailstones suggests that conclusions can be drawr 
as to how these have come into being and that eventually an explanation o: 
hail formation can be given which comes closer to reality. 

For the reproduction of atmospheric conditions in the laboratory the indi 
vidual parameters of the experiment have been specified as follows: Temperature 
down to —40°C; maximum free water content in the air 10 g/m3; size of drops 
between 1 and 200; highest wind speed 25 m/s; concentration of ice-forming 
nuclei of various substances: none, single, multiple. Although variations it 
pressure have not been employed, the laboratory conditions may be transformeé¢ 
to those of the open atmosphere by using the rules of similarity. i 

The construction of the hail tunnel (Figure 1) may be described as follows 
The blowing apparatus / drives the air into a temperature adjustor and filte: 
unit 2; details with an air-cooler 2, 19 kW heater 3 and electrostatic filter - 
are shown in Figure 2. After these units the air is led into a vertical sectiox 
where it is made humid and homogeneous. The ice-forming nuclei are also addec 
at this point. Immediately next to this comes the actual measuring section 4 
where the growth of a suspended or floating object is to be studied. The remaining 
parts of the tunnel 5 serve only to return the circulated air to the ventilator 
Figure 3 shows the refrigerating compressor, Figure 6 the control unit. The arrange 
ment of the drop-producing installations, such as the steam boiler 7, air-wate 
compressor 4, and rotor atomizer 3, can be understood from the plan of the 
aad _ bes ee ee the construction of the measuring section 

‘ e€ apparatus of the control unit is set up. For 


u 
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more detailed account of the way individual parts (particularly the control 
mechanism) of the equipment work, the reader iS referred to section 4. 

The performance of the hail tunnel has been measured as follows: Assuming 
in accordance with the values given in Figures 8 and 9 that injected humidity and 
added heat are interchangeable, various stationary operating conditions in terms 
of heat production were allowed to develop with the refrigerating compressor 
running at maximum performance. The maximum performance curves measured 
in this way can be derived from Figure 10 for the compressor revolution numbers I 
and II and different air speeds. These curves represent the effective energy avail- 
able for the experiment as a function of the air-cooler temperature ¢y, since all 
losses from the tunnel are involved in the readings and therefore included in the 
final result. On the basis of formula (8) and Figures 11 and 12 the fundamental 
temperature for the output of the refrigeration plant ¢, can be correlated with 
the temperature in the measuring section ¢,. The curve Lx, gives in every case 
the maximum work of the compressor as a function of ¢y. 

If we convert corresponding energy values and air-cooler temperatures into 
humidity and measuring-section temperature, the correlations which result are 

shown in Figure 13. The specification of a free water content of 10 g/m® is in fact 
exceeded over a wide range of temperatures. 

Since variations of the different factors play an important part in the imita- 
tion of hail-forming atmospheric conditions, the behaviour of the hail tunnel was 
also recorded during unstable operation. Figures 14, 15 and 16 show corresponding 
sequences in the governing temperature of the air cooler. Even more interesting, 
however, is the influence of the heat expenditure while the temperature of the 
vaporizer is kept constant. With a heat input of 18 kW there resulted, depending 

on the air speed, increases of temperature in the measuring section, as recorded 
for different time intervals in Figure 17. The observed time dependence is due to 
the inertia of the tunnel. 

As cooling-down takes as long as warming-up, we can calculate even forced 
‘periodic temperature fluctuations. Figure 18 gives an example of this for a period 
lasting 4 min and for air speeds of 10 and 25 m/s. At higher temperatures the 
largest fluctuation in temperature which can be produced is dependent only on 
the air speed; at lower temperatures it is limited by the capacity of the compres- 

-sor. If, however, humidity is introduced under these latter conditions — and this 
will be the case in the main experiments — then the total energy added increases; | 
adjustment of the cold compressor can only take place at higher temperatures, so 
that the limiting stability curves are correspondingly displaced. The extent of 
this displacement for various amounts of additional energy due to humidity QO, is 

shown in Figure 18. 
' The plant was constructed as part of the research programme of the Swiss 
Federal Commission for the Study of Hail Formation and Prevention. The results 
already obtained make it possible to expect that the tunnel will lead to a sub- 
stantial increase in our knowledge of hail formation. 
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Bemerkung zur Berechnung spezieller Potenzsummen 


Von Rupo_F HorMann, Berlin, Deutschland?) 


Die einschlagigen neueren?) Formeltafeln und Tabellen enthalten unter ande- 
rem, soweit dem Verfasser bekannt, nur die Potenzsummen der natiirlichen Zah- 
len in alilgemeiner Form. Die Summen der Potenzen der trigonometrischen Funk- 
tionen werden jedoch nur fiir ganzzahlige Exponenten M = 1, 2, 3, 4 angegeben. 
Es erscheint daher:der Hinweis angebracht, dass sich die entsprechenden allge- 
meinen Summenformeln explizit angeben lassen. Dabei wird die Methode zu 
ihrer Berechnung nur gestreift. 

Einen Beitrag zur Aufsummierung unendlicher Reihen in geschlossener Form 
mit Hilfe der Laplace-Transformation hat WHEELON®) geliefert. Dabei wird 
davon ausgegangen, dass die Summanden der vorgegebenen Reihe Bildfunktio- 
nen von bekannten Originalfunktionen sind. 

Es liegt nun nahe, auch den umgekehrten Weg zur Berechnung von Reihen 
einzuschlagen, und zwar dann, wenn die Summanden der Reihe Originalfunk- 
tionen bekannter Funktionen sind, die durch die komplexe Umkehrformel der 
Laplace-Transformation erhalten werden. é 

Sei namlich 


S= J’ A(h) (yy 
k 
gegeben*), ferner A(#) in der Form B(hk) - C(k) und 
t ¢+400 : 
Sey ee I 
Bk) = 5 “ NOW ek (2) 
c¢—100 


durch die komplexe Umkehrformel der Laplace-Transformation darstellbar, dann 
geht (1) iiber in 


c+100 


ies Dy C(k) on i b(s) es ds . (3) 


c—ico 


Das Integral ist dabei iiber eine Vertikale mit der Abszisse c in der Konvergenz- 
halbebene von O(s) zu erstrecken derart, dass alle Singularitaten von d(s) links 
der Vertikalen liegen. In allen Fallen, in denen es mdglich ist, die Reihenfolge von 
Summe und Integral zu vertauschen, was zum Beispiel bei endlichen Reihen 


.% 
] 


1) Osram Studiengesellschaft. 


2) I. M. Rysuix und I. S. GrapsteIn, Summen-, P i } 
PiepelG: ; -, Produkt- und Integraltafeln (Berlin 1957). 

8) A.D. WueEton, J. appl. Phys. 25, 118 (1954). : a ae e 
n 

iz 


rie Eee : . 
) Die Summanden seien als Funktionen der nur fiir ganzzahlige Werte vorkommende 


Grosse k aufgefasst. H 
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stets der Fall ist, geht dann (3) iiber in 


; € +400 
S = 50s) o(s) ds (4) 
c—t0o 


y(s) = os C(k) ek, 
k 


mit 


Damit ist fiir die Reihe S — unter den oben angegebenen Voraussetzungen — eine 
geschlossene Darstellung gefunden. 

Im folgenden soll die Methode anhand einiger Beispiele demonstriert werden; 
es werden mit ihrer Hilfe die Summen der Potenzen der trigonometrischen 


Funktionen 
nm nn 
y sin@k xy und y cosVk x 
k=1 k=1 


sowie die alternierende Potenzsumme der natiirlichen Zahlen 


n 
“x: (—1)* kM 
k=1 
fiir beliebige ganze Potenzen M > 0 hergeleitet. 


1. Mittels der Korrespondenz®) 


2 ZN 
sin“” « to—e - a Maes UND==e 10 P2, 2.) 


s [fis*+ (2.24) 7] 


y : 
_folgt fiir die Reihe ¥’ sin?’ & x gemiss (4) 


k= 
x €+400 ihe i 
3 n : a at 
¢ > sin?’ k x = ee ieee 1 { ds . (5) 
% mt 
net e—ico ra) s [J+ (2 #)7) 


: (5) lasst sich leicht auswerten. Die Partialbruchzerlegung 


x2 {/2N\ 1 / ters! 
oe SSA Jr aay Sore Hes i fs? + (2k x)? 
s [] s+ + (2k x)*] 


und die Anwendung des Residuensatzes liefert 


- 
4 
> 
- 
“ 
, 
| 
f 


S 1 {/2N : AN eo oe [ | 
Sent a x= fe |?) oi + Sar (y wea + v4 2h xi) (6) 


_k=1 


SN 


2 5), Die Bezeichnungsweise ist G. Dortscn, Tabellen zur Laplace-Transjormation (Berlin 1947), 
entnommen. 
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mit (ens — 1) 
p(s) = “(e* — 1) ae 
Es ist 
sinh (2 ” + 1) > 
p(s) + y(-s)==1+ —, 


Ss 
He cial 
sin 2 


: s 
sinh (v7 + 1) > sinhn — 


2 2 
WAS) seul =.S) eae eee ee 
sinh — 
Damit folgt aus (6) 
N Dae 
a. 2N i be 2 N\ <aye( 2N \ sink (20+ 1) % 
pat: ote g2N+1 ie 1)( N )+ ae - N — 1) sink x Se 
und 
N 
2N 2N 
Se aes 
a Dien (va) = (Cw) 8) 


Auf dieselbe Weise bekommt man wegen 


N ' , 
: 2 
sia’)! oboe (—1)* ( nb) eben Nis 0,0 en 


AM Nice ues ae Cire Wea . 
n N 
2N+1 1 2N + 1 
sin?V+1 , y = —_ Ena ) 
> ne 2 N—k . 
sin(2 k + 1) (n+ 1) = gin (2 6 Bajo OM 
x = 2 —- j 
sin (2 & + 1). ; 
wegen 
ON |. fi F ‘i 1 N 2N 2 : 
cos”"> *t o—e —— Beale) ar. 
‘| Mena 2 (eae [sta (2ka)] fo) oe 


(10) 


N—k sink x 


DO ee ee a kas | 


(2+ 1) eels br) a) eae 


mit 


1 ( 2N 2N 
Pe ee. (11 
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und wegen 


N 
2N+1 1 Caan s 
cos xt o—e —— cr SS ee Nees J 
aw Ne Eee SN Oe Baths Fae 


n a Nae sin(2k + 1) (2+ 1) wig 
ESS te be es i 2N+1 ‘baet T. kon (12) 
k= 2 22N +1 — Mes k é xv 
2 » hed sin(2k + 1) > 
mit 


= men 1} 


k=1 
folgt wegen 
! 
tN o—e mks ? N= 0, L 2 
und 
re ea ews 
spl e+ 
“aus (4) die Beziehung 
n on 
N! [(—1)” es —1] es 
k Bl 
» (-1) RN = sai? et eC L (14) 
k=l x (e8 + 1) s 


da wegen der fiir N gemachten Voraussetzung der Integrationsweg auf einen 
Kreis K um den Nullpunkt verlegt werden kann. Zur Auswertung von (14) wird 


[(-1y" ens — 1] 


ae Eo = f(s) 


entwickelt. Es sei 


“ 


z IIe 


“wobei die A, der Rekursionsformel 


4 p 
E ne ‘ee 2 
4 ee Sy ee (7) 4;= 9. p=1,2,..., 
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geniigen oder aus der Determinantendarstellung 


| 1 2 0 0 
1 Z Z 0 0 
ee es 3 3 2 0 
(2199 Ul oe vets. os poe DURE ine a ae 2 
A, Qp+1 1 (15) 


Ce) Cai NG 
OS a eae 


sehr einfach zu ermitteln sind. Die A, diirfen somit als bekannt vorausgesetzt 
werden. Mithin wird 


as! ( yy (hk A 1)” Lyre 
w= Ser ( XY (3) 4am + yt a), 


k=0 


so dass sich schliesslich aus (14) 


n N 
en) EN eA nits eal oe Ny ed A, (a 21h Nee (16) 
k=1 k=0 
ergibt. 
Summary 


Some series can be computed easily if the terms are original functions of known. 
Laplace transforms. The method is briefly outlined, as example the series : 


n 


n n 
a sinMk x , i cosMkw and Ds (—1)*§ kM 
k=1 kal 


k=l 
are computed in closed forms for integer exponents M> 0. 
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Minimum Weight Design of Circular Plates Under Arbitrary Loading!) 


By WILLIAM PRAGER®) and RicHarD T. SHiELp$), Providence, R.I., USA. 


To Professor F. K. G. Opevist on the occasion of his 60th anniversary 


1. Introduction 


The direct design procedures developed by Drucker and SuIELp [1, 2]4) have 
been used recently [3] to obtain the minimum weight design of circular plates 
under rotationally symmetric pressure distributions. The sandwich plate and the 
solid or homogeneous plate were considered and designs for both simply sup- 
ported and built-in edge conditions were derived. As an extension of this work, 
the minimum weight design is obtained in the following for a circular sandwich 
plate loaded by an arbitrary distribution of pressure. The deformation modes used 
to obtain the designs are the same as those used in [3] for rotationally symmetric 
pressure distributions. 


2. Definition of the Problem 


The ideal sandwich plate is composed of a core of constant thickness H 
between two identical face sheets of variable thickness h, where h < H. The core 
carries shear stress only while the face sheets carry direct stresses and so provide 
the bending moment across a section. The material of the face sheets is assumed 
to be elastic-perfectly plastic and to obey TreEsca’s yield condition, with yield 
stress o) in tension or compression. It is convenient to represent the principal 
bending moments M,, M, at a generic point of the plate by a point with rectan- 
gular coordinates (M,, M,) (see Figure 1). With Tresca’s yield criterion, the 
stress point must then be within or on the hexagon in the figure [4], the maximum ~ 
bending moment /, being given by 


M,=oHh. (1) 


For a point on the hexagon, plastic bending of the plate can occur. The principal 
curvature rates x,, x, associated with a mode of plastic deformation can be 
represented in Figure 1 by a vector with components proportional to .%, x2 
attached to the stress point in question. The flow rule requires the curvature rate 
vector to be normal to the hexagon for points on a side, and at the corners of the 
hexagon the vector must lie in the fan bounded by normals to adjacent sides. The 
rate of dissipation of energy D4 per unit area of the middle surface due to plastic 


1) The results presented in this paper were obtained in the course of research sponsored by the 


: Office of Ordnance Research, Department of the Army, under Contract No, DA-19-020-ORD-4564. 


2) Professor of Applied Mechanics, Brown University, Providence 12, Rhode Island. 

3) Associate Professor of Applied Mathematics, Brown University, Providence 12, Rhode 
Island. 

4) Numbers in brackets refer to References, page 425. 
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ion is given by 
action 1s gl ; Dz =M, 4+ M, x (2) 


is uniquely determined by the curvature rates 1, %9. 
see ee to determine ‘hie thickness h of the face sheets so that the plate 
can just carry a given distribution of pressure and so that minimum weight of 
material is involved. The material is assumed to be homogeneous and minimum 
weight coincides with minimum volume. The minimum weight design is achieved 


Figure 1 
Yield condition. 


by a plate designed to collapse in a mode for which the condition 


a = const (3)9 
is satisfied over the plate [1, 2], if body forces are neglected. Condition (3) is a 
condition on the rate of deflection w of the middle surface of the plate and the 
thickness / does not enter into condition (3) because of the linear dependence of — 
D4 on h. To this extent, the condition is independent of the pressure distribution 
over the plate but the form of the condition depends upon the position of the 
stress point on the hexagon, which in turn is influenced by the loading. 

For a circular plate with conditions of simple support or built-in support at 
the edge of the plate, the deflection mode w satisfying (3) is rotationally sym- 
metric. With polar coordinates (ry, 6), the edge of the plate being y = R, the 
curvature rate x, is zero and isotropy then requires the bending moment M9 
to be zero. For equilibrium, the bending moments M,,, Mg satisfy the equation 


0? 0M 4 1 027M . 
Or (, M,) Or a ye rts =p (4) 


where p= P(r, 9) is the pressure over the plate. For definiteness, the plate is 
taken to be horizontal. Positive values of M,, Mg stress the lower surface of the 


plate in tension and the pressure is applied to the upper surface. The shear forces — 
QO,, Og are given by 


0M, 1 1 0M 
Qr= ap + > (Mp— Me), Qo=—-— 8. (5) 
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3. Simply Supported Plate 


When the edge y = FR of the plate is simply supported, the minimum weight 
design is obtained by assuming that the whole of the plate is stressed to state 4 
of Figure 1. The deflection rate w (measured in the downwards direction), asso- 
ciated with point A and satisfying (3), and the condition w= 0 at r= R, is 
given by [3] 

w= «(R? —r?), (6) 


where « is a positive constant. The bending moments M,, Mg are both equal to 
the maximum moment M, and substitution in the equation of equilibrium gives 


aM, . 1 0M, 1 02M, 
ot ty ar + a oge  P- @) 


This equation and the condition M, = 0 at y= R are sufficient to determine M, 
over the plate. 

For a concentrated load P at the point (g, y), the solution of (7) satisfying 
wa — Oat ¥ = 1s 


Te, Y Q 
ih i — — Poy oe 
Oo t - ee in| aks (8) 
where 7,, ”, are the distances of the point (vy, 6) from the point (@, m) and its 
inverse with respect to the circle y = R, that is the point (R?/0, m), respectively: 


R4 2 
ri = 0? + r? — 2 orcos(6@— g), Wm De tr? — 2 1 c0s(0 = 9). (9) 


The solution for a distributed load p(7, 6) can be obtained from the funda- 
mental integral (8) by integration: 


22 
> 


{ R? + 7? 0?/R* — 2 ax cos(6 — g) 


R 
1 2 
oo Hh= | P(e, y) In) o? + 72 — 2 or cos(6 — 9) 
o 0 


| ¢ de dp. (10) 


4. Built-in Plate 
For the plate with a built-in edge, stress point C of Figure 1 applies near the 


edge of the plate and stress point A applies in the center of the plate. The velocity. 
field 


waa(lat—r) (0srsa), wa2a(R-y? (@SrSR), (1) 


where a = 2/3 R, satisfies the condition (3) over the plate, together with the 
conditions w = 0w/dvy = 0 at the clamped edge y= R and w and dw/dr con- 

tinuous at y = a [3]. The deformation mode (11) is associated with stress point 4 
for 0<yr<a and stress point C for axv< Rk. The moments M,, Mg are there- 


fore given by 
M,=M,=M, @s71<4);. M,=—My, Me =o (BES Ti) 6 (12) 
For equilibrium, it follows from equations (4) and (12) that M, must satisfy 
equation (7) for 0<y<a, and fora<r=sR 


= es ((M,)=rp. (13) 
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At the junction 7 = a, M, is continuous so that M, = 0 at ry = a. In addition, 
the shear force Q, must be continuous at 7 = a (in the absence of concentrated 
load) and this requires | 0M,/0r | to be continuous. Equation (7) and the condition 
M, = 0 atv = aare sufficient to determine M, for 0<7< a. Equation (13) and 
the condition M, = 0 at y = a together with the known derivative (0M,/07),-4 
are then sufficient to determine M, fora<7< R. 

We first consider a concentrated load P at a point (g, y) within the circle 
y = a, so that 9 < a. As in the previous section, the solution of equation (7) with 
p = 0 and having the appropriate singularity at (@, gy) and M,= 0 atr=ais 


eH eae a in| 7 2] (<a), (14) 
where now 
: ; ; at a2 
r= o+r?—2ercos(6—g), = mi Z 5 vy cos(6— 9). (15) 


For r= a, M, satisfies equation (13) with = 0 and as My= 0 at r= a, we- 
obtain 
M, = 1(6) (1 & <) (axrs Rj. 


Y 


The positive function /(@) entering into this equation is determined by the 
condition that the magnitude of the derivative 0M,/0r is continuous at 7 = a. 
This requires 


Gee iis Otte fa eee i (a<7r<R 16) 

°~ 20 a*%+ @?— 20a cos(0 — 9) 7) 4s 12 Fenn 

The design for a load f(r, 0) distributed inside the circle y = a and zero load — 

over the annulus a<Xv< FR can be obtained from the design (14), (16) by 
integration: 


1 2% @ ’ rehcois 
SN OF aarp BEN eben et Lobb 8 0) 
0 amen rm) al eo? + 72 — 2 er cos(6 — ¢) je. eee (17) 


fory<a,andforaxs7r<R 


2% @ 


1 a 2 2 
qghH— (1-2) ff a? — Q 
a r). A ?) att ght oa costa glee |: ne 


_ When there is no load inside the circle y = a, M, is harmonic for y< a and 
since My, = 0 at y= a, it follows that M, is zero identically inside v =a. We 
assume that the plate is loaded in the annular region by a line load & per unit. 
length along the arc element p< 6< p+ A@ of the circle 7 = eo, where @>a 
(Figure 2). The total load on the plate is then P= & 0 AO. Because p a5 in 
equation (13) except for 7 = 0, pS OS p+ AO, and also M, = 0M,/or = 0 atl 
Y S © oe : ane om equation (13) that M, is zero except in the tee ted sector 
pe edo ra Y Ae ; shown shaded in Figure 2. Thus the load is carried 
ene ae the point on the support circle y = R nearest to the region 

f Pp : ee of the load. In order to determine M o in the shaded region, we 
note that iy 1s zero at y = g and the shear force Q,, and therefore 0M,/or, has a 
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discontinuity of amount ® across the loaded element of arc. The result is 


% Hh=My= ea (1 e 2 (19) 


in the region ga<7v< Rk, pS OS qy+ AO. 


Figure 2 


Line load on arc element 4B. 


‘Integration of the design (19) leads to the design for a distributed load P(r, 6) 
which is zero inside the circle 7 = a: 


r 
a Hh = | ple, 6) e(1~ 8) ae (20) 
a 
_ for a<yv< R and the thickness fA is zero for y < a. The upper limit of integration 
_ in equation (20) is y because the pressure f(@, 6) at the point (g, 0) adds material 
only to points (v, 6) for which 9<y. 
When the pressure distribution (7, 0) is non-zero in both ySaandasx7rsR, 
the thickness distributions given by equations (17), (18) and equation (20) are 
added. 
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Zusammenfassung 


Das Dimensionierungsverfahren von DRUCKER und SHIELD [1, 2] wurde kiirz- 
lich auf die Dimensionierung rotationssymmetrisch belasteter Kreisplatten fiir 
Mindestgewicht angewandt [3], wobei sowohl gelenkig gestiitzte als auch ein- 
gespannte Sandwichplatten betrachtet wurden. In der vorliegenden Arbeit 
werden diese Untersuchungen auf kreisférmige Sandwichplatten unter nicht 
rotationssymmetrischer Belastung ausgedehnt. Es zeigt sich, dass auch fiir solche 
Lasten der Dimensionierung fiir Mindestgewicht diejenigen Verformungszustande 
zugrunde gelegt werden kénnen, welche schon zur Losung der entsprechenden 
Aufgabe bei rotationssymmetrischer Belastung konstruiert wurden. Fiir Bela- 
stung einer gelenkig gestiitzten oder eingespannten Kreisplatte durch eine Ein- 
zellast mit beliebigem Angriffspunkt werden explizite Dimensionierungsformeln 
gewonnen, fiir beliebig verteilte Belastung werden Integraldarstellungen der 
optimalen Dimensionierung gegeben. 
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A Note on Addition Theorems for Mathieu Functions 


By K. S#RMaARK, Copenhagen, Denmark?) 


It is the purpose of this note to point out the existence of an addition theorem 
for Mathieu functions besides the one given by MEIXNER and SCHAFKE [1-—3]?). 
The latter states — in the notation of [1] — the following. We consider two elliptical | 
cylinder coordinate systems A and B: 


A: #, +1%,=e¢cosh(z +71); 


Be a4 25 4 4s = 65 cosh (Zz, 22 4), 


a es 


connected by 
ecosh(z + it) = Je+t** 4 ¢, e~** cosh(z, + 7%,) . 


(We consider only teal, positive values of c, Cy, 2, 29, f and t); see figure.) Let 
z, and z, be the elliptical coordinates — in the A-system — of the focal points 
of the B-system and z,,,, the largest one of z,, 2), then 


} co 
MY (29; No) Me y(ty; Ab) = OE A, M9) (2; h) Me p+ s(t; h?) ( 
= 1) 


_ : s = 
0] = ies Px 3) Gye 4; Di integer = 0) 


will be valid in the region 
Eee ((St2 2 ax). 
In (1) MD (2; 
a (1) Mpi(z; h) me,(t; h?) and Mess ho) me, (ty; hi) are solutions of the wave- 
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equation in elliptical cylinder coordinates: 


Ag+ kh? p= 0 
andh=kc/2, hh=ke,/2. 

The expression for A, may be found in [1]. 

The formulation of this addition theorem given by MEIXNER and ScHAFKE 
is more general than the one given here, but for most physical applications the 
latter will be sufficient. However, the region z< z,,,, is — from the point of view 
of physical application — as useful as the region z > 2,,,,. We shall now show, by 
means of theorems given in [1], that in a part of this region another form of ad- 
dition theorem may be given, connecting the functions M')(z9; hy) me y(t; h3) 


(7 = 1, 2, 3 or 4) with the functions MMe; h) me,(t; h?). 


The mutual arrangement of the two elliptical cylinder coordinate systems A: (3, t) and B: (2p, tq) 
with centers at A’ and B’ respectively. J is the distance between A’ and B’. @ is the angle between r 
and /. (4) is valid inside the heavy line. 


Let (R, 0) and (r, «) be polar coordinate systems with centers at A’ and B’, 
respectively. Then according to [1] 


[ee 
j j i(p+s)a. 
DAT gist) ne gltg; he) = DS" a, ZO. (hr) atte "hs e 
s=— oo < 


day 5 = (—1)' CE, (hg); 4 


21, 0 (j = 1, 2, 3 or 4; p = integer = 0) 


Dy oS ee a 


will be valid for vy >c,. Here ZY) denotes, corresponding to Gia lela, 
respectively, Bessel, Neumann and Hankel functions of the first and second 
kind. Again according to [1] 


co 
Jy R) PPO — Da, Mp) slash) mepyalts H); : 
é s=—oo ( ) 
f +2 a ee =e = 
dens ts (ay Ch+2 (nd); Uae =0 (p= integer = 0) 


will be valid for z= 0. The quantities C{, occuring in (2) and (3) are the coeffi- 
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cients in the expansions of the periodic Mathieu functions: 


ee 2) 
ee a 3 


]=—=00 


By means of a well-known addition theorem for the cylinder functions: 


(21)% Z(k 1) Verte = Seis ZO, (ke l) , 
s=—co 
valid for R </, and (2) and (3) one then gets, observing that 06 = y+ u and 
at@?=atu-o, 
Mz; ho) me,,(to; Ab) = a, BY MIN h) me, (t; h?) (4) 
(j= 12, 3,4; m = integer = 0) 
where 


ie. 2) 

= Did cp OYE EY ge 

S, p=->00 
(4) will be valid outside a circle with center at B’ and r = Cp, but inside the circle 
with center at A’ and R = /. This region is shown in the figure. 

Comparing the expression for Bo with the expression given in [1] for A, in 

(1), it is interesting to note that the major difference is the occurrence of Z")(k 7) 
instead of /(k 1). This seems to indicate that a more fundamental derivation than ; 
the one given here is possible. 


Note Added in Proof: Actually the region of validity of the addition theorem (1) — 
is somewhat larger than stated by MEIXNER and SCHAFKE. It may be shown that 
(1) holds in the region outside the circle with centre at A’ and radius / and the 
circle with centre at B’ and radius Cp. 


REFERENCES 


(1] J. Merxner and F. W. ScHArKe, Mathieusche Funktionen und Sphéroid- 
funktionen, 1st ed. (Springer Verlag, Berlin 1954), p. 171-187. 

[2] F. W. ScHAFKE, Math. Z. 58, 436-447 (1953). 

[3] J. Merxner, Math. Nachr. 3, 14-19 (1949). 


Zusammenfassung 


Es wird gezeigt, dass ein Additionstheorem fiir Mathieusche Funktionen — 
verschieden von dem von MEIXNER und SCHAFKE gegebenen — existiert. 


(Received: December 1, 1958.) 
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Erganzung zum Aufsatz: 
«Der senkrechte Verdichtungsstoss am gekriimmten Profil»*) 


Von JURGEN ZiEREP, Aachen, Deutschland?) 


In der angegebenen Veréffentlichung wird das Verhalten eines senkrechten 
Verdichtungsstosses an einer gekriimmten Wand untersucht. Das Hauptergebnis 
ist, dass die Stosskriimmung an der Wand unendlich wird. Anschliessend werden 
einige Folgerungen im Hinblick auf einschlagige Versuchsergebnisse gezogen. 
Herr RANpDALL (Royal Aircraft Establishment, Farnborough) hat mich in dan- 
kenswerter Weise darauf hingewiesen, dass die Konsequenzen der mitgeteilten 
Ergebnisse weiter reichen. Wahrend in meiner Note im wesentlichen nur mit dem 
Absolutbetrag der Wandkriimmung gearbeitet wird, zeigt es sich namlich, wenn 
man das Vorzeichen derselben von Anfang an mitberiicksichtigt, dass man zu den 
folgenden Resultaten gelangt. Mit Hilfe der Darstellung (3.3) fiir die Stosskurve 
(jetzt giiltig fiir B= = 0!) 


4 o 
y= — — = 3/2 
z 24/3 Bz H 


erkennt man, dass zum Beispiel im Mach-Zahl-Bereich o(M) > 0 ein Stoss nur 
moglich ist, solange B= > 0 ist. Dies wiirde aber heissen, dass sich ein senkrechter 
Stoss (in diesem Mach-Zahl-Bereich) nur an einem konkaven Profil, dagegen 
nicht an einem konvexen Profil ausbilden kann. Dieses Ergebnis ist tiberraschend. 
Herr GARABEDIAN (Stanford University, Calif.) wies in einer miindlichen Mittei- 
lung jedoch darauf hin, dass bei diesem Schluss ganz wesentlich eine Stetigkeits- 
voraussetzung benutzt wird. Wir haben namlich in (2.6) angenommen, dass beim 
Grenziibergang 7 > 0 h, Bz gegen die Wandkriimmung konvergiert. Inwieweit 
diese Voraussetzung bei dem langs des Stosses durchgefiihrten Limes 7—> 0 
wirklich erfiillt ist und damit der oben beim konvexen Profil gezogene Schluss 
bindend ist, bleibt vorerst noch offen. Dariiber soll hier spater ausfiihrlich be- 
richtet werden. Zum Abschluss wollen wir noch folgende Bemerkung tiber das 
zu Beginn der Arbeit erwahnte orthogonale (&, 7)-Netz machen. Aus dem Text 
ist unmittelbar ersichtlich dass dieses Netz in der Note eigentlich gar nicht be- 
nétigt wird. Was gemacht wird, ist das Folgende: Die in den Bewegungsgleichungen 
-auftretenden Ableitungen 0u/dx, Ou/dy, Ov/0x, dv/dy werden lediglich auf der 
Stosskurve in Ableitungen langs des Stosses und senkrecht dazu umgeschrieben. 
Die Existenz eines krummlinigen orthogonalen Netzes ist hierzu nicht erforder- 
lich. Die Betrachtungen lassen sich genau so im kartesischen (¥, v)- bzw. (uw, v)- 
System durchfiihren, und man kommt zum selben Ergebnis. 


_ (Eingegangen: 2. Oktober 1958.) 


4 1) Z. angew. Math. Phys. 9b, 764-776 (1958). 
2) Institut fir theoretische Gasdynamik der DVL. 
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Friihjahrstagung der Schweizerischen Physikalischen Gesellschaft 
vom 2. und 3. Mai 1959 in Brugg 


Berichte iiber die Sitzungen fiir angewandte Physik und Mathematik 


Hochempfindlicher Transistorenverstarker mit einem Eingangs- 
widerstand von mehr als 5 M®. Von H. June, Ziirich’). 
Eine ausfiihrliche Veréffentlichung erfolgt demnachst in der ZAMP. 


Koinzidenzschaltung mit Schwingkreis. Von H. Guur, B. LopKowicz 
und P. Marnier, Ziirich?). 

Wird ein gleichstromdurchflossener Schwingkreis plotzlich aus dem Stromkreis 
abgeschaltet, so entsteht am Schwingkreis eine exponentiell abklingende Sinus- 
schwingung («ringing circuit»’), welche sich angenahert wie folgt darstellen lasst : 


Tne /z gr so t. 


Die Anfangsamplitude ist proportional zum Gleichstrom J, der durch den Schwing- 
kreis geflossen ist. Bei konstantem Verhaltnis L:C variiert die Anfangsamplitude 
fiir Schwingkreise mit verschiedenen Resonanzfrequenzen nicht. Figur 1 zeigt das 


Figur 1 
Prinzipschema einer Koinzidenzmischstufe mit Schwingkreis. 1. S; + S, gleichzeitig gedffnet: 
V(t) = (1y + Ig) A(t). 2. S; oder S, allein gedffnet: Stark gedampfte Schwingung (wenn R <7). 
3. S. nach S, gedffnet: V(t) = I, A(t). 


Prinzipschema einer Koinzidenzmischstufe mit Schwingkreis. Von den drei be- 
schriebenen Moglichkeiten, die zwei Schalter zu 6ffnen, lést nur das gleichzeitige 
Offnen der beiden Schalter eine Schwingung aus, deren Anfangsamplitude minde- 


1) ETH. Physikalisches Institut. 
*) ETH. Cyclotron-Laboratorium. 


8) G, R. GAMERTSFELDER and J. V. Hotpam, Waveforms (Mc Graw-Hill, 1949), S. 141. 
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stens doppelt so gross ist wie die Anfangsamplituden der Schwingungen der Nicht- 
koinzidenzfalle. Die Koinzidenzauflésungszeit 2 t wird durch die Resonanzfrequenz 
des Schwingkreises bestimmt und betragt rund 1/2 Periodendauer. 

Figur 2 zeigt das Schaltschema einer Koinzidenzschaltung mit Schwingkreis. 


Die beiden Pentoden E180 F haben die Funktion der Schalter und der Damp- 


E1B0F E 180F 


|7kQ 
1002 |] 


Figur 2 


Koinzidenzschaltung mit Schwingkreis. 


fungswiderstande FR der in Figur 1 gezeigten Schaltung. Der Dampfungswiderstand 
hat bei stromdurchflossener Pentode angendhert einen Wert von der Grdéssen- 
ordnung 1/S (S Steilheit der Pentode). Wird die Streukapazitat durch kurze Zu- 
fiihrungsleitungen zum Schwingkreis klein gehalten, so lassen sich in dieser Schal- 
tung konventionelle Schwingkreise verwenden mit einer Resonanzfrequenz bis 
gegen 100 MHz. Damit die Eingangsréhren méglichst rasch geschlossen werden, 
empfiehlt sich die Verwendung von « White-Kathodenfolgern» als Treiberstufen zu 
dieser Koinzidenzschaltung. Vom Schwingkreis gelangt das Signal iiber einen 
Kathodenfolger zur Amplitudenmodulation. Im Koinzidenzfall erscheint am Aus- 
gang der Schaltung ein Impuls, welcher die doppelte Amplitude besitzt gegentiber 
Impulsen von Nichtkoinzidenzfallen. Ein nachfolgender Amplitudendiskriminator 
hat die Aufgabe, nur Koinzidenzen weiterzuleiten. 


Nx in % 


7 
1 
t 
x 

2 


Ke ee we en exe 


Delay in0™’s 
-0-8-6-4-2 0 2 4 6 8 10 
Figur 3 
Anzahl der Koinzidenzen, die bei gegenseitiger zeitlicher Verschiebung zweier Eingangsimpulse 


erhalten werden fiir zwei Schwingkreise mit den Resonanzfrequenzen von 20 und 50 MHz. 
——0-— Schwingkreis f = 50 MHzj37-— = X — = Schwingkreis f = 20 MHz. 
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Die Messung der Kurven der Figur 3 erfolgte mit Hilfe eines Quecksilbereich- 
impulsgebers. Je kiirzer die Periodendauer des Schwingkreises ist, desto steiler 
werden die Flanken der Koinzidenzauflésungskurve. Beim Schwingkreis mit der 
Resonanzfrequenz 50 MHz erstreckt sich diese Flanke tiber eine Zeit von 2 - 1071s. 

Zwei Vorteile dieser Schaltung bestehen darin, dass exponentiell abklingende 
Eingangsimpulse nicht zuerst mit Verzégerungskabel in kurze Impulse umgewan- 
delt werden miissen und dass ein Koinzidenzereignis durch eine lang andauernde 
Schwingung gemeldet wird. 


Schweissen mit Elektronenstrahlen. Von E. B. Bas und G. CREMOSNIk, 
Zurich’). 

In den Elektronenstrahlen kann sehr hoher Energiefluss in Form der kinetischen 
Energie der Elektronen erzeugt und auf einen kleinen Bereich einer Festk6rperober- 
flache konzentriert werden. Diese Warmekonzentration fiihrt zu einer 6rtlichen 
Schmelzung des Festkérpers, die schon seit langem, vor allem in den Rontgen- 
rohren, als eine gefiirchtete Erscheinung bekannt ist. In der letzten Zeit hat man 
diesen Effekt an vielen Orten fiir die Verschweissung von zwei Korpern niitzlich 
auszuwerten versucht®-$). 

Die einfachste Anordnung fiir Schweissen mit Elektronen ist die der Roéntgen- 
rohre; das heisst, das zu schweissende Werkstiick wird in einem Rontgenrodhren- 
system als Anode verwendet, wobei durch geeignete Bewegung des Werkstiickes die 
Schweissnaht durch den Brennfleck des Elektronenstrahles gefiihrt wird. Diese 
von WyMANN?®) und Srour®’) angewendete Einkammeranordnung hat trotz ihrer 
Einfachheit verschiedene Nachteile. Wir sind deswegen zu einer Zweikammeran- 
ordnung iibergegangen, welche in Figur 1 schematisch dargestellt ist. Die bedeutendste 
Autgabe bei einem Elektronenschweissgerat ist die Erzeugung der Elektronenstrahl- 
leistung und der Leistungsdichte im Schweissfleck. Diese Aufgabe ist besonders 
kritisch bei einem Zweikammersystem. Eine einfache Abschatzung der Schweiss- 
leistung fiir Materialwanddicken von einigen Millimetern fiihrt auf eine bereitzu- 


HU-Pumpe 


Werkstlick. : H ee g 
H V-Pumpe Oe 


Figur 1 
Schematische Darstellung des Zweikammer-Elektronenschweissgerites. 


i i 


——— 


STH. Abteilung fiir industrielle Forschung (AFIF). 
) W. L. Wyman, Welding Journal 37, 41 (Februar 1958). 
) J. A. Stour, Le Vide 13, 163 (Mai-Juni 1958). 
7 ‘ 4 
? a fs pte Proc. Ist. International Congress on Vacuum-Technique 1958 (Pergamon Press). 
) E. B. Bas, Proc. 1st. International Congress on Vacuum-Technique 1958 (Pergamon Press). 
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Figur 2 
Hochstromelektronenkanone mit Bolzenkathode: 1 Wolframbolzen; 2 Einspannklemme fiir Wolf- 


rambolzen; 3 Kathode (Tantal, Renium, Wolfram); 4 Wolfram-Heizwendel; 5, 6 Haltestiften; 
7 Strahlungsschutzschirme; 8 Kathodenblende; 9 Anode. 
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Figur 3 
Raumladungscharakteristik der Hochstromelektronenkanone mit Bolzenkathode (Durchmesser der 
Emissionsflache 1,55 mm). 
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60 mm o 


Figur 4 
Stumpfgeschweisste Cr-Ni-Stahlrohre (Aufsicht). 


Figur 5 
Querschnitt der Schweissnaht in Cr-Ni-Stahl. 


Figur 6 
Schweissnaht im Kupferrohr; Rohrdurchmesser 60 mm, Wanddicke 2,5 mm, Schweissleistung 1 kW 
(vorgewarmt), Vorschub 2 mm/s. 
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stellende Strahlleistung von einigen Kilowatt; das-heisst, bei einer Beschleunigungs- 
spannung von 20 kV muss die Elektronenkanone Stréme von 100 mA und mehr. 
liefern kénnen. Diese Str6me miissen auf einen Bereich von etwa 1mm Durch- 
messer fokussierbar sein. In unserem Gerat wurde eine Elektronenkanone mit gross- 
flachiger Bolzenkathode®) angewendet. Diese ist in Figur 2 schematisch dargestellt. 
Die Emissionsstromcharakteristik dieser Elektronenkanone ist in Figur 3 wieder- 
gegeben. Wie ersichtlich, ist die Elektronenkanone in der Lage, bei einem Durch- 
messer der Emissionsflache von nur 1,5mm Stréme von 100 mA und mehr zu liefern. 

In Figur 4 ist eine Schweissnaht in Cr-Ni-Stahl 18/8 zu sehen. Ein Querschliff 
durch diese Schweissnaht ist in Figur 5 dargestellt. Cr-Ni-Stahl benétigt sehr 
geringe Schweissleistung (kleine Warmeleitfahigkeit). Bei der Schweissung in Fi- 
gur 4 betrug diese nur 150 W bei einer Vorschubsgeschwindigkeit von etwa 4 mm/s. 
Kupfer benétigt nach Wolfram die grésste Schweissleistung. Figur 6 zeigt die 
Schweissnaht zwei aufeinander stumpfgeschweisster Kupferrohre von 60 mm 
Durchmesser und 2,5mm Wanddicke. Figur 6a stellt den Querschliff durch die 
Schweissnaht dar. In Figur 6) ist die Schweissnaht geatzt und in Aufsicht abge- 
bildet. Der Vorschub erfolgte von rechts nach links (im Aufsichtsbild) und betrug 
2 mm/s bei einer Schweissleistung von ungefahr 1 kW. Das Schweissgut wurde bei 
defokussiertem Strahl vorgewarmt. 

Ein ausfiihrlicher Bericht erfolgt demnachst an einer andern Stelle. 


Uber eine Fernfokus-Elektronenkanone und ihre Anwendungen. Von 
E. B. Bas und F. Gaypovu, Ziirich?*). 


Zusammenfassung 


Es wurde ein Strahlerzeugungssystem (Figur 1), bestehend aus einer Aquipoten- 
tiel-Bolzenkathode®) mit planer kreisf6rmiger Emissionsflache (A), einer Wehnelt- 
Blende besonderer Geometrie (W) und einer Anode (A) auf seine fernfokussierende 
Eigenschaft untersucht. Bei optimaler Wahl der Geometrieparameter kann man den 
Elektronenstrahl in 150 mm Entfernung von der Kathode in einem Brennfleck von 


Automatische Zentrierung des Elektronenstrahlers auf die Spitze eines dachkant- oder kegelf6r- 

_ migen Targets. K Kathode; W Wehnelt-Blende; A Anode; T Dachkant- oder Kegeltarget; Si, So 

MeBsonden fiir Sekundarelektronen; B Rohrenbriicke fiir die Steuerung des magnetischen Zentrier- 
feldes; M magnetisches Ablenkelement. 


9) E. B. Bas, Z. angew. Phys. 7, 337 (1955). 
10) ETH, Institut fiir technische Physik. 
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150 u Halbwertsdurchmesser fokussieren. Dabei kénnen Strahlstrome von einigen” 
Milliampere bei einer Beschleunigungsspannung von 50-60 kV erreicht werden. Der 
kritischste Parameter ist die Lage der Vorderbohrung der Wehnelt-Blende und 
deren Durchmesser. Das Wegriicken der Anode von der Kathode vermindert zwar 
auch den Strahldurchmesser im Fokus, es wird aber dabei auch die Feldstarke vor 
der Kathode reduziert, was sich auf die Raumladungsverhaltnisse ungtinstig aus- 
wirkt. Je nachdem der Strahldurchmesser im Fokus grésser oder kleiner als der 
Kathodendurchmesser ist, erscheint der Elektronenstrahl als ein divergierendes 
oder auf den Fokuspunkt konvergierendes Strahlenbiindel. 

In einer besonderen, in Figur 1 dargestellten Anordnung wurde der Elektronen- 
strahl auf die Spitze eines Kegeltargets fernfokussiert. Der Brennfleck von etwa 
150 wu Durchmesser fallt gewohnlich nicht zentrisch auf die Spitze des Targets. Um 
einen zentrierten Betrieb zu erméglichen, wurde folgender Weg eingeschlagen: Die 
von dem Target ausgelésten Sekundarelektronen werden von zwei MeSsonden S, 
und S, aufgefangen und die Stréme J, und J, einer Rohrenbriicke zugeftihrt, welche 
ein magnetisches Ablenkelement speist. Trifft der Strahl die Spitze nicht zentrisch, 
so ist J, + J,, und durch die Ablenkspule fliesst ein Strom, welcher den Strahl auf 
die Spitze ablenkt, bis der Strahl wieder die Spitze zentrisch trifft und bis auf eine 
kleine Regelrestgrésse 1, = I, wird. Die Selbstzentrierung des Strahles eliminiert 
auch die Einfliisse der Storfelder. Fiir die Zentrierung auf ein Kegeltarget ist die an- 
gegebene Schaltung fiir zwei aufeinander senkrecht stehende Richtungen auszufihren. 

Die beschriebene Fernfokus-Elektronenkanone mit Selbstzentrierung erméglicht 
uns die Realisierung einer neuartigen Weitwinkelrontgenréhre fir Panorama- 
Gebissaufnahmen in der Dentaldiagnostik™). Die prinzipielle Anordnung dieser 
neuen Methode der Radiographie ist in Figur 2 dargestellt. 

Ein ausfiihrlicher Bericht iiber die Untersuchungen erfolgt demnadchst an einem 
andern Ort. 


Figur 2 
Eine Weitwinkel-R6ntgenréhre fiir Dental-Panoram 
3 kegelférmiges Target (Antikathode) 
satzes der Rohre; 5 Film; 6 Filmtrager; 7 vakuumaufg 


aaufnahmen. 1 Réntgenrdhre; 2 Oltank; 
; 4 Aluminium-Abschirmrohr fiir den Schutz des Glasfort- 
d edampfte MeBsonden fiir Sekundarelektronen 
auf der Innenseite des Glasfortsatzes; 8 magnetisches Zentrierquadrupol. 

1) E.B. Bas und F.Gaypou: Proc. Ist. In 
(Pergamon Press). 
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Uber die relativistische Elektronenoptik eiektrostatischer Beschleu- 
niger. Von E. B. Bas, L. Preuss und W. ScHNeEIpER, Ziirich!), 

Ein elektrostatischer Partikelbeschleuniger besteht aus einer Kette von Be- 
schleunigungsstufen, zum Beispiel in der-Form von Lochblenden- oder Rohr-Im- 
mersionslinsen. Die erste Aufgabe-dieser Linsenkette ist die Beschleunigung der 
mit einer Energie U, eingeschossenen Partikel auf eine Endenergie Up. Sie muss 
auch eine zweite Aufgabe erfiillen, die darin besteht, den engsten Querschnitt P 
des Einschu8strahles auf ein Target Q am Ausgang des Beschleunigers abzubilden. 
Die Berechnung der Elektronenoptik des Beschleunigers muss also eine Beziehung 
zwischen den Achskoordinaten 2, und 2, liefern, und dabei ist noch die Lateralver- 
grosserung f und die Angularvergrésserung y anzugeben. 

Fir die achsennahen Elektronenbahnen in einem elektrischen Feld mit dem 
Achsenpotential @(z) gilt allgemein folgende Differentialgleichung : 


Sy eer ae 
a & = — aches ig =, = 4 pe — 
Vi+eo) b © Vare®) | tg +209) o" 0 
mit 
ee ad 
eal EAE 


Die Multiplikatoren in runden Klammern beriicksichtigen den relativistischen 
Effekt. Fiir ein Potential d = 5-106 V erhalt man zum Beispiel 


fiir Elektronen (1+ ¢@) = 5,88, 
fiir Protonen (1 + e¢ @) = 1,0027. 
Wir ersehen daraus, dass fiir lonenbeschleuniger der relativistische Effekt in den 
elektrostatischen Beschleunigern nicht beriicksichtigt werden muss. Hingegen bei 


Elektronenbeschleunigern ist es angebracht, sich iiber den relativistischen Effekt 
Rechenschaft zu geben. 


| 
45cm-—=— 
EF 405m ~— 20000 
605m 


U, 
228. 
4 1 hs 


Figur 1 
Schema eines elektrostatischen Beschleunigers. 


; Wir haben diese Rechnung fiir einen 10stufigen Elektronenbeschleuniger mit 
_ der Endenergie von Ug = 1,5- 10° V durchgefiihrt. Die Geometrie der Beschleuni- 
gungsstrecke ist in Figur 1 dargestellt. Die einzelnen Beschleunigungsstufen sind 
- Zweirohr-Immersionslinsen mit Spaltbreite S = 4cm, Rohrradius R = 5 cm. Die 
Potentialverteilung in solchen Immersionslinsen wurde in einem Widerstands- 


ar 
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12) ETH. Institut fiir technische Physik. 
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Netzwerkmodell ausgemessen. Die Berechnung der Kardinalelemente der Linsen 
erfolete mit der elektronischen Rechenmaschine ERMETH. Die Rechnungsergebnisse 
fiir de in Figur 1 angegebenen Daten sind in Figur 2 zusammengestellt. Es sind 


0 
0 x=Uy/U, ———— se Q5 10 
0 U, kv] 50 100 150 
Figur 2 
10 Stufen; 1500 kV; S/R = 0,8; R =5cm; —relativistisch; -—-—-— nichtrelativistisch. 


dort in Abhangigkeit von der Einschussenergie U, die Lage der Elektronenquelle z, 

(bezogen auf die Mittelebene der ersten Linse), die laterale Vergrésserung der Elek- 

tronenquelle auf dem Target # und der Winkelvergrésserung y angegeben, wobei 

die Rechnung sowohl relativistisch wie nichtrelativistisch durchgefiihrt wurde. 
Ein ausfiihrlicher Bericht erfolgt demnachst in der ZAMP. 


- 


Reaktoroszillator zur Bestimmung makroskopischer Absorptions- 
querschnitte von thermischen Neutronen. Von T. HURLIMANN, P. SCHMID 
und H. WINKLER, Wiirenlingen?’). 

Am Swimming-Pool-Reaktor SAPHIR ist ein Reaktoroszillator eingerichtet 
worden. Er wird hauptsachlich fiir Reinheitsuntersuchungen von Materialien ein- 
gesetzt. 

Figur 1 gibt eine schematische Darstellung der Einrichtung. Die gewahlte Spalt- 
stoffanordnung zeigt an der Messposition ein Maximum des Produktes aus thermi- 
schem Fluss und thermischer Einflussfunktion. Dort wird das die Probe enthaltende 
Schiffchen periodisch gegen ein zweites Schiffchen ausgetauscht, das lediglich als 
«Absorptionsgegengewicht» dient. Ein solches Gegengewicht ist bei wassermoderier- 
ten Reaktoren vorteilhaft, da eine einzelne Probe infolge ihrer Wasserverdrangung 
grosse Absorptions- und Moderationssignale hervorruft. Durch ein Auswechseln 
der Probe an der MeB8stelle gegen eine ahnliche Probe kann das Bruttosignal 
; 


F 


13) Reaktor AG, 
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relativ klein gehalten werden, so dass die Apparatur in einem empfindlicheren 
Bereich arbeiten kann. Bei geeigneter Wahl des gegenseitigen Abstandes der beiden 
Schiffchen wird ausserdem erreicht, dass das wiihrend des Austauschvorganges 
hervorgerufene kombinierte Transfersignal auf den 50. Teil des einzelnen Transfer- 
signals sinkt und daher die Messung nicht mehr stért. Allerdings sind dann die 
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Figur 1 


Anordnung des Reaktoroszillators. Das obere Probenbiichslein befindet sich eben in der 
Messposition. 


Proben in der Position ausserhalb des Reaktors nicht mehr ganz ohne Einfluss auf 
die Reaktivitat des Reaktors. Aus diesem Grunde wird das untere Schiffchen 
nicht gerade mit dem fiir den Bezug verwendeten Standard beladen. 

Die Messproben bestehen gewohnlich aus zwei Scheiben von 26mm Durch- 


_ messer und 12 mm Hohe. Sie werden jeweils wahrend 19 s in einer Position gehalten. 


Der Austauschvorgang dauert eine ganze Sekunde, wird aber dafiir sehr prazis 
sinusf6rmig und ohne Phasenstreuung durchgefiihrt. Ein Messzyklus besteht aus 
10 Halbperioden und ist deshalb in 200 s beendet. Die erzeugte Leistungsmodulation 
AN wird in Einheiten der mittleren Reaktorleistung Ny gemessen, wie auch das 
Signal eines in einem weiteren Messzyklus oszillierten bekannten Standards. Der 
Bezug dieser Einheit auf den Absorptionsquerschnitt erfolgt mittels eines Eich- 


" satzes. Er besteht aus 3 lithiumhaltigen Aluminiumscheibchen, welche zwischen die 


Probenhalften bzw. dariiber und darunter eingebracht werden. Die relativen 


_ Lithiumgewichte (3,64: 1:1) sind so gewahlt, dass der Eichsatz wirkt, als ob das 
_ Messgut moglichst gleichmassig mit Lithium versetzt ware. Diese Wirkung wird 
~ dann vollkommen erreicht, wenn sich die allgemeine axiale Abhangigkeit des Pro- 


duktes von thermischem Fluss und Einfluss innerhalb der Probe durch ein Polynom 


- vierten Grades ausdriicken lasst. Diese Eigenschaft des Eichsatzes erlaubt die Mes- 
_ sung von Absorptionsquerschnitten in Proben, die merkliche Flussanderungen im 
- Innern aufweisen. Insbesondere ist es méglich, bedeutend gréssere Konzentrationen 


‘\% 
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yon den zu untersuchenden Substanzen in Losungen zu dulden. Bei normalen 
Messungen ist der Eichsatz durch einen Nullsatz aus Reinaluminium ersetzt. Die 
Differenz der Absorptionsquerschnitte von Eich- und Nullsatz ist auf denjenigen 
von Bor (Standard aus Harwell) bezogen. 
Der Verlauf der Reaktorleistung wird mit einer Ionisationskammer verfolgt. Die 
Auswertung des Kammerstromes geschieht mit der in Figur 2 dargestellten Appa- 


Demodulator Integrator X 


DC-Verstarker 
A U; 


lonisations- 
kammer 


Fichen Demodulator ‘Integrator Y 
7 


RovMessen 97, 9 
Figur 2 
Messapparatur. Der Kammerstrom i(t) = i) + Ai entspricht der Reaktorleistung N(t) = No + AN. 


Kompensationsspannung: y U, = R ig. Verstarktes Modulationssignal: U, = R Ai/a«. Im Integrator 
X wird das Signal registriert. Der Integrator Y dient der Kontrolle (siehe Text). 


ratur. Sie enthalt Integratoren mit unbeschrankter Integrationskapazitat, prakti- 
scher Tragheitslosigkeit und ausserordentlich guter zeitlicher Stabilitat (Drift pro 
Stunde < 2-10-§ des maximalen Eingangssignals). Die Eichung des Signals in 
Einheiten von N, erfasst Ungenauigkeiten und Drifts der gesamten elektronischen ~ 
Apparatur mit Ausnahme der Integratoren und der eigentlichen Eichwiderstande. 
Die rechteckige Demodulation des Reaktorsignals (Figur 3) erfolgt im einen Schal- 
tungszweig phasenrichtig, so dass im Integrator X das Signal gemessen wird. In 
einem zweiten Schaltungszweig geschieht die gleiche Demodulierung um 90° pha- 
senverschoben, so dass bei einwandfreiem Funktionieren im Integrator Y nur das 
Rauschen registriert wird. 

Um den Einfluss von Reaktorleistungsdrift auf das Messergebnis méglichst zu 
eliminieren, wird die Demodulation der 10 Halbperioden eines Messzyklus mit ver- 


schiedenen Gewichten g durchgefiihrt. Falls der Binomialkoeffizient ( pace als Ge- 


wicht der k-ten Halbperiode gewahlt wird, werden die ersten neun Glieder der 
Potenzreihe genau kompensiert, in welche sich der Reaktordrift entwickeln lasst. 
Auf diese Weise ist das Messergebnis in praktisch allen Fallen von langsamen 
Leistungsanderungen ungestort. 

_ Die mittlere Schwankung von Absorptionsmessungen, die bei einer Reaktor- 
leistung von 1 W durchgefiihrt werden, betragt 0,1 mm?. Dieser Wert entspricht 
genau dem erwarteten Beitrag des Reaktorrauschens. Die Ionisationskammer ist 
mit ihrer Mefistatistik dank der giinstigen Plazierung im Swimming-Pool-Reaktor — 
und der hohen Neutronenempfindlichkeit ohne Einfluss auf das gemessene Rau- 
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schen™). Die Messgenauigkeit entspricht somit dem theoretisch optimalen Wert. 
Sie kann natiirlich noch weiter verbessert werden durch mehrfache Wiederholung 
der Messung oder durch Erhéhen der Reaktorleistung. Die Reproduzierbarkeit 


xan jeweiligem Probenwechsel ist von der gleichen Grésse wie die mittlere Schwan- 
ung. 


Reaktorleistung Nt) 


0 16 Zz 
Figur 3 


Gp 


Leistungs- und Demodulationssignal. Das Reaktorleistungssignal entspricht einem Absorptions- 
querschnitt von 4 mm?. 


Die Permeabilitat von Pulverbetten. Von W. WIELAND”). 
Eine ausfiihrliche Publikation erscheint demnachst in der ZAMP. 


Isotherme Strémungen in Raketendiisen. Von V. STINGELIN, Ziirich"*). 

Die Aufgabe soll sein, mit einer gegebenen Warmequelle (Reaktor), deren 
Héchsttemperatur beschrankt ist, einem Gase eine extrem hohe Austrittsgeschwin- 
digkeit zu erteilen. 

Bei den auf chemischer Basis arbeitenden Raketendiisen ist der maximalen 
Ofenenthalpie durch die Verbrennungsenergie des Treibstoffes zum voraus eine 
obere Grenze gesetzt. Haben wir einen Reaktor zur Verfiigung, so kénnen wir 
prinzipiell den gleichen Prozess durchfiihren, nur dass sich jetzt das Gas von aussen 
her, durch Warmeiibertragung an der Wand, erhitzt. Die maximal zuldssige Reaktor- 


_ temperatur begrenzt auch hier die Austrittsgeschwindigkeit. 


= 


RAY 


Aus diesen Erwagungen haben wir die alte Idee der isothermen Expansion 
(Zwicky, WINTERBERG) wieder aufgegriffen und naher untersucht. 

14) O. R. Friscu und D. J. Litrrer, Phil. Mag. 45, 126 (1954). 

15) Zementfabrik Holderbank-Wildegg AG, technische Stelle. 

16) ETH. Institut fiir Aerodynamik. 
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Unter der vereinfachenden Annahme von konstanter Wandtemperatur erwar- 
men wir das Gas anfanglich isobar, bis die Differenz zwischen. Wand- und mitt- 
lerer Gastemperatur klein ist. Unter weiterem Autheizen expandieren wir anschlies- 
send das Gas derart, dass seine mittlere Temperatur konstant bleibt. Die Enthalpie 
des Mediums gibt uns zum Schluss durch adiabatische Expansion auf den Aussen- 
druck nochmals einen Geschwindigkeitszuwachs. Im leeren Raum, mit dem Druck0, 
kénnen wir theoretisch auf diese Weise beliebig viel Energie zufiihren und somit 
die gewiinschten grossen Geschwindigkeiten erreichen. Selbst mit Reibung, welche 
ja bei realen Strémungen immer auftritt, ware der Prozess realisierbar, allerdings 
mit entsprechenden grésseren Dimensionen der Diisen (siehe Figur). 


< 
wa 
5 2 PLE 
0, 
: 
b fe ? 


ANANANAAAANNNY 


Figur 1 


Zustandsdiagramm fiir konventionelle Erwarmungs- und Expansionsprozesse: a Gasgemisch durch 
Verbrennung geheizt; 6 Gas durch Reaktor geheizt. 


‘Figur 2 
Prinzipskizze (nicht maBstablich) und T-S-Diagramm der isothermen Expansion. 


Viel schlimmer wirkt der Temperaturriickgewinn in der Grenzschicht. Die fiir 
den Warmeiibergang massgebende Temperaturdifferenz AT bleibt nimlich mit zu-_ 
nehmender Mach-Zahl nicht konstant. Unter der Beriicksichtigung der Kompres- 
sibilitat kénnen wir fiir die laufende Temperaturdifferenz schreiben: ; 

aes 2 
Ow = 0, —Y% ‘i 1 . Os 35 
y Zz 
0. die durch die mittlere Gastemperatur dividierte laufende Temperaturdifferenz 
zwischen — Wand und Grenzschicht; 
o =w//RT,, I, mittlere Gastemperatur, R Gaskonstante, y = Coleus 
vy  Temperaturriickgewinnfaktor (0,85 fiir laminar, 0,9 fiir turbulent). 
Der Index 0 bedeutet am Eintritt der Diise. 


Reg 5 eA > ey ewe ee 
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Daraus berechnet sich sofort die maximale isotherme Geschwindigkeit bzw. 
Omax DENN © = Mygy, WENN Ow = 0, also 


ae 
Omax = \/2 a : y 5 : | 


Im schlimmsten Fall, y = 1, wird die totale Geschwindigkeit, welche wir mit unserem 
Prozess erhalten, 


» . rc » rag 
pop 1S 2 
Cot = |/ Wimax a Wad 


max 


wu 2% 2, 


und mit 


FAN == 
Wot =|/2R7, ae 


max 


Dies ist aber die maximale Geschwindigkeit eines Gases, das von der Temperatur 
T,, (Wandtemperatur) auf den Druck 0 adiabatisch expandiert wird. 

Wir erhalten also im Extremfall bei unendlich langey Diise nicht mehr, als wenn 
wiry im unendlich langen Rohr isobar erwarmen und anschliessend adiabatisch auf den 
Druck 0 expandieren. 


_Ein optischer Tiefentaster. Von W. Lormar, Aarau’). 
Bei der dreidimensionalen Ausmessung von K6rpern wie zum Beispiel Turbinen- 
schaufeln oder Windkanalmodellen mit Hilfe von Koordinatenmessmaschinen be- 
steht die Aufgabe, die Tiefenlage von Objektpunkten zu ermitteln, welche im all- 
gemeinen auf schief zu den Koordinatenachsen stehenden Flachen liegen. Da die 
mechanische Antastung auf solchen Flachen wegen des endlichen Radius der Tast- 
spitze prinzipiell mit Messfehlern behaftet ist, wird vielfach ein Mikroskop als 
optischer Tiefentaster verwendet. Man stellt dabei an einer durch eine Marke be- 
-zeichneten Stelle des Gesichtsfeldes das Objekt auf beste Scharfe ein. Um eine genii-’ 
gende Messgenauigkeit zu erreichen, ist die Anwendung einer relativ starken Ver- 
grésserung (etwa 200mal) nétig, was wiederum eine kraftige Lichtquelle fiir die 
Auflichtbeleuchtung bedingt. Wegen der Unsicherheit bei der Einstellung auf 
«beste Scharfe» sowie der thermischen Einwirkungen einer starken Lichtquelle auf 
die Messmaschine befriedigt jedoch diese Anordnung nicht ganz. 
Es ist bekannt, dass man mit einem Mikroskop viel genauer auf die Lage einer 
Messmarke quer zur optischen Achse einstellen kann als auf ihre Lage in Richtung 
der Achse. Die Suche nach einer Methode, bei welcher die Anderung der Tiefen- 
koordinate eine seitliche Bewegung einer Messmarke bewirkt, fiihrte auf die in der 
‘Figur gezeigte Anordnung, welche im Prinzip dem Lichtschnittmikroskop von 
‘ScuMatz ahnlich ist. Ein unter etwa 30° zur Messtichtung stehendes Mikroskop 
projiziert auf das Objekt eine Marke, deren Lage durch ein symmetrisch dazu an- 
-geordnetes zweites Mikroskop beobachtet wird. Bei Anderung der Tiefenkoordinate 
-wandert die Marke durch das Gesichtsfeld und kann mit einer im Betrachtungs- 


17) Kern & Co. AG. 
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mikroskop befindlichen festen Marke zur Deckung gebracht werden. Mit der Ver- 
schiebung der Marke im Gesichtsfeld ist allerdings auch eine Anderung ihrer 
Scharfe verbunden, doch lasst sich die Anordnung natiirlich so justieren, dass mit 
Einstellung der Koinzidenz auch die beste Scharfe erreicht wird. Durch Anbringung 
einer Skala im Ableseokular kann die Anordnung ferner als Mikrometer ausgebildet 
werden. 


Messrichtung 


Messungen an einer Versuchsausfiihrung!*) haben ergeben, dass bereits mit 
50facher Vergrésserung eine Messgenauigkeit von + 1p erreicht wird. Bei Verwen- 
dung von Objektiven mit einer numerischen Apertur von etwa 0,15 geniigt dann 
bereits eine kleine Gliihlampe von einigen Watt Leistung zur Beleuchtung vollauf. 

Das Okularbild kann auch auf eine Mattscheibe projiziert werden, doch ist dann 
eine kraftigere Lichtquelle erforderlich. 
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The Hypercircle in Mathematical Physics. A Method for the Approximate 
Solution of Boundary Value Problems. Von J. L. Syncr (Cambridge University 
Press,)1957): 424.5. 134 Ris.; 70s4 > 

In den letzten zehn Jahren sind von SyNGE und anderen (J. B. Diaz, W. PRAGER, 
A. WEINSTEIN) eine ganze Reihe von Arbeiten erschienen, welche bei gewissen 
Randwertaufgaben der mathematischen Physik, vor allem bei Problemen aus der 
Mechanik, mit Hilfe der Hypercircle-Methode Fehlerabschatzungen gaben. Es ist 


sehr zu begriissen, dass der ganze Ideenkreis im vorliegenden Buche nun eine ge- 
schlossene Darstellung gefunden hat. ' 


18) Ausgefiihrt von der Société Genevoise d’Instruments de Physique. 
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Viele der angedeuteten Randwertaufgaben lassen sich als Extremalprobleme 
formulieren, das heisst, bei der Interpretation in einém Funktionalraume entspricht 
die Lésung demjenigen Punkt in einer Ebene mit dem kleinsten Abstand von einem 
festen Punkt ausserhalb der Ebene, oder,-was dasselbe ist, dem Schnittpunkt von 
zwei total senkrechten Ebenen. Bei-einem vorliegenden Problem dieser Art gilt es 
nun zuerst, die Aufspaltung in die beiden senkrechten Teilraume iiberhaupt zu 
finden. Sodann kann man aber, mit Hilfe von elementargeometrischen Betrach- 
tungen, die Lésung auf einem Hyperkreis (= Schnitt Hyperkugel-Hyperebene) von 
bekanntem Zentrum und Radius lokalisieren und damit zum Beispiel, falls das 
Zentrum als Naherungslésung verwendet wird, Abschatzungen fiir den Fehler im 
Sinne der mittleren Quadrate gewinnen. Diese Abschatzungen sind, auch wenn das 
Problem mit einem Differenzenverfahren behandelt wird, theoretisch exakt. 

Es werden keine Vorkenntnisse aus der Funktionalanalysis vorausgesetzt, viel- 
mehr wird auch der Nichtmathematiker anhand von instruktiven Beispielen und 
Ubungsaufgaben griindlich in die Geometrie des Hilbert-Raumes eingeftihrt. Die 
praktische Durchfiihrung der Methode wird im Falle der positiv definiten Metrik 
am Beispiel des Dirichlet- und Neumann-Problems (vor allem fiir verschiedene 
Torsionsaufgaben) ausfiihrlich demonstriert. Im Falle der indefiniten Metrik werden 
Schwingungsprobleme behandelt. 

P. LAUCHLI 


Matrizen (2. Auflage). Von Dr.-Ing. R. ZuRMUHL (Springer-Verlag, Berlin 
1958). 467 S., 76 Abb.; DM 33.-. 

In grossen Ziigen ist der Inhalt der zweiten Auflage gleich geblieben, aber der 
methodische Aufbau der einzelnen Kapitel — insbesondere im theoretischen Teil — 
wurde gegeniiber der ersten Auflage bedeutend verbessert. Insbesondere wirkt sich 
vorteilhaft aus, dass die Theorie der linearen Gleichungen (Kapitel II) jetzt auf dem 
GauBschen Algorithmus aufgebaut wird; der Rangbegriff ergibt sich damit in 
zwangloser Weise. Kapitel VI (numerische Methoden) ist durch Aufnahme des Ver- 
fahrens von STIEFEL-HESTENES (konjugierte Gradienten) bereichert worden, da- . 
gegen verbleibt der Abschnitt iiber Eigenwertprobleme im wesentlichen leider auf 
dem alten Stand von 1950, der heute als iiberholt bezeichnet werden muss. Beispiels- 
weise erlaubt keines der vom Verfasser propagierten Verfahren die saubere Bestim- 
mung aller Eigenwerte einer symmetrischen Matrix grosser Reihenzahl, obwohl es 
mehrere solche Methoden gibt. Das Kapitel VII (technische Anwendungen) wurde 
unter Aufnahme von Beispielen aus der Mechanik erweitert. 

H. RUTISHAUSER 


Light Scattering by Small Particles. Von H. C. van bE Huzst (John Wiley 
_ & Sons, Inc., 1957). 470 S., 103 Fig.; $ 12.00. 

Die Erfolge der Kolloidchemie, durch Heranziehung der Lichtstreuung Auf- 
schliisse iiber die Teilchenstruktur zu gewinnen, und in jiingster Zeit La MERs eé¢ al. 
Entdeckung eines Tyndall-Linienspektrums an monodispersen Aerosolsystemen 
haben dem Problem der Wechselwirkung zwischen Licht und dispersen Stoft- 
teilchen wieder vermehrte Aufmerksamkeit gebracht. Die Erfindung des Radars 
‘und insbesondere die neuen Méglichkeiten, die damit fiir die Meteorologie entstanden 

sind, belebten zusehends auch von der Seite der Mikrowellen her das Interesse 
an den Lichtstreuungsvorgangen. Die Bedeutung einer mathematischen Erfassung 
der Streuungserscheinungen ist insofern noch erheblich gesteigert worden, als auch 
die Kernphysik sich mit ahnlichen Vorgangen, wie zum Beispiel bei der Streuung 
von Neutronen, auseinanderzusetzen hat. 
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So ist ein Teilgebiet der theoretischen Physik, das den fiihrenden Forschern des. 
19. Jahrhunderts grosse Fortschritte zu verdanken hatte und zu denen im gegen- 
wartigen Jahrhundert nur noch die Arbeiten von DEBYE ausschlaggebend bei- 
steuerten, in jiingster Zeit wieder zu neuem Leben erwacht und hat, besonders ent- 
facht durch die Erfordernisse der Kernphysik, in der Entwicklung mathematischer 
Verfahrenstechniken erhebliche Friichte getragen. Es muss daher vAN DE HULsT 
sehr angerechnet werden, in einem umfassenden Werk die Grunderscheinungen 
der Lichtstreuung an kleinen Stoffteilchen zur Darstellung gebracht zu haben. 

Das Buch besteht aus den drei Hauptabschnitten: 

1. Allgemeine Stérungstheorie. Diese wird entwickelt, ohne auf besondere For- 
men der streuenden Teilchen Bezug zu nehmen. 

2. Spezielle Typen der Streuteilchen. Dieser Abschnitt beansprucht den grossten 
Teil des Werkumfanges. 

3. Anwendungen. Sie erstrecken sich auf ausgewahlte Gebiete der Chemie, 
Physik, Meteorologie und Astronomie. Im Falle der Sternkunde wird auch die Frage 
der Lichtstreuung an interstellarer Materie berihrt. 

Das Buch beschrankt sich auf die Behandlung der einfachen, kohadrenten 
Streuung an Teilchen, die voneinander unabhangig sind; trotzdem verliert es aber 
dadurch nicht an Bedeutung. Eines der Kernprobleme der Schrift bildet die mathe- 
matische Erfassung der Streuung an einer homogenen Kugel. Erfreulich ist der 
klare Aufbau und die konsequente Darstellungsweise. Jedes einzelne Kapitel schliesst 
mit einem Literaturverzeichnis, dem die einschlagigen Arbeiten entnommen werden 
kénnen. Auch die Aufmachung von seiten des Verlages verdient Anerkennung. 

R. SANGER 


Theorie der Relaissysteme der automatischen Regelung. Von Ja. S. 
ZyPxkIn (R. Oldenbourg, Miinchen; Verlag Technik, Berlin 1958). Ubersetzung aus 
dem Russischen: W. Haun, Braunschweig, und R. HERSCHEL, Ulm. 472S., 249 Abb., 
17 Tab.; DM 52.-. 

Relaissysteme nehmen in der Regelungstechnik einen bedeutenden Platz ein, doch 
sind die Eigenschaften der geschlossenen Systeme infolge der Nichtlinearitat wesent- 
lich schwieriger iiberblickbar als Regelsysteme mit linearen Gliedern. Dem Ver- 
fasser gelingt es, mit exakten Methoden die Wege anzugeben, auf welchen man unter 
Zuhilfenahme der Laplace-Transformation, bzw. mittels der Ubertragungsfunktion 
(Frequenzgang) oder der Ubergangsfunktion (im Buche mit Zeitcharakteristik be- 
zeichnet) fiir den linearen Teil des Regelkreises, zu den allgemeinen Stabilitats- 
kriterien gelangt. Der gesamte Regelkreis wird dabei aus einem linearen Teil und 
dem Relaisteil gebildet. Der Relaisteil kann eine Totzone, Hysterese, keines von 
beiden oder beide kombiniert, zwischen zwei festen Ausgangswerten in Funktion 
der Eingangsgrésse aufweisen. Zuerst wird immer der Fall ohne Totzone und ohne 
Hysterese untersucht. Die notwendigen Modifikationen fiir die ibrigen Falle sind 
dann relativ leicht ableitbar. 

Eine besondere Charakteristik von Relaissystemen ist das Auftreten von stabilen 
Schwingungen mit begrenzter Amplitude (zum Beispiel Temperaturregelung mit- 
telst Thermometerkontakt). Die Ubertragungseigenschaften des linearen Systems 
beeinflussen diese Erscheinung, so dass von der Aufgabenstellung her geeignete 
Massnahmen getroffen werden miissen, um eine Verbesserung der Regeleigenschaf- 
ten zu erreichen. Graphische (Phasendiagramme, Hodographen) und analytische 
Verfahren werden diskutiert und an Beispielen erlautert. Solche stabilen Schwin- 
gungen k6nnen vom System selbst herriihren oder von einer periodischen Eingangs- 
grosse erzwungen sein (Mitnahmeeffekt). Im allgemeinen Fall treten mehrere Eigen- 
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schwingungen auf, so dass ein periodisches Verhalten beobachtet werden kann. Die 
Untersuchung der Stabilitat solcher periodischen’Zustande braucht schon erheb- 
liche Mittel, mit welch»n besondere Stabilititskriterien herausgearbeitet werden. 
Schliesslich werden Massnahmen zur Linearisierung von Relaissystemen und dabei 
auftretende optimale Vorgange diskutiert. Ein mathematischer Anhang, Tabellen 
und ein ausfiihrliches Literaturverzeichnis, nach Autoren geordnet, in dem die rus- 
sischen Quellen stark iiberwiegen, beschliessen das Buch. 

Der Vorzug des Werkes liegt in der einheitlichen Verarbeitung des weit verstreuten 
Materials. Wer mit der Theorie der linearen Regelvorgange vertraut ist, wird mit 
Genuss dieser bei aller Konzentration sehr vollstandigen Darstellung der Regelung 
mittelst Relaissystemen folgen. Etwaige ungewohnte Fachausdriicke, durch die 
Ubersetzung bedingt, kann der Leser leicht selbst korrigieren. H. WERER 


Les calculateurs numériques automatiques et leurs applications. 
Von A. SESTIER (Editions «Hommes et Techniques», Neuilly 1958). 183S., 
fFr. 3000.-. 

Das Buch von SESTIER vermittelt einen allgemeinverstandlichen Uberblick iiber 
die digitalen Rechenautomaten. Sowohl deren charakteristische Merkmale wie 
auch ihre Anwendungsméglichkeiten in Industrie, Handel und Wissenschaft werden 
besprochen und ein Ausblick auf die technische Entwicklung und die méglichen 
Fortschritte in der nachsten Zukunft gegeben (wobei solchen Extrapolationen 
natiirlich immer ein grosses Mass an Unsicherheit eigen ist). Ebenso werden weiter- 
gehende Anwendungen, wie Ubersetzungen, Dechiffrage, Dokumentation, kurz er- 
wahnt. Die Darstellung hat vor allem die Verhaltnisse und die Modglichkeiten 
Frankreichs im Auge und wendet sich wohl in erster Linie an leitende Persénlich- 
keiten industrieller und kommerzieller Unternehmungen, fiir die gelegentlich auch 
praktische Hinweise verschiedener Art gegeben werden. Einige bibliographische 
Angaben beschliessen das Buch. E. RotH-DESMEULES 


Theorie der Limitierungsverfahren. Von Kari ZELLER ([Ergebnisse der 
Mathematik und ihrer Grenzgebiete. Neue Folge, Heft 15] Springer-Verlag, Berlin 
1958). 242 S.; DM 36.80. 

Die Theorie der Limitierungsverfahren, speziell fiir divergente Zahlenfolgen, ist 
zu einem umfangreichen Gebiet der Analysis angewachsen, so dass eine knappe 
und iibersichtliche Zusammenfassung der Ergebnisse sehr erwiinscht ist, auch wenn 
heute noch keine einheitliche und abschliessende Darstellung méglich ist. Das Buch 
von ZELLER kommt diesem Wunsche in hervorragender Weise entgegen und ist 
vor allem geeignet, den Leser in die Hauptprobleme sowie in die vorhandene 
Literatur einzufiihren. Zur Hauptsache werden die Matrixtransformationen und 
die zugehérigen Limitierungsverfahren betrachtet, wahrend andere Methoden und 
Anwendungen nur beilaufig erwahnt werden. 

Die ersten Kapitel des Buches bringen die Grundbegriffe, die Hilfsmittel aus 
der Funktionalanalysis, die Struktur der Wirkfelder und weiter, im Rahmen der 
allgemeinen Theorie, die direkten Satze und die Umkehrsatze. Die folgenden Kapi- 
tel befassen sich mit speziellen Verfahren, wie die vom Cesaro-Abelschen Typ, die 
funktionentheoretischen Verfahren sowie andere Klassen von Verfahren. Die Be- 
weise sind im allgemeinen nur angedeutet. Von besonderer Wichtigkeit ist das bei- 


_gefiigte Literaturverzeichnis von nicht weniger als 64 Seiten, das bis 1955 nachge- 


fiihrt ist und kaum eine bedeutende Arbeit vermissen lassen diirfte. Das Buch 


- wendet sich in erster Linie an den Mathematiker. E. RotH-DESMEULES 


ray =e 
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An Introduction to Fourier Methods and the Laplace Transformation. 
Von Puiip FRANKLIN (Dover Publications, London and New York 1958). 289 S., 
TAM SATO. 5 ; 

Diese Einfiihrung ist eine korrigierte Neuauflage des bereits 1949 erschienenen 
Buches Fourier Methods und befasst sich mit den Fourier-Reihen und der Laplace- 
Transformation, und zwar vor allem im Hinblick auf ihre Anwendungen auf phy- 
sikalische Probleme. Nach einem einleitenden Kapitel iiber komplexe Zahlen und 
elektrische Schwingungskreise werden die Fourier-Reihen und Fourier-Integrale 
behandelt. Alsdann folgt eine Einfiihrung in die partiellen Differentialgleichungen 
(insbesondere Warmeleitung, Schwingungsprobleme, Maxwellsche Gleichungen 
usw.) und die zugehérigen Randwertprobleme. Schliesslich wird die Laplace- 
Transformation eingefiihrt mit Anwendungen auf gewohnliche Differentialglei- 
chungen sowie auf gewisse partielle Differentialgleichungen. An Vorkenntnissen 
wird im wesentlichen nur die Differential- und Integralrechnung vorausgesetzt. Die 
Darstellung ist sorgfaltig und enthalt iiber 800 Aufgaben mit den Losungen sowie 
einige bibliographische Hinweise. Das Buch ist gut geeignet fiir Studierende, wo- 
bei auch der recht giinstige Preis nicht ausser acht gelassen werden darf. 

E. RotH-DESMEULES 


Elementary Statistical Physics. Von C. Kitrer (John Wiley & Sons, 
New York 1958). 228 S., 28 Fig.; $ 8.-. 

Der Biichermarkt wurde in der letzten Zeit mit Neuerscheinungen tiber stati- 
stische Mechanik férmlich tiberschwemmt. Eine gewisse Aktualitat ist dem Gebiet 
wegen den zahlreichen neueren Anwendungen vielleicht nicht abzusprechen. Davon _ 
ist in dem Buch von KITTEL erfreulicherweise manches zu finden; wir erwahnen 
stichwortartig: Negative absolute Temperatur, Streuung extrem hoch energeti- 
scher Teilchen, Strahlungsrekombination von Elektron-Loch-Paaren in Halb- 
leitern. Leider hat sich der Leser dabei meist mit Andeutungen und Hinweisen auf 
die Originalliteratur zu begniigen, denn das Buch will «elementary sein. Aber 
KitteL hat die Gabe, schwierige Fragen auf einfache zu reduzieren (siehe zum 
Beispiel die Behandlung des Overhauser-Effektes als Ubungsaufgabe). Das Buch 
wendet sich an die Physikstudenten und will eine Einfithrung vermitteln vor allem 
in die Methode der Ensembles, die stochastischen Prozesse und die Transport- 
theorie. C. ENZ 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLtcce, Band 26, 2. Teil: Licht und Materie (Springer-Verlag, Berlin 1958). 
965 S., 754 Abb.; DM 168.-. 

Der 26. Band enthalt 4 Kapitel sehr verschiedenen Umfanges und Stils iiber 
Lumineszenz (G. F. J. Garuick, 128 S.), Temperaturstrahlung fester Kérper (G. A.N. 
Rutcers, 42S.), Raman-Effekt (S. MizusHima, 72 S.) und Infrarotspektvoskoptie 
(J. LecomrTeE, 695 S.). 

( Der Artikel tiber Lumineszenz ist fiir Festkérperphysiker und (zum Teil auch 
fiir physikalische) Chemiker geschrieben und geht von den alten einfachsten 
Theorien der Kinetik (Reaktion 1. und. 2. Ordnung fiir die Lumineszenzzentren, 
JABLONSKIS Termschema mit metastabilem Zwischenzustand und thermischem 
Aktivierungsprozess) zu den heute tiblichen Vorstellungen der Energiebander (Zu- 
stande der Lumineszenzzentren und der «electron traps» im «gap» zwischen Va- 
lenz- und Leistungsband) iiber. Das Kapitel iiber Beobachtungen an verschiedenen 
Phosphor- Klassen (Frank-Condon-Prinzip) wird vom Standpunkt dieser Theorien 
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und des Hippelschen Modells vorgenommen, die Struktur der Lumineszenzzentren 
wird dabei nicht eingehend behandelt. Bei der Bésprechung einzelner Phosphore 
wird die Theorie verfeinert, speziell fiir Thermolumineszenz bei metastabilen 
Zentren, bei photohalbleitenden Phosphoren, mit und ohne strahlungslose Pro- 
zesse und bei der wichtigen Frage der Polarisation der Lumineszenzstrahlung. 
Eine relativ ausfiihrliche Behandlung findet ferner die Elektro- und die praktisch 
wichtige Radiolumineszenz. 

Der Rutgersche Artikel iiber Temperaturstrahlung scheint unseres Erachtens 
hauptsachlich fiir Ingenieure geschrieben, indem er sich auf die Beschreibung der 
Temperaturstrahlung anhand phaénomenologischer Gesetze beschrankt und eine 
Reihe von Messungsergebnissen fiir Metalle und technisch interessante Oxyde und 
dergleichen zitiert. 

Der Artikel iiber den Raman-Effekt wiederum scheint hauptsachlich fiir den 
Chemiker gedacht, die klassische Theorie des Effektes, ihr Zusammenhang mit der 
Symmetrie der Molekel und deren Grundschwingungen sowie ihre Stellung im 
Rahmen der Strahlungstheorie ist angedeutet, ebenso die (gruppentheoretische) 
Ableitung der Auswahlregeln des Schwingungs-Raman-Effektes. Als Beispiele sind 
einige experimentelle Resultate angegeben. Man kann sich fragen, warum die 
Herausgeber nicht einfach eine der bewahrten Darstellungen des Raman-Effektes 
abgedruckt haben. 

Obgleich nach Angabe des Autors hauptsachlich fiir Physiker gedacht, diirfte 
der als monumental anzusprechende Artikel iiber Infrarotspektroskopie am besten 
den Bediirfnissen des Chemikers dienen. Die Theorie der Infrarotspektren ist fast 
vollig beiseite gelassen. Dafiir findet man zundchst eine (im wesentlichen qualita- 
tive) Besprechung der experimentellen Hilfsmittel (Spektrographen, Materialien 
usw.). Ferner ist in umfangreichen Kapiteln mit rund 500 (in vielen Fallen unndtig 
uneinheitlichen und gelegentlich unvollstandig beschrifteten) Figuren eine bis 1956 
reichende Sammlung experimenteller Ergebnisse iiber Infrarotspektren organischer 
und anorganischer Molekeln und Kristalle zusammengestellt. Einzelne Kapitel 
zeigen eine bisher nicht bekannte Vollstandigkeit hinsichtlich der gegenwartigen 
empirischen Kenntnisse und Literaturzitate, sind aber weitgehend qualitativ ge- 
halten und auch nicht durchweg fehlerfrei. Speziell die Diskussion der Kristall- 
spektren ist unseres Erachtens zu vereinfacht. Ein Kapitel iiber die bisherigen Er- 
gebnisse der Intensitatsmessung der Infrarot-Interferometrie und die photometri- 
sche qualitative Analyse schliesst den Artikel ab. Hs. H. GUNTHARD 


Technische Hydro- und Aeromechanik. Von W. KaurMann (Springer- 
Verlag, Berlin 1958). 386 S., 266 Abb.; DM 37.50. 
Diese zweite, verbesserte Neuauflage stellt fiir den jungen Ingenieur und Studen- 
ten auf dem Gebiet der Stromungsmechanik ein dusserst wertvolles und praktisches 
Lehrbuch dar. Die rasch vergriffene erste Auflage bestatigte die ausserordentlich 
giinstige Wahl der behandelten Gebiete und bedingte nur wenige, eher formale 
Abianderungen des alten Textes. Zugleich sind aber dartiber hinaus einige Ergan- 
_zungen neu hinzugenommen worden. el 
In seinen drei Hauptabschnitten, Eigenschaften der Fliissigkeiten und Gase, 
“Hydro- bzw. Aerostatik und der Lehre von den Bewegungen der Flissigkeiten 
(Dynamik), ist das Hauptgewicht auf das stoffbedingt grosste Gebiet der Dynamik 
gelegt. Das Werk zeichnet sich aus durch eine klare Beschreibung der physikali- 
-schen Grundlagen und folgert die Theorien in enger Anlehnung an die Praxis. 
Viele Gebiete wurden durch neueste Forschungsergebnisse bereichert, ohne dass da- 
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bei der Charakter eines Grundlagenwerkes verlorenging. Mancherorts sind zum 
tieferen Verstandnis die Erlauterungen durch Rechenbeispiele erganzt. Es wiirde 

hier zu weit fiihren, die behandelten Gebiete im einzelnen zu besprechen, es sei 

lediglich hervorgehoben, dass sémtliche dem heutigen Stand der Technik der 

Strémungslehre entsprechenden Grundgebiete behandelt werden. Angefangen bei 

den Bewegungsgleichungen von EvuLer, dem Energiesatz fiir inkompressible 

Fliissigkeiten von BERNOULLI bis zum Charakteristikenverfahren von PRANDTL 

und BuSEMANN wird iiberall, durch wertvolle Hinweise auf die einschlagige Litera- 

tur, der Weg zur vertiefteren Ausbildung erleichtert. 

Die oben erwahnten Erginzungen betreffen vor allem die Arbeiten iiber die 
Grenzschichttheorie, Tragfliigeltheorie zur Berechnung des Auftriebes an gepfeilten 
und schiebenden Fliigeln. Vermehrte Aufmerksamkeit wurde dem Verhalten der 
Unstetigkeitsflache hinter Tragfliigeln geschenkt. Die Gasdynamik hat auch eine 
der Aktualitat entsprechende Erweiterung erfahren. Das im Gas- und Dampf-_ 
turbinenbau wichtige Verhalten der Grenzschicht an rotierenden Scheiben hat in 
einem neuen Abschnitt seine Wiirdigung gefunden. Das Buch kann in dieser neuen 
Fassung wiederum bestens empfohlen werden. V. STINGELIN © 


Numerische Mathematik. Zeitschrift herausgegeben von R. SAUER (Min- 
chen), E. StreFeL (Ziirich), J. Topp (Pasadena) und A. WALTHER (Darmstadt) 
bei Springer (Heidelberg). Bd. 1, Heft 1 (1959), 60 Seiten, Gr. 8°; brosch., 
DM 16.80. 

Der Springer-Verlag bietet den an der numerischen Analysis interessierten 
wissenschaftlichen Kreisen eine neue Zeitschrift an («Numerische Mathematik»), 
deren 1. Heft nun in vorziiglicher Ausstattung vorliegt. Die Zeitschrift soll auf 
breiter internationaler Grundlage Arbeiten verdffentlichen, die sich mit allge- 
meinen Problemen des digitalen Rechnens und der Entwicklung von numerischen 
Verfahren beschaftigen. . 

Das wissenschaftliche Redaktions-Komitee besteht aus den Herren R. SAUER 
(Miinchen), E. STIEFEL (Ziirich), J. Topp (Pasadena) und A. WALTHER (Darm- 
stadt). Das erste Heft enthalt Aufsatze von E. STIEFEL, A. S. HoUSEHOLDER und 
F. L. Baver, M. WEISFELD, sowie einen interessanten Bericht von A. J. PERLIS- 
und K. SAMELSON iiber die neue Formelsprache aLGor, die dazu dienen soll, die 
Methoden der automatischen Programmierung fiir Rechenautomaten auf inter- 
nationaler Grundlage einheitlich zu entwickeln. H. RUTISHAUSER 


Kerntechnik, Physik, Technologie, Reaktoren. Von W. RizzLEeR und 
W. WatcHeR (B. G. Teubner, Verlagsgesellschaft Stuttgart 1958). 1002 S., 
493 Fig., 154 Tab.; DM 150.-. ; 

Gemessen am Werdegang der uns bereits vertrauten Wissenschaften, hat sich 
die Reaktortechnik dusserst sprunghaft entwickelt. Kaum zwei Jahrzehnte nach 
der erstmaligen Verwirklichung einer sich selbst unterhaltenden Kettenreaktion 
mit Kernbrennstoff sind Nutzenergie erzeugende Kernkraftwerke in Betrieb, und 
die Errichtung weit grésserer Einheiten steht in naher Zukunft bevor. Diese stir: 
mische Entwicklung ruft gebieterisch nach Fachleuten aller Ausbildungsstufen. 
Hierbei ist besonders die dem Akademiker zugedachte Aufgabe eine mannigfal- 
tige, verlangt doch der Reaktorbau Kenntnisse der mathematischen, physikali- 
schen, chemischen, metallurgischen und technologischen Wissenschaften. Es ist 
gerade diese Vielseitigkeit, welche einerseits diesem modernen Betatigungsfeld 
einen so sonderlichen Anreiz gibt, welche aber auch dem Entstehen einer Lehr- 
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tradition, wie sie zum Beispiel fiir die Thermodynamik vorhanden ist, Hindernisse 
in den Weg legt. Die moderne Literatur itiber Reaktortechnik kénnte man daher 
in zwei Gruppen scheiden. Die erste umfasst Werke, welche die Anwendungen der 
einzelnen Wissenschaften fiir die Probleme der Reaktortechnik zum Thema 
haben. Die andere Gruppe enthalt Sammlungen der verschiedensten Gesichts- 
punkte, welche beim Reaktorbau beriicksichtigt werden miissen. Je nach Vorbil- 
dung und Veranlagung der Autoren werden dabei die einen oder anderen Probleme 
griindlicher, klarer oder mit mehr pidagogischem Geschick behandelt. Solche 
Biicher sind naturgemass oft sehr umfanglich, wie zum Beispiel das vorliegende 
erste deutschsprachige Werk Kerntechnik von RirzLER und WALCHER. 

Die Herausgeber haben sich der Mitarbeit einer grossen Zahl von Fachleuten 
versichert, so dass die einzelnen Kapitel mit der nétigen Wissenschaftlichkeit 
abgefasst werden konnten. Ein solches Vorgehen bringt zwar notwendigerweise 
eine zusatzliche Heterogenitat mit sich, indem die verschiedenen Autoren an ihre 
Leser unterschiedliche Anforderungen stellen. So erfordert unter anderem der 
Abschnitt tiber Reaktortheorie eine mathematische Vorbildung, iiber welche der 
heute ausgebildete Ingenieur nicht ohne weiteres verfiigt. Andererseits hatte die 
Behandlung der reaktorkinetischen Vorgange durch Verwendung der beim Regel- 
fachmann als bekannt anzunehmenden Methoden der Laplace-Transformation an 
Anschaulichkeit gewonnen. Obwohl die Herausgeber sich um eine einheitliche 
Bezeichnungsweise bemiihen, wird beispielsweise der Ausdruck Reaktivitat von 
den Teilautoren in unterschiedlicher Weise gebraucht (vgl. Definition S. 500 mit 
Bedeutung im Kapitel Wdrmeprobleme, zum Beispiel S. 598), 100 Seiten des 
Buches dienen der Beschreibung von Forschungs- und Leistungsreaktoren, wobei 
anhand ausgewahlter Beispiele die besonderen Eigenschaften und Aufbauprin- 
zipien der verschiedenen Systeme erlautert werden. Uber Sicherheitsfragen, die 
mé6glichen Griinde fiir ein Abweichen eines Systems von vorgeschriebenen Be- 
triebszustanden und die notwendigen organisatorischen Massnahmen sowie Richt- 
linien zur Standortwahl des Reaktors handelt ein weiterer Abschnitt. Diese 
Probleme werden anhand der Erfahrungen einiger Reaktorunfalle diskutiert. Die 
einleitenden Abschnitte des Buches vermitteln eine Einfiihrung in die kernphysi- 
kalischen sowie die zum Versténdnis des Reaktorbaues notigen technologischen 
und verfahrenstechnischen Grundlagen. In die Zukunft weisen Bemerkungen 
iiber Kernfusion. 

Der ganze Band enthalt eine ungewoéhnlich grosse Zahl von Figuren und 
Tabellen. Ein ausfiihrliches Sachverzeichnis erleichtert den Gebrauch als Nach-’ 
schlagewerk wesentlich. 

Die Autoren haben sich bemiiht, die durch die angelsachsische Originallite- 
ratur eingefiihrten Bezeichnungen ins Deutsche zu tibertragen. Dieser Versuch 
ist ihnen weitgehend gelungen. Einige Ausdriicke, wie Stromkammer fiir [onisa- 

-tionskammer oder Wanderflache fiir Migrationsflache, werden jedoch kaum 
_allgemeine Zustimmung finden. W. HALG 


Tabellen der Atomkerne. Teil I: Eigenschaften der Atomkerne. Band 1: 
_ Die Elemente Neutron bis Zinn. Von W. Kunz und J. SCHINTLMBEISTER (Akademie- 
Verlag, Berlin 1958). 465 S., 305 Abb.; DM 105.-. © 
Bei dem vorliegenden Band handelt es sich um ein Tabellenwerk, das die 
-charakteristischen Daten der stabilen und unstabilen Isotope der Elemente 
Neutron bis Zinn erfasst. Den Tabellen geht eine in deutscher, englischer und 
_russischer Sprache abgefasste Erlauterung voraus, in der die in einem Tabellen- 
werk unumganglichen Abkiirzungen sorgfaltig interpretiert werden. Begriffe wie 
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zum Beispiel Konversionskoeffizienten, die in der Literatur zuweilen verschieden 

definiert erscheinen, werden eindeutig festgelegt. — Samtliche tiber einen Kern bis 

zum 1, Januar 1958 bekannten Daten wurden kritisch durchgesehen und gemass 

ihrer Zuverldssigkeit klassifiziert. An Stellen, wo dies zulassig erscheint, haben 
die Autoren aus vorliegenden Messresultaten einen sorgfaltig gewogenen Mittel- 

wert berechnet. Die Zahlenwerte werden, soweit dies méglich ist, mit Fehler- 

grenzen, Messmethode und Autor angegeben. 

Die stabilen Isotope sind in diesem Werk durch folgende Grossen charakteri-_ 
siert: Genaue Atommasse, relative Haufigkeit des Isotops, Spin, magnetisches 
Dipol- und Oktupolmoment sowie elektrisches Quadrupolmoment. Isomere Zu- 
stande mit Halbwertszeiten, die 0,1 us iiberschreiten, erscheinen separiert. Fiir die 
radioaktiven Kerne ist das Zerfallsschema — sofern es mit einiger Sicherheit an-_ 
gegeben werden kann — graphisch dargestellt. Der Zerfall wird tiberdies durch 
folgende tabellierte Werte beschrieben: Halbwertszeit, Energie der Teilchen und 
Quanten, Spin und Paritat der auftretenden Zustande sowie Konversionskoeffi- 
zienten und Multipolordnung der Gammaiibergange. Die prozentuale Intensitat 
der angegebenen Strahlungskomponenten, die ft-Werte fiir Betazerfalle und die 
méglichen Methoden zur Herstellung des betreffenden Isotopes vervollstandigen 
die Angaben. 

Bei der Ausseren Gestaltung des Buches wurde weniger Wert auf eine gedrangte 
Darstellung gelegt; dafiir gewinnen die Tabellen an Ubersichtlichkeit. Mess- 
groéssen, die zurzeit noch nicht zur Verfiigung stehen, konnen leicht in die ent- 
sprechenden Kolonnen eingetragen werden. Eine merkliche Erleichterung fiir den 
Bentitzer der Tabellen bedeutet es, dass jedem Element ein separates Literatur- 
verzeichnis zugeordnet ist. Auf diese Weise wird das so lastige Nachschlagen von 
Zitaten vermieden. 

Der dauerhafte Einband und die vorziigliche Ausstattung des Buches lassen 
dieses Werk als wertvolles Hilfsmittel nicht nur des Kernphysikers und Radio- 
chemikers, sondern auch all derjenigen erscheinen, die sich mit der Anwendung 
von Isotopen in Medizin, Biologie und Technik befassen. F. HEINRICH 


Introduction to Advanced Dynamics. Von S. W. Mc CusKry (Addison 
Wesley Publ. Co., Inc., Reading, Mass., 1959). 263 S., 120 Fig.; $ 8.50. 

Eine hibsche und saubere Einfiihrung in die héhere Dynamik. Sie behandelt 
ike die Bewegung (und Relativbewegung) des Massenpunktes, 2. das Hamiltonsche 
Prinzip- und die Lagrangeschen Gleichungen, 3. die Zentralbewegung, 4. die 

Dynamik des starren K6rpers mit einer Einfiihrung in die Probleme des Kreisels, 
5. Schwingungen von konservativen und nicht konservativen Systemen mit zwei 
Freiheitsgraden, 6. die Hamiltonschen Gleichungen und den Phasenraum, 7. die 
Hamilton-Jacobischen Gleichungen. Das Buch wendet sich in erster Linie an 
Studenten der Physik, ist aber 4usserst klar und einfach geschrieben und damit 


fiir jeden verstandlich, der einen elementaren Mechanikkurs der bei uns iiblichen 
Art hinter sich hat. H. ZIEGLER 


Statics and the Dynamics of a Particle. Von W. Duncan Mac MILLAN 
(Dover Publications, London and New York 1958). 450 S.,.183 Fig.; $'2,00, 

Das Buch gibt eine breite und elementare Darstellung der Punktmechanik, 
erganzt durch einzelne Kapitel der Statik starrer Korper und der Festigkeits- 


lehre. | 
H. ZIEGLER 
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NEUERSCHEINUNGEN 


Petrochemische Berechnungsmethoden 


auf aquivalenter Grundlage 
a: (Methoden von P. Nice) 


Von Professor Dr. Conrap Burrt 
Professor an der Eidgendssischen Technischen Hochschule und an der Universitat Ziirich 


334 Seiten mit 65 Figuren (1959). Preis gebunden Fr./DM38.—. Mineralogisch-geotechnische 
Reihe, Band 7. palning «Lehrbiicher und Monographien aus:dem Gebiete der exakten 
Wissenschaften». 


Inhalt: Einleitung — Methoden der graphischen Darstellung — Hilfsmittel — Die Niggli- 
Werte — Beziehungen zu Gewichts- und Kationenprozenten — Anwendung zur Veranschau- 
lichung fundamentaler petrochemischer Bezichungen — Grenzen der Anwendung — Die 
Aquivalentnorm — Prinzipielles — Basis und Standard-Katanorm — Beziehung zu Gewichts- 
und Aquivalentprozenten sowie Niggli-Werten — Anwendung zum Studium der Hetero- 
morphie und der metamorphen Gesteine — Beriicksichtigung der Anionen nach Bartu und 
Esxoia — Stoffbilanzen isovolumetrischer metasomatischer Prozesse — Tabellen der 
1000fachen atomaren und molekularen Aquivalentzahlen. 


ELEKTRIZITATSLEHRE 


Von HeEtnz ScHILT 
Privatdozent an der Universitat Bern 


Mit einem Geleitwort von W. Pautt. 216 Seiten mit 187 Figuren (1959). In Ganzleinen 
Fr./DM 24.—. Physikalische Reihe, Band 3.. Sammlung «Lehrbticher und Monographien 
aus dem Gebiete der exakten Wissenschaften». 


Aus dem Inhalt: 1. Grundbegriffe und Elektrostatik. 2. Das stationdre elektrische Feld. 
8. Das stationdire magnetische Feld und langsam veranderliche Felder. 4. Rasch ver- 
, anderliche Felder. 5. Anhang. ~ 


| Es besteht heute noch ein gewisser Gegensatz zwischen einer elementaren Elektrizitats- 

lehre und der wissenschaftlich mustergiiltigen Maxwellschen Theorie der Elektrizitat. Das 

- Buch yon Scuirr fiihrt den Leser mit einfachen mathematischen Methoden in die Maxwell- 

sche Auffassung der Elektrizitat ein und stellt anhand vieler Beispiele den Zusammenhang 

mit der Praxis her. Es gibt der elementaren Lehre neue Impulse, die zum Uberwinden des 

tac ie Gegensatzes notwendig sind, und hilft damit dem Anfanger iiber unn6dtige 
Schwierigkeiten hinweg. 
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